Departamento de Fisica Teorica I 
Universidad Complutense de Madrid 



Doctoral Thesis 

Bipartite entanglement of 
localized separated systems 



Author: Supervisor: 
Carlos Sab in Lestayo Dr. Juan Leon Garcia 

Instituto de Fisica Fundamental 
Consejo Superior de Investigaciones Cienti'ficas 




November 2011. Committee: M. A. Martm-Delgado, A. Luis, I. Cirac, E. 
Solano, A. Cabello. 



Para Mencha. 
Para Nicolas. 



i 



Agradecimientos 



"People are afraid to face how great a part of life is dependent on 
luck" (Woody Allen, Match Point) 

Dcjenme que les cuente una historia sobre el azar. 

El azar llevo a un estudiante sin ambiciones investigadoras de los ultimos 
cursos de Fi'sica a una conferencia divulgativa de Guillermo Garcia Alcaine sobre 
Mecanica Cuantica, dirigida a los alumnos. Al final de la charla, Guillermo hizo 
publicidad sobre un Trabajo Academicamente Dirigido que le gustaria dirigir, 
con un tenia relacionado con el entrelazamiento. Puesto que mis companeros es- 
taban en clase (dondc sin duda deben'a estar yo tambicn), practicamente nadie 
de los ultimos cursos recibio esa propaganda, y fuc la unica que se hizo. Asf 
que no tuve ninguna competencia cuando al dia siguicntc apareci por su despa- 
cho. Ese fue el comienzo de mas de dos ahos de investigation, que dieron lugar, 
ademas del mencionado trabajo, a una Tesis de Master y un arti'culo sobre clasi- 
ficacion y medida del entrelazamiento. El recuerdo de Guillermo y su sentido 
del humor me acompaha siempre, con su ejemplo de rigor, honestidad y perfec- 
cionismo, cualidades que en mi se dan en tan pequehas dosis que son solo torpes 
rcmcdos de las suyas. Solo tengo para el admiration y gratitud. 

Por motivos que no vicnen al caso, Guillermo ya no dirigia Tesis Doctorales. 
Por el camino yo me habi'a ido creyendo que podia dedicarme a eso de la inves- 
tigation, pero el ticmpo iba pasando sin que se me abriera ninguna puerta. El 
dia que hice el ultimo examen de la carrera, mi costumbre de leer los cartclcs 
de las paredes (que tanto solia exasperar a mis acompahantes) me llevo a leer 
un anuncio de Tesis Doctoral de un tal Juan Leon del Instituto de bla, bla, bla 
sobre... jteoria del entrelazamiento! ^Ven lo que les deci'a sobre el azar? 

La expresion "padre cientifico" se usa muchas veces con cierta alegrfa para 
designar la relation de cierto investigador con su director de Tesis. Pocas veces 
es tan cierta como en mi caso. Juan confio en mi sin tener demasiados motivos 
para hacerlo y decidio hacer frente, por pura bondad, a todas las dificultades 
que desde el principio se vio que yo le procuraria. Tuvo varias ocasiones para 
deshacerse elegantemente de mi, pero no lo hizo. En lugar de eso se desvivio 
por ayudarme. Me enseno rigor, disciplina intelectual, ambition investigadora, 
y despues comprendio que podia, que debia dejarme volar con libertad. Hoy 
yo le hablo ya como los hijos adolescentes hablan a sus padres, mezclando la 
indiferencia con el sarcasmo. Es pura mascara: jamas olvidare todo lo que Juan 



ii 



ha hecho por mi. 

Guillcrmo, Juan... a hombros de gigantes tales, ^quien no acertaria a ver 
con claridad? Pero recuerden que yo querfa hablarles del azar. Yo ya llevaba 
un tiempo en el Consejo, aguantando cada manana la mirada cntre ironica ("£a 
donde va este?") y admonitoria ("jhaz algo, idiota!") de D. Santiago Ramon y 
Cajal y D. Severo Ochoa, (cuyas estatuas se erigieron a los pocos di'as de llegar 
yo en sustitucion del busto de algun olvidado politico de la dictadura), cuando 
Juanjo Garcia Ripoll saco plaza en nuestro grupo. Un dfa que comimos los dos 
juntos en el comedor, surgio la idea de intentar ver que aplicaciones podn'an 
terner los resultados que habfamos ido obtcnicndo Juan y yo en el campo emer- 
gente de los circuitos superconductores. Ahi empezo nuestra colaboracion con 
Juanjo. Su asombrosa intuition, cimentada sobre un abrumador conocimiento 
de toda la Fi'sica experimental y teorica de los siglos XX y XXI, permitieron 
dar a esta Tesis otro salto de calidad. Gracias a el ademas pude entrar en 
contacto con Kike (asf quiere el que lo escribamos) Solano, cuyo entusiasmo 
y atrevimiento, su vision poetica y romantica del oficio, son un cstfmulo para 
todos los investigadores que nos cruzamos en su camino, quienes secretamente 
aspiramos a ser, cuando seamos jovenes, como el. 

Pero basta ya de hacer la pelota a los maestros. Durante la elaboration de 
esta Tesis, como es natural, conoci a muchas personas. La mejor de todas ellas 
es mi buen amigo Borja Peropadre. Su bondad, sensatez, perspicacia y sabidurfa 
me hacen seguir confiando en la condition humana y tratar de ser mejor. jQuicn 
le iba a decir a un sociopata como yo que iba a disfrutar tanto compartiendo 
un dcspacho! Tambicn son buenos amigos cl inimitable Marco del Rey, cuya 
exigencia intelectual y espiritu critico elevaron considerablemente la calidad del 
capitulo 5 de esta Tesis, el Dr. Edu Martin Marthez, el SuperManager Emilio 
Alba, Andrea Cadarso y el resto de gente de QUINFOG, asi como Alejan- 
dro Bcrmudez, Juan Manuel Perez Pardo, Fermando Lledo, Luis Garay, Lucas 
Lamata, Pol Forn, Jorge Casanova y quien sabe cuanta gente que olvido, pero 
con la que habre compartido al menos el cafe aguado de algiin congreso. Dire 
tambien que sin Paco y sin David, y el resto del Centro de Estudios Cuzco, no 
habria podido terminar esta Tesis. 

Muchos padres intentan educar a sus hijos a su imagen y semejanza, con- 
fiando en que repitan sus doctrinas y consignas. Mis padres no me ensenaron 
otra ideologfa que el amor y los buenos sentimientos, y les estoy muy agrade- 
cido por ello. Si hubiera heredado la mitad de la originalidad y libertad de 
pensamiento de mi padre Jose Ramon, de su inquietud y su curiosidad, de 
su insobornable integridad c individualidad, de su amor por la lectura y el 
conocimiento, de su profundidad para entender lo esencial, me daria por satis- 
fecho. Si algo hubiera aprendido del amor incondicional de mi madre Rosalia, 
de su entrega desinteresada por la familia, de su capacidad para seguir adelante 
por amargo que sea lo que le ofrezca la vida, para tirar de recursos y mantenerlo 
todo en pie y en orden con originalidad y creatividad, ya serfa mucho. Espcro 
que cuando tengan esta Tesis en sus manos sientan, aunque sea por un instante, 
que en el fondo merecio la pena. 

Todo lo que soy, para bien o para mal, tiene que ver con mi hcrmano Mon, 



in 



pocta, profcsor, ingcnicro y filologo (si mcncionamos solo lo mas evidentc). El 
me enseno, entre otras muchas cosas, a ser libre, a amar la cultura, a no tener 
miedo a pensar y a seguir tu propio camino, y aunquc 61 no se acuerde, me 
convencio sin pretenderlo para que estudiara Fi'sica en lugar de Filosoffa. Su 
influencia sobre mi ha sido tan grande que lo he tenido que asesinar en sen- 
tido Freudiano muchas veces, para resucitarlo inmediatamente. Y mi hermana 
Maria, a la que llamaba cuando lloraba desconsolado dc nino, completo mi for- 
mation como hombre, con su profundo sentido de la generosidad y la justicia. 
No puedo evitar acordarme ahora de mi tfa Carolita, que ya no podra felicitarme 
desde el cielo. 

Durante la elaboration de esta Tesis me case, mi mujer tuvo dos cmbarazos, 
nacio un hijo al que he visto crecer, aprender a andar y a construir sus primeras 
frases. Esta Tesis es para ellos. Para Menchita, que ha confiado en mi capacidad 
y en mis posibilidades mucho mas que yo mismo, que ha sacrificado tanto para 
que llegara este dfa, que con su sola presencia da sentido a la extrana comedia 
de estar vivo, de poner un pie delante de otro y caminar, uno, dos, inspirar, 
espirar, para Menchita que es el unico motivo por el que he seguido en pie tras 
esa noches en las que despertaba de pronto y sentfa tanto miedo. Sin ella, "las 
estrellas, a pesar de su lampara encendida pcrdcrian el camino/ ^que seria del 
Universo?" Para Nicolas, que pronto cumplira dos ahos, a quien deseo que 
no entienda jamas ni una palabra de lo que dice esta Tesis, que encuentre la 
felicidad mas facilmcnte que su viejo, su pobre viejo que no podra enseharle mas 
que algun poema y aquello que dice esa vieja cancion de los 60, la misma que 
escuchaba Desmond Hume en el bunker: 

"You have to make your own kind of music/ sing your own special 
song/ make your own kind of music/ even if nobody else sings along" 

Madrid, Septiembre del 2011. 



iv 



List of publications 



Chapter 3 

• (In this chapter, authors in alphabetical order) 

• "Entanglement swapping between spacelike separated atoms" , Juan Leon, 
Carlos Sabi'n, Phys. Rev. A 78, 052314 (2008). 

• "Photon exchange and correlation transfer in atom-atom entanglement 
dynamics", Juan Leon, Carlos Sabi'n, Phys. Rev. A 79, 012301 (2009). 

• "Generation of atom-atom correlations inside and outside the mutual light 
cone", Juan Leon, Carlos Sabi'n, Phys. Rev. A 79, 012304 (2009). 

• "A case of entanglement generation between causally disconnected atoms" , 
Juan Leon, Carlos Sabi'n, Int. J. Quant. Inf. 7, 187-193 (2009). 

• "Atom-atom entanglement generated at early times by two photon emis- 
sion", Juan Leon, Carlos Sabi'n, Phys. Scr. T135 014034 (2009). 

Chapter 4 

• "Dynamics of entanglement via propagating microwave photons" , Carlos 
Sabm, Juan Jose Garcfa-Ripoll, Enrique Solano, Juan Leon, Phys. Rev. 
B 81, 184501 (2010). 

• "Detecting ground-state qubit self-excitations in circuit QED: A slow 
quantum anti-Zeno effect", Carlos Sabm, Juan Leon, Juan Jose Garcfa- 
Ripoll, Phys. Rev B 84, 024516 (2011). 

Chapter 5 

• "The Fermi problem with artificial atoms in circuit QED" , Carlos Sabm, 
Marco del Rey, Juan Jose Garcfa-Ripoll, Juan Leon. Phys. Rev. Lett. 
107, 150402 (2011). 



vi 



• "Short-time quantum detection: probing quantum fluctuations" , Marco 
del Rey, Carlos Sabm, Juan Leon. Submitted to Phys. Rev. Lett. arXiv: 
1108.0672. 

Chapter 6 

• "Quantum simulation of the Majorana equation and unphysical opera- 
tions", Jorge Casanova, Carlos Sabi'n, Juan Leon, liiigo L. Egusquiza, 
Rene Gerritsma, Christian F. Roos, Juan Jose Garcfa-Ripoll, Enrique 
Solano. Accepted to Phys. Rev. X. |arXiv:1102.1651| 

• To be submitted "On Majorana Hamiltonians" , Carlos Sabm, Jorge 
Casanova, Juan Jose Garcfa-Ripoll, Juan Leon, Enrique Solano, Ihigo L. 
Egusquiza. 

• To be submitted "Quantum simulation of relativistic potentials without 
potentials" , Carlos Sabm, Jorge Casanova, Juan Jose Garcfa-Ripoll, Lucas 
Lamata , Enrique Solano, Juan Leon. 

Other publications not included in this Thesis 

• "A classification of entanglement in three-qubit systems" , Carlos Sabm, 
Guillermo Garcia- Alcaine, European Physical Journal D, 48, 435-442 (2008) 



vii 



Contents 



Agradecimientos [I] 

List of publications [v] 

1 Introduction [TJ 

2 Preliminaries [8j 

2.1 Matter-radiation interaction. Hamiltonian [3J 

2.1.1 Introduction [5J 

2.1.2 Two level-atoms. Interaction picture [TT] 

2.1.3 Rotating wave approximation and beyond [TTJ 

2.1.4 Rabi and Jaynes-Cummings models [12] 

2.2 Circuit QED [JJ 

2.3 Entanglement measures in two-qubit systems [T3J 

3 Dynamics of entanglement in matter-radiation interaction 1161 

3.1 Generation of entanglement inside and outside the mutual light 
cone [TO] 

3.1.1 Introduction [TC] 

3.1.2 The model [TJ 

3.1.3 The case with n = [20] 

3.1.4 Photon emission [5JJ 

3.1.5 Tracing over the field [M] 

3.1.6 Different initial states [35] 

3.1.7 Conclusions [30] 

3.2 Entanglement swapping between spacelike separated atoms . . . [31] 

3.2.1 Introduction [3T] 

3.2.2 Entanglement swapping between spacelike separated atoms [33] 

3.2.3 Conclusions [37] 

3.3 Entanglement Sudden Death and Sudden Birth. Photon ex- 
change and correlations transfer in atom-atom entanglement dy- 
namics [35] 

3.3.1 Introduction [33J 

3.3.2 Hamiltonian and state evolution [TO] 

viii 



3.3.3 Sudden death and revival of atom-atom entanglement . . 2JJ 

3.3.4 Atom-field entanglement EH 

3.3.5 Tripartite entanglement [48j 

3.3.6 Conclusions 

Circuit QED as a quantum simulation of matter-radiation in- 
teraction beyond RWA 1511 

4.1 Detecting ground state qubit self-excitations in circuit QED: 

slow quantum anti-Zeno effect [52] 

4.1.1 Introduction E2] 

4.1.2 Detecting ground state qubit self-excitations [53] 

4.1.2.1 The Rabi model [53] 

4.1.2.2 Detecting excitations with one measurement . . [55] 

4.1.2.3 Repeated measurements: survival probability . . [57] 

4.1.2.4 Numerical experiments [55] 

4.1.2.5 Weak measurements [53] 

4.1.2.6 Relaxation and dephasing [62] 

4.1.3 Conclusions 12 

4.2 Entanglement dynamics via propagating microwave photons . . . [55] 

4.2.1 Introduction [55] 

4.2.2 Superconducting qubits coupled to a quantum field .... [57J 

4.2.3 Entanglement dynamics and single photons [70] 

4.2.4 Experimental implementation [71] 

4.2.5 Conclusions 1731 

Causality in matter-radiation interactions 1741 

5.1 The Fermi problem with artificial atoms in circuit QED [74] 

5.1.1 Introduction [JJ] 

5.1.2 There are no causality problem in Fermi's two-qubit system [75] 

5.1.3 Experimental proposal [75] 

5.1.4 Conclusions [82] 

5.2 Quantum fluctuations and short-time quantum detection [83] 

5.2.1 Introduction [S3] 

5.2.2 What does a detector's click mean? [531 

Quantum simulations of relativistic dynamics. Majorana physics 1871 

6.1 On Majorana Hamiltonians [55] 

6.1.1 Introduction [55J 

6.1.2 Majorana Hamiltonians [55J 

6.1.3 Hamiltonization [50] 

6.1.4 Examples [92] 

6.1.5 Conclusions and outlook [M] 

6.2 Quantum simulation of Majorana equation and unphysical oper- 
ations insi 

6.2.1 Introduction [55] 

6.2.2 Quantum simulation of antiunitary operations [96] 



ix 



6.2.3 Simulating antiunitary operations [551 

6.2.4 Conclusions ITUU1 

6.3 Quantum simulation of relativistic potentials without potentials . 1101 1 

6.3.1 Introduction 1 101 1 

6.3.2 One particle systems 1 101 1 

6.3.3 Examples H051 

6.3.4 Bipartite systems 11061 

6.3.5 Conclusions 11061 

a una 
b run 

C HT31 

D Majorana equation in 1+1 dimensions 11171 

E Implementation of the quantum simulation of the Majorana 
equation in trapped ions 11201 



Conclusions 11221 



Bibliography H~26l 



x 



Chapter 1 

Introduction 



In their famous paper of 1935, Einstein, Podolsky and Rosen (EPR) tried to 
show that Quantum Mechanics couldn't be a complete theory due to the pre- 
diction of strange correlations appearing among the parties in some states of 
composite systems [T]. Later that year, in a letter to Einstein, Schrodinger 
referred to the phenomenon with the German word "verschrankung" , that was 
translated by himself to "entanglement" in [2j . Considered unphysical and para- 
doxical by EPR - "no reasonable definition of reality could be expected to permit 
this"-, Bohr [3] and Schrodinger [2] immediately acknowledged entanglement as 
an essential feature of the theory. Today, after the series of experiments based 
on Bell inequalities 4 ruling out local hidden variable theories [5], Quantum 
Mechanics and its relativistic field version Quantum Field Theory are not con- 
sidered just as mere mathematical models that fit with outstanding precision 
a huge bulk of experimental data, but also a complete theory of Nature at the 
atomic and subatomic level. Being at the very heart of the theoretical building 
of Quantum Mechanics, entanglement is thus an essential feature of Nature. 
Besides its importance from this foundational viewpoint, entanglement is also a 
key resource for Quantum Information and Quantum Computation tasks. 

In the early days after its discovery entanglement was considered as a con- 
sequence of the interaction. For instance, Schrodinger wrote [2]: 

"When two systems, [...] enter into temporary physical interaction 
due to known forces between them, and when after a time of mutual 
influence the systems separate again, then they can no longer be 
described in the same way as before. [...] By the interaction, the two 
representatives (or ^-functions) have become entangled" 

Therefore, the "magic" of entanglement relies completely on the quantum state, 
which is generated by a direct interaction between the parties. Today, we know 
that this is not a complete picture. For instance, in the "entanglement swap- 
ping" protocol [5] the parties that get finally entangled do not interact directly 
with each other. Instead if we want to entangle A and B with this protocol, A 
has to share an entangled state with some third party C and also B indepen- 
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dently with the same C. Then, a measurement of the state of C leaves A and B 
in an entangled state. But in this case the origin of the A-B entanglement can 
be traced back to the origin of the A-C and B-C entanglements which, following 
Schrodinger would be again the interaction. 

In some sense, the entangling interaction can be regarded as a"black box" 
that generates the desired initial entangled state. In typical experiments involv- 
ing entanglement, the black box may be a nonlinear crystal exhibiting Sponta- 
neous Parametric Down-Conversion [7] . A "pump" laser beam is directed to the 
crystal and some of the photons split into entangled pairs, which are then used 
in the experiment. But what it is inside the box? Where it is exactly the magic? 
Can we have a deeper look onto the origin and the generation of entanglement? 
These are some questions that we will try to analyze in this Thesis. 

Actually, the above is deeply connected with modern investigations on the 
intriguing notion of "vacuum entanglement" . That the vacuum of the field is an 
entangled state was discovered by Summers and Werner in the 1980's [5], but it 
was considered as a mere formal theoretical result until it attracted some atten- 
tion from a different perspective [3]. In [pj, vacuum entanglement is exploited 
to generate entanglement between two spacelike separated detectors. After a 
finite time of interaction with a scalar field initially in the vacuum state, the 
two-level detectors state is shown to evolve from an initial separable state to an 
entangled one, even if the detectors remain spacelike separated. The standard 
interpretation of this effect is that the entanglement initially contained in the 
vacuum state of the field is transferred to the detectors, like in the standard 
entanglement swapping described above but replacing C by two spacelike sep- 
arated regions of the field and without measurements. An alternate view was 
given in [TU], where the entanglement between spacelike separated detectors is 
attributed to the exchange of virtual photons between them. Whatever the 
point of view, the mathematical structure underlying the phenomenon is the 
same, namely the Feynman propagator of the field. The fact that this object is 
not restricted to the lightcone was noticed by Feynman himself 

An important part of this Thesis is devoted to going a step further the 
previous results on entanglement generation between spacelike separated objects 
[TU1 [5J [T2] by moving to the physical framework of matter-radiation interaction 
and performing a thorough analysis of the phenomenon, keeping also in mind 
possible applications in realistic Quantum Information protocols. In a nutshell, 
one of the main goals was to propose the first realistic experimental proposal to 
test these effects in the lab. Unfortunately, as we will explain with more detail 
in the main text, a lot of experimental difficulties arise when dealing with real 
atoms and photons. This is the main reason to come to the different but related 
framework of circuit QED. 

In circuit QED, superconducting qubits can play the role of artificial two- 
level atoms that can be coupled to transmission lines along which photons propa- 
gate. Thus, these systems mimic the standard matter-radiation interaction with 
the advantage of experimental amenability. Although they are interesting by 
themselves, circuit QED setups can be regarded as ID Quantum Simulations of 
Quantum Optics setups. In this thesis, we will exploit the analogy to propose 
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the desired experimental test of "outside the lightcone" entanglement . Besides, 
we will explore some features of the novel physics showing up in this emerg- 
ing field, as a consequence of the possibility of achieve very large values of the 
qubit-field coupling strength and the breakdown of the standard Rotating Wave 
Approximation. 

A natural question showing up when dealing with effects as those described 
above is: What about causality? Is it violated? Indeed the relationship between 
entanglement and causality has been present since the EPR paper. At first sight 
it may seem that the perfect correlations appearing for instance in the paradig- 
matic singlet state represent by themselves an instantaneous transmission of 
information. This seems to be also the viewpoint of Einstein himself when 
he write to Born the famous "spukhafte Fernwirkung" -commonly translated as 
"spooky action at a distance". Today, we know that Quantum Mechanics is 
actually a nonsignaling theory p~3l [14] , that is, the statistics of the measure- 
ments performed by B is completely independent of the measurements realized 
by A if A and B are spacelike separated. This entails that information between 
A and B cannot be transmitted at superluminal rates. In the singlet state, B 
can only get information on the state of A if both parties share some previous 
information on the shared entangled state and on the choice of the observable 
measured by A. In general, is well known in Quantum Information that quan- 
tum correlations have to be assisted with Classical Communication -which is of 
course subluminal- in order to transmit information. Therefore, the fact that 
a separable state can evolve to an entangled one "faster than light" does not 
represent a violation of Einstcnian causality, since cannot be used by itself to 
transmit information. 

Although this general principle is clear, the translation to the particular 
setup consisting of a pair of neutral atoms interacting with the electromagnetic 
field -the one that we extensively deal with in this Thesis- has been an open 
theoretical problem -the so-called Fermi problem- since the 1930's. In partic- 
ular, in 1932 Fermi [15] proposed a "gedanken" experiment to check that the 
probability of excitation of one of the atoms behaves in a causal fashion (see 
Chapter 5 for more details). Although Fermi's conclusions were conceptually 
right about the causal nature of the model, a mathematical flaw of his compu- 
tations opened a Pandora box that couldn't be closed neither by a theoretical 
consensus nor by an experimental test [161 117j . An important part of this The- 
sis is devoted to show that there are no causality problems in Fermi two-atom's 
system and to propose a feasible experimental test, taking advantage again of 
the possibilities offered by circuit QED. The results of such experiment could 
close an eight-decade old controversy on the foundations of Quantum Mechanics 
and Quantum Field Theory. 

As we have explained above, circuit QED setups play a central role in this 
Thesis, as ID Quantum Simulators of matter-radiation interaction. In a nut- 
shell, a Quantum Simulator, as Feynman envisioned |18| is an experimentally 
amenable quantum device that mimics the dynamics of an unaccessible quan- 
tum model. The emerging field of Quantum Simulations is developing very 
fast a number of such devices in several branches of Physics. In particular, in 



3 



this Thesis we pay attention to Quantum Simulations with trapped ions of Rel- 
ativistic Quantum Mechanics systems. [THj. Relativistic Quantum Mechanics 
can be understood as a Quantum Field Theory with a fixed number of particles, 
thus it amounts to Quantum Field Theory below the energy threshold of pair 
creation. In this low-energy regime is the right model for relativistic particles. 
Surprisingly, its equations -such as Klein-Gordon or Dirac ones- come with a 
wealth of striking predictions, like "Zitterbewegung" or Klein paradox, which 
are very far from a direct detection in the lab [20] [21] . In this sense, recent 
experiments with trapped ions simulators of the Dirac equation [22] have shed 
light on the nature and experimental implications of the mentioned phenomena. 
An important goal for the future would be to simulate the fully relativistic ver- 
sion of the two-atom model used in most part of this Thesis where the atoms 
would be replaced by relativistic particles following the Dirac equation. Steps 
in this direction are given at the end of this Thesis, where bipartite systems 
of relativistic interacting particles are considered. In the meanwhile, we take 
a little diversion from the main course and propose the quantum simulation of 
a crucially different but related one- particle relativistic dynamics, namely the 
Majorana equation- see details in chapter 6-. The very interesting properties of 
this equation -which takes its roots in the prematurely lost overwhelming genius 
of Ettore Majorana [23] but was actually introduced in [23] [25]- lead us to the 
introduction of a new class of mathematical objects with physical content, that 
we refer to in this Thesis as "Majorana Hamiltonians" . 

Besides the Majorana equation, this new class of generalized Hamiltonians 
is connected with a wider range of applications, going from the implementation 
of operations like time reversal and charge conjugation to the partial transpose- 
a crucial formal operation in Quantum Information due to its relationship with 
entanglement measures-. These operations are "unphysical" , in the sense that 
they are not associated to operations that can be implemented directly in the 
laboratory. Instead, they entail a mathematical formal computation that has 
to be performed for instance in a classical computer and require the knowledge 
of the quantum state. With the techniques developed in this Thesis these op- 
erations can now be simulated with real physical operations, without stopping 
an ongoing experiment. From a fundamental viewpoint the similarities and dif- 
ferences between Majorana Hamiltonians and standard textbook Hamiltonians 
allows us to explore important questions lying at the foundations of Quantum 
Mechanics and Quantum Field Theory. 

Structure of the Thesis. Brief summary of objec- 
tives and results. 

The structure of this thesis is the following. 

• In chapter 2 we will just provide the reader with some theoretical back- 
ground on the formalism used in the main text. 
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• In chapter 3 we will analyze different aspects of the entanglement dynam- 
ics in a system consisting of a pair of two-level neutral atoms with an 
electric dipole coupling to the electromagnetic field. We mainly focus in 
the relationship of entanglement with the existence or lack of a causal 
connection between the atoms. 

* In section 3.1 we analyze whether a pair of neutral two level atoms 
can become entangled in a finite time while they remain causally dis- 
connected. The interaction with the electromagnetic field is treated 
perturbatively in the electric dipole approximation. First, we start 
from an initial state in which the field is in the vacuum and only one 
atom is excited. We obtain the final atomic correlations for the cases 
where n = 0, 1, or 2 photons are produced in a time t, and also when 
the final field state is unknown. Our results show that correlations 
are sizable inside and outside the mutual light cone for n = 1 and 
2, and also that quantum correlations become classical by tracing 
over the field state. For n = we obtain entanglement generation by 
photon propagation between the atoms, the correlations come from 
the indistinguishability of the source for n = 1, and may give rise 
to entanglement swapping for n = 2. Finally, we consider a similar 
scenario but starting with both atoms excited as initial state. 

* In section 3.2 we show a mechanism that projects a pair of neutral 
two-level atoms from an initially uncorrelated state to a maximally 
entangled state while they remain spacelike separated. The atoms 
begin both excited in a common electromagnetic vacuum, and the 
radiation is collected with a partial Bell-state analyzer. If the inter- 
action time is short enough and a certain two-photon Bell state is 
detected after the interaction, a high degree of entanglement, even 
maximal, can be generated while one atom is outside the light cone 
of the other, for arbitrary large interatomic distances. 

* In section 3.3 we analyze the entanglement dynamics of a system 
composed by a pair of neutral two-level atoms that are initially en- 
tangled, and the electromagnetic field, initially in the vacuum state, 
within the formalism of pcrturbative quantum field theory up to the 
second order. We show that entanglement sudden death and revival 
can occur while the atoms remain spacelike-separated and therefore 
cannot be related with photon exchange between the atoms. We 
interpret these phenomena as the consequence of a transfer of atom- 
atom entanglement to atom-field entanglement and viceversa. We 
also consider the different bi-partitions of the system, finding similar 
relationships between their entanglement evolutions. 

• In chapter 4 we move to the framework of circuit QED with artificial 
atoms, considering these setups as Quantum Simulators of a ID version 
of matter-radiation interaction. 
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* In section 4.1 we explore a consequence of the breakdown of the Ro- 
tating Wave Approximation in the ultrastrong coupling regime of cir- 
cuit QED. In particular, we study an ultrastrong coupled qubit-cavity 
system subjected to slow repeated measurements. We demonstrate 
that even under a few imperfect measurements it is possible to de- 
tect transitions of the qubit from its free ground state to the excited 
state. The excitation probability grows exponentially fast in analogy 
with the quantum anti-Zeno effect. The dynamics and physics de- 
scribed in this section is accessible to current superconducting circuit 
technology. 

* In section 4.2 we propose a simple circuit QED experiment to test the 
generation of entanglement between two superconducting qubits that 
arc initially in a separable state. Instead of the usual cavity QED 
picture, we study qubits which are coupled to an open transmission 
line and get entangled by the exchange of propagating photons. We 
compute their dynamics using a full quantum field theory beyond 
the rotating-wave approximation and explore a variety of regimes 
which go from a weak coupling to the recently introduced ultrastrong 
coupling regime. Due to the existence of single photons traveling 
along the line with finite speed, our theory shows a light cone dividing 
the spacetime in two different regions. In one region, entanglement 
may only arise due to correlated vacuum fluctuations, while in the 
other the contribution from exchanged photons shows up. 

• In chapter 5 we show that the nonlocal quantum correlations phenomena 
considered in previous chapters are perfectly compatible with Einstenian 
causality by showing explicitly that Fermi's two-atom system behaves in 
a causal way, and explore the theoretical consequences of going beyond 
RWA in a quantum detection. 

* In particular, in section 5.1 we propose a feasible experimental test 
of a 1-D version of the Fermi problem using superconducting qubits. 
We give an explicit non-perturbative proof of strict causality in this 
model, showing that the probability of excitation of a two-level arti- 
ficial atom with a dipolar coupling to a quantum field is completely 
independent of the other qubit until signals from it may arrive. We 
explain why this is in perfect agreement with the existence of nonlocal 
correlations and previous results which were used to claim apparent 
causality problems for Fermi's two-atom system. 

* In section 5.2 we study the information provided by a detector click 
on the state of an initially excited two level system. By computing 
the time evolution of the corresponding conditioned probability be- 
yond the rotating wave approximation, we show that a click in the 
detector is related with the decay of the source only for long times 
of interaction. For short times, non-rotating wave approximation ef- 
fects like self-excitations of the detector, forbid a naive interpretation 
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of the detector readings. These effects might appear in circuit QED 
experiments. 

• Finally, in chapter 6 we consider different quantum simulations, namely 
of relativistic quantum mechanics setups. 

* In section 6.1, we introduce the notion of Majorana Hamiltonians, 
objects that are in general neither linear nor antilinear but with 
an associated dynamics which conserves the norm. The Majorana 
equation, that is, the relativistic quantum mechanical equation of a 
fermion with Majorana mass term, is the archetypal example of this 
family. We analyze some amazing properties like the relevance of ini- 
tial global phases in the dynamics that makes them detectable and, 
as a consequence, the inadequacy of the notion of mixed states. We 
provide an alternative description in terms of standard Hamiltonians, 
showing that a Majorana Hamiltonian can be mapped to a Hamil- 
tonian in a real Hilbert space. In this way, Majorana Hamiltonians 
and antiunitary operators can be experimentally demonstrated. 

* In section 6.2 we show how to design a quantum simulator for the 
Majorana equation, a pseudohamiltonian - Majorana-hamiltonian- 
rclativistic wave equation that might describe neutrinos and other 
exotic particles beyond the standard model. The simulation demands 
the implementation of charge conjugation, an unphysical operation 
that opens a new front in quantum simulations, including the discrete 
symmetries associated with complex conjugation and time reversal. 
Finally, we show how to implement this general method in trapped 
ions. 

* Finally, in section 6.3 we show that a potential can be simulated with 
a free Dirac or Majorana Hamiltonian since there is a local phase 
transformation between solutions of the equation with a potential 
and solutions of the free equation. The transformation depends on 
the potential, which is codified in the phase. In some cases, the 
probability density is unchanged, entailing that the particle behaves 
under the potential as if it were free. It is valid for a large class 
of potentials for the Majorana equation and a different class in the 
Dirac case, which explains why Dirac and Majorana particles exhibit 
different behavior under the same potentials. We extend the results 
to two-body equations with interaction potentials. 

This thesis concludes with five appendices. The first three develop with more 
detail some involved computations related with chapters 3, 4 and 5. Finally, ap- 
pendix D is devoted to a theoretical explanation on the 1-D Majorana equation 
and appendix E sheds light onto the implementation of the quantum simulation 
of the Majorana equation with trapped ions. 
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Chapter 2 

Preliminaries 



2.1 Matter-radiation interaction. Hamiltonian 
2.1.1 Introduction 

Throughout this work we will deal with a system consisting of one or more two- 
level systems (qubits in the language of Quantum Information) interacting with 
the electromagnetic quantum field. This section is devoted to provide a brief 
summary of the theoretical grounds of the employed model, with no pretension 
of completeness. Wider textbook treatments on these topics can be found in 

H>], [17], [in, |2g. 

We will start with the well known minimal-coupling Hamiltonian of an atom 
interacting with a quantum field: 

// = E^(P^* A ( X ^) 2 + 1/ (I X - X *I) + E / d 3 kh^(at x a kX + ^), (2.1) 

i 1 A 

Pi, x.j, being the momentum and position, respectively of each particle of mass 
TOj and charge qi bound by a potential V to a force center (nucleus) located 
at x. The quantum electromagnetic field can be written down in terms of the 
creation and annihilation operators for each mode k (with frequency ujk = c |k|) 
and polarization A: aj^, a^A as: 



M^) = ] j 2c£o h {27r)3 T, J ^(e ikXi e (k,A)a fcA + e - ikx - e *(k,A)4 A ) (2.2) 

e(k, A) being the polarization vectors for each A and assuming the commutation 
relationships 

[okA.oLvl =S 3 (k-k / )6 xy . (2.3) 
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(2.1| would be the quantized version of the classical minimal-coupling Hamilto- 



nian, which can be derived from a Lagrangian whose Euler-Lagrange equations 
give rise to the correct equations of motion of the system, namely the Maxwell 
equations for the field and the Lorentz equations for the particles. Notice that 
we can split (2.1) into a "free" part Hq and an interaction part Hj : 



H = H + Hi 



(2.4) 



with 

i 

and 

But then: 



Pf 

2 TO," 



g, 2 A 2 (^) 
2to 



Hi 



f d 3 kftw k (4 A a kA + ~), (2.5) 
Qi Pi ■ M x i) 



E 



Ho ^ Hqr + H p, 



(2.6) 



(2.7) 



that is, the "free" part of the Hamiltonian cannot be decomposed into the free 
Hamiltonian of the radiation, 



H r = ^ / d 3 k/iu>k(a kA akA + \) 



and the Hamiltonian of the particles without field: 

^2 



H 



OP 



^-Bi- + y(|x-x 2 

^— ' 2 to 



(2- 



(2.9) 



There is an extra, "self-interaction" term: iiAJ^A^ Thus the eigenstates of 



Hq (2.5 1 are not product states of eigenstates of Hqr and Hop. To avoid this, 



instead of using the Hamiltonian in (2.1 1 we will move to a different equivalent 
description of the system. 

To this end we start by assuming the so-called "long wavelength approxima- 



tion", which amounts to replacing A(xj) (2.2) by A(x) in the Hamiltonian (2.1 ). 
This approximation is only valid for those modes verifying k • |x — Xi| << 1, 
that is, for wavelengths much larger than the atomic size. Then, we can apply 
the following unitary transformation: 



U 



e » 



d(x).A(x) 



where we have introduced the electric dipole moment d: 

d(x) = ^?i(xi -x). 



(2.10) 



(2.11) 
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This unitary transformation is an example of displacement operator, whose 
properties are well known in the framework of coherent states of the electro- 
magnetic field. Using these properties and with a little algebra it is possible to 
arrive to the transformed Hamiltonian: 



H' = UHU' 1 , 

and assuming that we arc dealing with a globally neutral atom, 



(2.12) 



(2.13) 



H' takes the form: 

H' = H 0R + H QP + V / + H' I . (2.14) 
The first two terms are the free Hamiltonians of the particles and the field 



respectively p^8|)-(|2^9|. The third term V is: 

d 3 k- 



V' = 



A 



1 



2 e (2^) 3 



( £ -d(x)) 2 



(2.15) 



a divergent contribution which cancels out with one of the terms coming from 
the energy shift caused by emission and reabsorption of photons, as we will see 
in the next chapter. Thus, we can let this term drop. The last term is an 
interaction Hamiltonian: 

H'j = -d(x) • E(x) (2.16) 
where E(x) is the electric field vector operator: 



E(x) = l V2^W^ / d3fc ^( e4kX£ ( k ' A ) a ^- e " kX£ * ( k > A )4A)- 



So, finally we can split the new Hamiltonian H' (2.14) into: 

H' =H' + H'j 



but now 



OR 



H, 



OP, 



(2.17) 
(2.18) 
(2.19) 



in contrast with (2.7). Obviously, H and H' are equivalent descriptions of the 



system since they are related through a unitary transformation (2.10), (2.12) 



but states have to transform also according to the unitary. Throughout this 



work, we will the use the Hamiltonian given by (2.16), (2.18), (2.19), whose 



'bare" eigenstates are product states of eigenstates of the free Hamiltonians 
for the particles and the field. Equivalent results would be obtained with the 



Hamiltonian in (2.1 ), provided that the states are correctly transformed through 
(2~10|. 
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2.1.2 Two level-atoms. Interaction picture 

Now we will particularize to the case of two-level atoms and move to the inter- 
action picture. We consider that we have an atom with a ground state \g) and 
an excited state |e) such that: 



H \g)=E g \g) ,H \e)=E e \e) 



and 



E e Eg — hVi. 



Under that conditions, we can use the completeness relation 

1 = |.9> (9\ + \e) (e|, 



(2.20) 
(2.21) 

(2.22) 



in order to write the electric dipole moment operator in a simple form. That is: 



d(t) 



(\9) (Sl + |e) (e|)e^de^(| 5 ) (g\ + |e) (e|) 

t n t 



= e * d ge a~ 



snt , 



with 
and 



(9\d\e) 
a \a <<\ .. + = \e){g\. 



(2.23) 

(2.24) 
(2.25) 



Thus the interaction Hamiltonian (2.161 in the interaction picture takes the 
form: 

flJ(i) = -d(x,t)-E(x,t), (2.26) 



with d(x, t) given by (2.23) and the standard expression for the electric field 
vector operator in the interaction picture: 

E(x,i) = i 



h c 



2e (2 



— Y,f d3k (e J ( kx -^*)e(k,A)a fcA - 

n ' A 

e - l (kx- Wfc t) £ * ( k;A ) a t^ (2.27) 



In the next chapter, we will extensively use the Hamiltonian ( |2.26 1, and a TD 
version of it will be considered in chapters 4 and 5. 



2.1.3 Rotating wave approximation and beyond 



The Hamiltonian in (2.26) can be written as a sum of four terms using (2.23) 
and (2.27), each of them containing two operators: an atomic ladder operator 
cr^ and a creation or annihilation Fock operator ak\, <xt\. Thus we have two 



OfcA, 



terms a al\, <r + ak\ associated to the frequency ft — to and two terms a 
a + a k)v associated to the frequency fl + u>. The rotating wave approximation 
(RWA) let the latter so-called counter-rotating terms drop, retaining only the 
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first two ones. The usual argument is that the counter-rotating terms oscillate 
very quickly and average to zero in observable time scales. This is generally true 
in the Quantum Optics realm, due to the weak atom-field coupling. As we will 
see in this thesis (chapters 4 and 5), the argument fails for stronger couplings, as 
those appearing in circuit QED, where counter-rotating terms can be relevant 
at experimental times. Even for weak couplings, a theoretical analysis of the 
short-time behavior like the one in chapter 3 must take the non-RWA terms into 
account. Therefore, throughout this thesis the RWA will not be assumed. 



2.1.4 Rabi and Jaynes-Cummings models 



When the boundary conditions of the system are such that the electromagnetic 
field is restricted to a given region of space, engineered in such a way that the 
fundamental frequency of the field loq approximately matches with the frequency 
of the qubit transition, then the multimode expansion of the electric field in 



the Hamiltonian (2.16) (or (2.261 in the interaction picture) can be restricted 



to only one mode wq. In that case, the Rabi model is obtained. If the RWA is 
applied to the Rabi model, the resulting hamiltonian is the celebrated Jaynes- 
Cummings model, describing a two-level system in an optical cavity. We will 
deal with the Rabi model in circuit QED in chapter 4. 



2.2 Circuit QED 

Although the oscillations in electrical circuits are usually described by classical 
mechanical laws, at micro and nano scales and ultra-low temperatures, the ther- 
mal noise can be negligible compared to the spacing and width of the energy 
levels, unveiling the quantum nature of the system. Moreover, non-linear ele- 
ments such as Josephson junctions do not obey the principle of correspondence, 
that is, the averages of the momentum and position operators do not follow the 
classical equations of motion, making apparent the necessity of a fully quantum- 
mechanical description of such elements. Thus non-trivial quantum effects can 
be observed in systems that can be considered as macroscopic, in the sense that 
they contain a huge number of elementary particles -the quantum numbers being 
collective degrees of freedom such as flux or charge-. Among all the possibilities 
opened up by this novel approach to Quantum Mechanics, in this thesis we are 
concerned with just one: circuit Quantum Electrodynamics. 

A systematic quantization procedure of the usual elements of a circuit, such 
as inductances, capacitances and resistances is given in |30|, 13 1] within the ap- 
proach of lumped-element circuits, valid when the relevant wavelengths are much 
larger than the element's size. This can also be extended to waveguides -also 
called transmission lines- which can be regarded as a chain of infinitesimally 
small LC lumped-element circuits, obtaining a quantized Hamiltonian which is 
the 1-D version of |2.8| With this lumped-element approach is also possible to 
obtain a quantum description of several types of superconducting qubits, us- 
ing different combinations of superconductors and Josephson junctions [521 155] . 
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Engineering the length of the transmission line to match the fundamental wave- 
length of the qubit transition, a cavity QED system like the ones described in 
section 2.1.4 is obtained |34j|35]. Without this restriction the full hamiltonian 
of the matter-radiation interaction emerges j3SJ GH] , but in a range of experi- 
mental parameters different from the typical one in Quantum Optics. In this 
artificial ID version of atomic QED, the frequency of the qubit transition is 
in the order of Ghz, and thus the relevant wavelengths of the field are in the 
microwave regime. 

In chapters 4 and 5 of this Thesis, we will exploit the possibilities of circuit 
QED as a quantum simulator of matter-radiation interaction, taking advantage 
of its experimental amenability, which includes the possibility of achieve very 
strong qubit-field couplings, entering into a regime in which non-RWA effects 
can be observed in the lab. 

2.3 Entanglement measures in two-qubit systems 

In this thesis, we will analyze the dynamics of the entanglement between qubits, 
like the two-level atoms described in the previous sections. In this section we 
give a brief description of the notion of entanglement and discuss different ways 
of quantify it. 

Let us start with pure states of two-qubit systems, which take the following 
general form: 

ihj = e 7<?) being a basis of the Hilbert space Ha ®Hb for qubits A and 
B. These pure two-qubit states are called separable if it can be written as: 

l*> = l*>A l*>m (2-29) 

\*b) B being states of the Hilbert spaces Ha, Hb, respectively. Otherwise, 
the state is referred to as entangled, and exhibit correlations between the parties 
that cannot be achieved classically. For the non-pure case, a state p is separable 
if they can be written as 

p = Y.p^®p^ ( 2 - 3 °) 

i 

where pf and pf are state operators of subsystems A and B respectively, Pi > 
Vi, J2iPi = 1- Otherwise p is entangled. 

In general, an entanglement measure is some function of the state that is zero 
for separable states and non-zero for entangled ones. Besides, some technical 
properties must be fulfilled, the most standard one being the so- called mono- 
tonicity under Local Operations and Classical Communication (LOCC): these 
operations cannot create entanglement, thus any entanglement measure must 
non-increase under them. There are several measures of the entanglement of 
pure bipartite states; we shall cite only three of them: von Neumann's entropy 
of reduced states, concurrence (38], and negativity. [39] . 
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The Von Neumann's entropy of a state p is defined in Information Theory 

as 

S(p) = -T /Pj \og 2 p j , (2.31) 



3 



being pj the eigenvalues of p. Reduced states pa, p B are the result of tracing p 
over qubits B and A respectively: 

PA =Tr B p,PB =Tr A p. (2.32) 

Von Neumann's entropy of reduced states is the simplest measure of bipartite 
entanglement for pure states, but its extension to non-pure states does not 
correctly discriminate between separable and entangled states. Therefore, it is 
not a good measure of entanglement for general mixed states. 

Concurrence is the more usual measure of entanglement for non-pure two- 
qubit states. It is defined as 



C{p) = max{0, VAi - V\ 2 - ~ V M> ( 2 -33) 
where {Aj}, (i = 1...4) are the eigenvalues of 

R = pp (2.34) 

in decreasing order, with 

P = (<Ty (g) (Ty)p*(<Jy ® (J y) , (2.35) 

and p* being the complex conjugate of the state operator p. In arbitrary di- 
mensions, operational generalizations of concurrence are known only for pure 
states. Therefore, it is a good measure of entanglement for general two-qubit 
states (pure or not) and for pure states in any dimension, but not for mixed 
bipartite states of arbitrary dimensions. 

Finally, the negativity of a bipartite state p, defined in [3D] as: 

N{p) = -2Y,°i{p T % (2.36) 

i 

where {&i(p TA )} are the negative eigenvalues of the partial transpose p TA of 
the total state p respect to the subsystem A, defined as 

(iAjB\p TA \kA,lB) = {k Al jB\p\iA,lB), (2.37) 

A and B denoting the two subsystems: this is twice the value of the original 
definition of the negativity in [39 J . It can be proved that for pure bipartite states 



of arbitrary dimensions the negativity (2.371 is equal to the concurrence 



The main advantage of the negativity is that it can be evaluated in the same 
way for pure and non-pure states in arbitrary dimensions, although there are 
entangled mixed states with zero negativity in every dimensions except 2x2 
(two qubits) and 2 x 3 (a qubit and a qutrit) [41], [42]. Thus, negativity is 
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not in general an ideal measure of bipartite entanglement of non-pure states: 
no measure discriminating separable from entangled states in the general non- 
pure case is known. However, non-zero negativity is a sufficient condition for 
entanglement and for distillability [IT] (a state is said to be distillable if a 
maximally entangled state can be obtained from it through (LOCC)), a very 
important property for Quantum Information tasks. 
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Chapter 3 

Dynamics of entanglement 
in matter-radiation 
interaction 



"A devil with merely local powers like a parish vestry would be too 
inconceivable a thing" (Sir Arthur Conan Doyle, The hound of the 
Baskervilles) 

In this chapter we will analyze exhaustively the entanglement dynamics in a 
system consisting of a pair of two-level atoms interacting through a quantum 
electromagnetic field. In Section 3.1 we start from different initial separable 
states and consider different projections onto final states of the field. In Section 
3.2 a protocol for generating highly entangled states between spacelike separated 
atoms is considered. Finally, in section 3.3 we start from an initial entangled 
state and study the phenomena of entanglement sudden death and entanglement 
sudden birth. Throughout all the chapter, entanglement measures are used to 
compute the evolution of entanglement and relate its behavior to the spacetime 
region in which the atoms are placed. 

3.1 Generation of entanglement inside and out- 
side the mutual light cone 

3.1.1 Introduction 

Quantum superposition and entanglement are the cornerstones lying at the foun- 
dations of quantum information and the principal support of the new quantum 
technologies which are at different stages of conception and development at 
present. Putting entanglement to work, enabling its use as a resource, is the 
key to the success of these technologies. Therefore, a complete understanding 
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of entanglement, necessary at the fundamental level, is also important for these 
developments to occur. 

As we have explained in Chapter 1, entanglement can be envisaged in very 
different forms; it originally appeared in quantum mechanics PQ as a direct 
connection between distant particles, a residue of past direct interaction between 
them [2] . In quantum field theory entanglement can be traced back to the non- 
locality of the vacuum state [5J 143) or, simply, to field propagation. Similar 
arguments operate for a lattice of coupled oscillators [44] . 

In this section we analyze some features of entanglement generation closely 
related to the microscopic causality of quantum field theory. Put in simple 
words, this work attempts to ascertain whether a pair of spatially separated 
parties (say, a pair of neutral two-level atoms A and B) can get entangled in 
a finite time while they remain causally disconnected [HI H3 [46] . Each party 
interacts locally with the electromagnetic field, the carrier of the interaction. We 
stress here that, as shown in perturbation theory produces non-signalling 
P3| results for this system and that the apparent causality violations come from 
the nonlocal specification of some final states. At first sight, the question can 
be answered in the negative; if the parties remain causally separated from each 
other, they can not entangle. However, the propagator D(x,y) is finite even 
when c{x — y)° < |x — y|, and perhaps some correlations could be exchanged 
between both parties [TU] ■ Alternatively, the correlations could be blamed on the 
preexisting entanglement between different parts of the vacuum [SJ I12j , which 
could be transferred to the atoms. Whatever the point of view, correlations 
are exchanged through (time ordered) products, while only commutators are 
restricted to be causal [47]. Our analysis can not sidestep that the role of the 
field goes beyond that of a mere carrier, quanta could be absorbed from the field 
or escape in the form of photons [3H1 SS] • How does the entanglement between 
A and B depend on the state of the field? This question shapes our discussion 
below. 

We will include in the final state all the perturbatively accessible field states, 
analyzing for each of them the correlations in the reduced atomic state. We 
compute the entanglement measures for different values of (x — y)° and |x — y|, 
that lie inside the atoms mutual light cone and beyond. The atomic state that 
results after tracing over the states of the field is separable, which means, in 
the scheme of [9], that there is no transference of vacuum entanglement, only 
classical correlations. In [5], these correlations become entanglement when a 
suitable time dependent coupling with the scalar field is introduced. As pointed 
out in [12] , this would require an unrealistic control of the atom-field interaction 
in the electromagnetic case that we are dealing with here. As an alternative way 
to achieve entanglement between the atoms we consider a post-selection process 
of the field states with n = 0,1,2 photons. This is a nonlocal operation and 
therefore entanglement generation is allowed. In [ID] , only the vacuum case 
when |x — y| 3> c(x — y)° was analyzed, and no entanglement measures were 
considered. We get quantum correlations for all the different field states. We 
also get useful hints on the nature of the correlations, whether they come from 
photon exchange, source indistinguishability, etc. 
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3.1.2 The model 



We will consider the field initially in the vacuum state, including the cases with 
0, 1 and 2 final photons to analyze perturbatively the amplitudes and density 
matrices to order a. We assume that the wavelengths relevant in the interaction 
with the atoms, and the separation between them, are much longer than the 
atomic dimensions. As seen in the previous chapter, under these conditions the 
Hamiltonian can be split into two parts 

i? = #o + (3.1) 

that are separately gauge invariant. The first part is the Hamiltonian in the 
absence of interactions other than the potentials that keep A and B stable, 

Ho = H a pA + Hqpb + H r- (3.2) 



HqpAi Hqpb being the corresponding "free" Hamiltonians (2.9) of particles A 



and B, respectively, and H r is given by (2.8). The second contains all the 
interaction of the atoms with the field 

Hj = - d n (i)E(x n ,t), (3.3) 

n=A,B 



where E is the electric field ( 2.27[ ), and 



d„ = ^2 q t (xj - x„) (3.4) 





= dAge 


d B 


= dsge 


d A \ 


= |d B | 



is the electric dipole moment (2.11) of atom n, whose matrix elements (2.24) 
we will take as real and of equal magnitude for both atoms 

dAeg 
dBeg 

(3.5) 

| e ) and | g ) being the excited and ground states of the atoms, respectively. 
In what follows we choose a system given initially by the product state, 

|^(0)) = \eg0) (3.6) 

in which atom A is in the excited state | e ) , atom B in the ground state | g ) , 
and the field in the vacuum state | ) . The system then evolves under the effect 
of the interaction during a lapse of time t into a state: 

I V (*)> = So dt'Htit')/*]^ (3.7) 

(T being the time ordering operator) that, to order a, can be given in the 
interaction picture as 

latom^atoma, field) t = ((1 + A) \ eg) + X \ ge)) \ 0) 
+(U A \gg)+V B \ee))\l) + (F\eg)+G\ge))\2) (3.8) 
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where 



.4 
U A 
F 



-(0|T(5+57+5 B 5+)|0), A 
V B = (115+10) 



(Q|T(5+5^)|0) 



;i|5i|o; 



(3.9) 



= -(2|T(5+S7 +5 B 5+)|0), G= (2|T(5+57)|0>, 



being 



S n = -- J dt'H I {t')=S+ +S- 
S n=-\ f Q dt'e^a* d„E(x,i), 



(3.10) 



(the sign of the superscript in 3.10 is associated to the energy difference between 
the initial and final atomic states of each emission) and \n), n = 0, 1, 2 is a 
shorthand for the state of n photons with arbitrary momenta and polarizations, 
i.e. 



1) (II = 



A 



d 3 k|k, e(k,A))(k, e(k, A) |, 



(3.11) 



etc. fi is the frequency of the atomic transition, as defined in (2.21 1 



Among all the terms that contribute to the final state (3.8) only X corre- 
sponds to interaction between both atoms, which is "real" interaction only if 
ct > r (r being the interatomic distance). This would change at higher order 
in a. Here, A describes intra-atomic radiative corrections, Ua and Vb single 
photon emission by one atom, and G by both atoms, while F corresponds to 
two photon emission by a single atom. Details on the computations of these 
quantities would be given in Appendix A. 

In Quantum Optics terms like Vb, F and G are usually neglected by the 
introduction of a RWA. But here we are interested in the short time behavior, 
and therefore all the terms must be included, as in [46l [50l [51] . Actually, only 
when all thesenon-RWA effects are considered, it could be said properly that 
the probability of excitation of atom B is completely independent of atom A 
when r > ct [50J [5T] (r being the distance between the atoms) . 



Note that the actions in (3.10) depends on the atomic properties f2 and d r , 



and on the interaction time t. In our calculations we will take (f2|d n |/ec) = 
5 • 10~ 3 , which is of the same order as the Is — > 2p transition in the hydrogen 
atom, consider Qt > 1, and analyze the cases (r/ct) ~ 1 near the mutual light 
cone, inside and outside. 

Given a definite field state | n ) the pair of atoms is in a pure two qubits 
state as shown in (3.8). We will denote these states by | AB) n , then 



p%l = \AB) nn (AB | 



(3.12) 



and the corresponding reduced states by 



(n) 

Pa 



rp „(«) An) 
J TB P AB , P B 



Tr A Paw 



(3.13) 
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in the following, and will compute the entropy of entanglement S^ n ^ 

sj(W) rr, (n) n (n) 
§ v ' = Tr py log py 

and the concurrence [551: 



max 



Ai being the largest of the eigenvalues Aj (j = 1, ...,4) of 

-RAS = [(cTy ® Oy)PAB * (o-y ® CT; / )]pAB 

for them. 

3.1.3 The case with n = 



1.x 1CT 6 



8.X10 



(3.14) 
(3.15) 

(3.16) 




Figure 3.1: Concurrence of the atomic state in the e.m. vacuum p^ B as a 
function of a; = (r/ct) for three values of z = (fir/c) — 5 (solid line), 10 
(dashed line) and 15 (dotted line). The height of the peak is C*- ^ = 1. 

We first consider the case n = 0, where the field is in the vacuum state and, 



after (3.8), the atoms are in the pure state 

|AB) = ((1 + A)\EG)+X\GE))/co, 



(3.17) 



where Co = */fl + A\ 2 + \X\ 2 is the normalization, giving a concurrence 

= 2\X\\1 + A\/c 2 . (3.18) 
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It is interesting to note that at lowest order the concurrence arises as an ef- 
fect of the mutual interaction terms X mediated by photon exchange or, in a 
different language, by the vacuum fluctuations. As expected, at higher orders 
the radiative corrections described by A dress up these correlations. Analytic 
lowest order calculations ([53]) showed that they can persist beyond the mu- 



tual light cone, vanishing for x = (r/ct) oo. We sketched in Fig. 3.1 the 
concurrence C 1 - ' for x around 1. Our computations were done for the illus- 
trative case where both dipoles are parallel and orthogonal to the line joining 
A and B. We will adhere to this geometrical configuration for the rest in the 
following. It would correspond to an experimental set up in which the dipoles 
are induced by suitable external fields. C^ ^ shows a strong peak (of height 1) 
inside a tiny neighborhood of x = 1. The features outside the mutual light cone 
are $(|d|/er) 2 ~ 1CP 6 here, and could be larger if Qt < 1 entering into the 
Zeno region (incidentally, \X\ oc t A for very small t [53]). Notice the change of 
behavior between the region where the atoms are spacelike separated (x < 1) 
and the region where one atom is inside the light cone of the other (x > 1). 
This quantitative treatment complements the qualitative one given in [10 . 
The entropy of entanglement written in terms of the small quantity 

% = (|A|/c ) 2 e(0,l) (3.19) 

is 

§ (0) = -(1 - m) Iog(l - Vo) - Vo log/*) , (3.20) 

this is a positive quantity in (0, 1), which attains its maximum possible value 
§(°) = i when the state is maximally entangled at 770 = 0.5. This is well within 
the small neighborhood of x = 1 mentioned above. Radiative corrections would 
shift the maximum to \X\ = 1 1 + A| , so the entropy is sensitive to the Lamb 
shift when this contributes to the dipole radiative corrections. 



3.1.4 Photon emission 

We now come to the case n = 1, where the atoms excite one photon from the 
vacuum, jumping to the state 

\AB) 1 = (U A \gg) + V B \ee))/c u (3.21) 

(with ci = a/| Ua | 2 + I Vb | 2 ), during the time interval t. The density matrix 
for this case contains the term 

l = V B U* A =Tr 1 {l\S+\0){l\S^\0r = (0|5+5+|0), (3.22) 

producing a concurrence 

C (1) = 2\l\/c\ , (3.23) 

so, even if this case only describes independent local phenomena attached to 
the emission of one photon by either atom A or B, the concurrence comes from 
the tangling between the amplitudes u and v which have different loci. The 
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Figure 3.2: Concurrence for one photon final state (3.23) as a function of x — 
r/ct for three values of z = Qr/c = 5 (solid line), 10 (dashed line) and 15 
(dotted line). Entanglement vanishes as t — > oo (x — ¥ for a given r) and is 
sizeable for x > 1. 



state of the photon emitted by A and the state of A are correlated in the same 
way as the state of the photon emitted by B with the state of B are. These 
independent field-atom correlations are transferred to atom-atom correlations 
when we trace out a photon line with different ends, A and B, when computing 
V U* . In fact, while \U\ 2 and \V\ 2 are independent of the distance r between 
the atoms, 

I = -SMl m _ h ?. )M » [r) + {Si . + (3.24) 



where 



which depends explicitly on r. Above we use d S 1 ^) — sin( to t/2)/(nuj), which 



M(r)= I dfc^^<5*(Q + cfc)<5*(fi-cfc) (3.25) 



becomes S(lo) in the limit t — > oo. In Fig. 3.2 we represent in front of 
x = r/ct for some values of z = flr/c. As the Figure shows, there may be a 
significative amount of concurrence for all x, indicating that is an entangled 
state inside and outside the mutual light cone. The peak at x = 1 comes from 
the term with phase k(r—ct) that can be singled out from the linear combination 



of phasors in the integrand of (3.25) 
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Here we have a lone photon whose source we cannot tell. It might be A or B, 
with the values of I and C^ 1 - 1 depending on their indistinguishability. Eventually, 
RWA will forbid the process g — > e + 7 for large interaction times. Therefore, 
Vb, I and C' 1 - 1 will vanish as t grows to infinity (x — > for each value of z in 
Fig. 2), as can be deduced from the vanishing of <5* (O + ck) for t 00. 

The entropy of entanglement gives an alternative description of the situation. 
Its computation requires tracing over one of the parts A or B, so no information 
is left in §W about r, but it still gives information about the relative contribution 
of both participating states |ee) and \g g) to the final state. In terms of 

m = \v\ 2 / Cl 2 e (0,1), (3.26) 

we have 

= -(1 - 7?!) io g (i - m) - m logm (3-27) 

Would not be for the difference between £l + ck and £1 — ck, Vb should be equal 
to Ua, "Hi = 0.5, and would attain its maximum value. Not only this is not 
the case but, as said above, Vb will vanish with time and only \g g ) will be in the 
final asymptotic state. Notice the result, indistinguishability was swept away 
because for large t we know which atom (A) emitted the photon. Therefore, the 
entropy will eventually vanish for large interaction times. 

There are two possibilities with n — 2; one (with amplitude F) when both 
photons are emitted by the same atom, the other (with amplitude G) when each 
atom emits a single photon. The final atomic state 

\AB) 2 = (F\eg) + G\ge))/c 2 , (3.28) 

with c 2 = Vl^l 2 + \G\ 2 , is in the same subspace as for n = 0. The normal- 
ization C2 is 0(a) like the expectation values F, G, so that all the coefficients 

(2) 

in may be large. The concurrence is 

C (2) = 2\FG* |/4 (3.29) 

Due to the tracing over photon quantum numbers, F G* is a sum of products 
containing not only factors Ua and Vb, but also r dependent factors like /. The 
entropy §' 2 ^ is now given in terms of a parameter 

m = \G\ 2 /cl (3.30) 

Notice that 

|G| 2 = |{/4| 2 |V B | 2 + |;| 2 . (3.31) 

Hence, both C^ 2 ) and S^ 2 \ depend on r. This is different from the single photon 
case, where the only r dependence was in the coherences of Pab> which did not 
feed into p^ . The correlations came in that case from the indistinguishability of 
the photon source. The case n — 2 resembles that of the entanglement swapping 
paradigm [6], where there are two independent pairs of down converted photons. 
Here we have two independent atom - photon pairs. The swapping would arise 
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in both cases from detecting one photon of each pair. But with the initial state 
we are considering here, both F and G eventually vanish. More interesting 
would be the case with the initial atomic state | e e ) , that will be considered in 
section 3.1.5. 



3.1.5 Tracing over the field 



We have seen that if the state of the field is defined, the atomic state is entangled 
inside and outside the light cone. But what happens if the field state is ignored, 
that is, if we trace over the field degrees of freedom? Then the atomic state is 
represented by the following density matrix (in the basis {\ee),\eg),\ge),\g g)}): 



PAB 



( \Vb? I \ 

|1 + A| 2 + |F| 2 (1 + X)*+FG* 

X(l + A)*+F*G |X| 2 + |G| 2 

V i* o o \u A \ 2 J 



N- 1 (3.32) 



where I = (0\S^S^ |0) was used again, and N = |1 + A\ 2 + \X\ 2 + \U\ 2 + \V\ 2 ■ 
\F\ 2 + \G\ 2 . 
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Figure 3.3: Mutual information of p^g as a function of x — r/ct for z — 
ttr/c = 5 (solid line), 10 (dashed line), 15 (dotted line). The inset shows the 
finite concurrences that are possible only for small values of x. 



The state in (3.32) is an example of the so-called X-states, and the concur- 
rence of such a state is given by: 

2 

C{pab) = - max{|p 2 3| - y/pnpu , \pu\ ~ \/p22P33 , 0} . (3.33) 
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Numerical computations show that the concurrence associated to this density 
matrix always vanishes except for a bounded range of small values of x. Beyond 
this range pab is a separable state with no quantum correlations, either inside 
or outside the light cone. But the atoms A and B are mutually dependent even 
for zero concurrence. Their mutual information 

Kpab) = Hpa) + S(pb) - S(pab), (3-34) 

which measures the total correlations between both parties, is completely clas- 
sical in this case, but may be finite. We show this quantity in Fig. |3.3| for 
different values of z with an inset with the concurrence for small values of x. 

3.1.6 Different initial states 
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Figure 3.4: Concurrence of the atomic state in the e.m. vacuum p^ B as a 
function of x — {r/ct) for z — (ilr/c) = 5 (solid line), 10 (dashed line) and 15 
(dotted line). The height of the peak is — 1. x —> (t — >• oo) is the region 
usually considered in Quantum Optics. 

In what follows we choose a system given initially by the product state, 

|^(0)) = \ eeO) (3.35) 

in which atoms A and B are in the excited state | e ) and the field in the vacuum 
state | ) . The system then evolves under the effect of the interaction during a 
lapse of time t into a state: 

| V- (t)> = e~ l Si dt ' H ' | V> ) (3.36) 
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Figure 3.5: Concurrence of the atomic state in the e.m. vacuum p AB if 
\ee) is the initial state as a function of z — (fir/c) for z/x = Q,t = 6 (solid 
line), 9 (clashed line) and 12 (dotted line), x > 1 amounts to z > fit in each 
case. 



that, to order a, can be given in the interaction picture as 



|atomi,atom 2 , field) t = ((1 + A') \ ee) + X' | g g)) \ 0) 
+(U A \ge) + U B \eg))\l) + (F'\ee)+G'\gg))\2) 



(3.37) 



where 



.4' 

17a 
F' 



= -(0\T(S+S A +S+S B )\Q), X 1 = (0\T(S B S A )\0) 



i|^|o; 



U B 



MSb\0) 



(3.38) 



2\T(S+S A +S+S B )\Q), G> = {2\T{S B S' A )\Q) 



The objects has been defined in (3.10). As in the previously analyzed case, 
in (3.37) only X' corresponds to interaction between both atoms, and now is a 
completely non-RWA term. A' describes intra-atomic radiative corrections, Ua 
and Ub single photon emission by one atom, and G' by both atoms, while F' 
corresponds to two photon emission by a single atom. 

We begin with the case n = 0, where the field is in the vacuum state and, 
following (3.37), the atoms are in the projected pure state 

\AB) a = {{\ + A')\ee)+X'\GG))/c 0) (3.39) 
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Figure 3.6: Concurrence for one photon final state if |e e ) is the initial state 
as a function of a; = r/ct for three values of z = Vlr/c. The values of C for x > 1 
are of the same order as those displayed in Fig. . 

where c = + A'\ 2 + \X'\ 2 is the normalization, giving a concurrence 

C (0) = 2\X'\\1 + A'\/c 2 . (3.40) 

The computation of A' and X' can be performed following the lines given in 
Appendix A where they were computed for the case of a initial atomic state 
\eg). We will consider that the dipoles are parallel along the z axis, while the 
atoms remain along the y axis. Under that conditions, using the dimensionless 
variables x = r/ct and z = £lr/c: 

A' = AtICz l (i n h_ ^1 + 2t7r), 
3 a; z 

X' - ad 4 L {-V 2 5 V] +V l V 3 )I, (3.41) 
with K = a | d | 2 /(e 2 r 2 ) and 1 = 1+ + where: 

I± = ~ % l ' J [±2 cos( -) e ± 1 z Et(Tiz)+ e' lz{1± ^ 
2 z x 

Ei(iz(l±-)) ~ e lz{1± * ] Ei{-iz{l± -))] (3.42) 

x x 

for x > 1 with the additional term — 2nie lz ^ 1 ~ 1 ^ x ' > for x < 1. We use the 
conventions and tables of [54] . 
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Figure 3.7: Concurrence for one photon final state if |ee) is the initial 
state as a function of z = ilr/c for three representative values of the time 
Sit = 8, 10, 12 with peaks at z = 8, 10 and 12 respectively. 



We show in Fig. 3.4 the concurrence CW (|3.40l for x near 1 for given values 
of z. Like the case where | e g ) is the initial atomic state, jumps at x = 1 and 
has different behaviors at both sides. In Fig. |3.5| the concurrence is sketched as 
a function of z for given values of f2 1 — z/x. The tiny values of the concurrence 
for the region z > fit (which corresponds to x > 1), diminish as t grows and 
will eventually vanish, since X' is a non-RWA term. 

In the case n = 1 the final one photon atomic state would be 



\AB) 



(U A \ge) + U B \eg))/c' 1 , 



(3.43) 



where U a is the same function of x^ as Ub is of x^, and c[ — ^/2| U a | 2 • Now, 
the indistinguishability of the photon source commented in Section 3.1.3 persists 
for large t and so does entropy and concurrence. In particular, the r-dependent 
concurrence is 

C (1) = 2|/'|/(c' 1 ) 2 , (3.44) 



where 



/' = 



0) 



(3.45) 



We represent it in Fig. 3.6 In Fig. 3.7 we have represented the concurrence 
for the case where the initial atomic state was |e e ) in terms of the inter-atomic 
distance for three fixed values of time. What we obtain is a shift of the concur- 
rence features to longer r as t grows (so that they appear at the same (r/ct)), in 
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Figure 3.8: Concurrence of the atomic state with n — 2 photons p^ B in 
front of x = r/(ct) for z = Qr/c = 5 (solid line), 10 (dashed line) and 15 
(dotted line). 

such a way that, even if t is just the duration of the interaction, it plays the role 
of propagation time for the generated correlations. They are negligible small 
for large r, peak at the "light cone" but, on the other hand grow, as we would 
expect, for larger interaction times. 

Now we will focus on the two photon case. The final atomic state 

\AB) 2 = (F'\eg) + G'\ge))/c 2 , (3.46) 

with c 2 = y/\ F' | 2 + | G' | 2 , is in the same subspace as for n = 0. The 
normalization c 2 is 0(a) like F', G", so that all the coefficients in may 
be large. Therefore, although the probability of attaining this state is small, 
the correlations are not. The concurrence is 

C (2) =2\F'G'*\/cl (3.47) 

Wc find that: 

F' = 0(h - t 2 )( V A (h) U' A (i 2 ) + U A (*i) V'a fe) 

+ V B (ti) U' B (t 2 ) + U B (ti) V'b (h) ) (3.48) 
G' = U B U' A + U A U' B 

with Va = ( 1 1 S A | ) and Vb = ( 1 1 Sg | ) . The primes account for the two 
single photons, i.e. 

|2) = |ke A , k'e v ). (3.49) 
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Figure 3.9: Concurrence of the atomic state with n = 2 photons p^ B in 
front of z = Qr/c for fit = z/x = 6 (solid line), 9 (dashed line) and 12 (dotted 
line) . 



The quantities | U A | 2 = | U B | 2 = |[/| 2 , \V A \ 2 = \V B \ 2 = \V\ 2 ,l = 
XJaV'b = U B VX, UV* = U* A VX = f/sVg, f7 B C^l and V^Vg can be com- 
puted follo wing the lines inAppendix A. 



In Fig. 3.8 we show C^ 2 ' in front of x for given values of z. When x — > 
(t — > oo, i.e. the Quantum Optics regime), F' vanishes and the final atomic 
state would be the separable state \gg), with zero concurrence. Entanglement 
is sizable for x > 1, and could be maximized if a particular two photon state 
was detect ed, a s we will explain in the next section. 



In Fig. 3.9 , C (2) is sketched as a function of z for given values of O t = z/x. 
Again, the concurrence for the region z > Ctt (x > 1), diminish as t grows 
and will eventually vanish, since it is due to /, which is a non-RWA term. 
Interestingly, as we noted for the single photon emission, x — 1 is a singular 
point that divides the spacetime into two different regions. This occurs even if 
t is not the propagation time of any physical signal between the atoms. This 
effect comes from the appearance of effective interaction terms like I, that would 
be missing if we could discriminate the source of emission of each photon. 

3.1.7 Conclusions 

In this section we have studied the correlations between a pair of initially separa- 
ble neutral two-level atoms that are allowed to interact with the electromagnetic 
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field, initially in the vacuum state. We have computed the concurrences that 
arise when the final state contains n — 0, 1, or 2 photons. They may be sizable 
for x small (t — > oo for a given r) and also around x = 1. Only in the case n = 
there are interactions between both atoms, generating an entanglement that 
may persist asymptotically in the case that \eg) is the initial atomic state. We 
have carefully taken into account all the terms contributing to the amplitude for 
finite time (they are oc t 4 , not oc t 2 as is sometimes assumed). A small amount 
of entanglement can be generated between spacelike separated parties due to 
the finiteness of X when x > 1, but a change of behavior appears for x < 1. 
For n — 2 the final atoms are in the same subspace than for n = 0. There are 
similar correlations that in this case can give rise to entanglement swapping, by 
measuring both photons in a definite state for instance. Naturally, in this case 
entanglement may be sizable for spacelike separated parties, as here this is not 
related to any kind of propagation. Entanglement in the case with only a final 
photon (n — 1) comes from the indistinguishability of the photon source. It will 
vanish asymptotically in the case that the initial states is \eg) when eventually 
only one atom [A in the present case) may emit the photon, and it is also sizable 
when x > 1. It is interesting how these correlations become classical (except for 
small x) when the states of the field are traced over. We have shown through 
the mutual information the residues of what were quantum correlations in the 
individual cases analyzed before. 

3.2 Entanglement swapping between spacelike 
separated atoms 

3.2.1 Introduction 

There are mainly two different known ways of generate entanglement: by in- 
teraction between the atoms (for instance, [55]) or by detection of the emitted 
photons 48, 56, 57, 58, 49 . Some of these proposals have been realized experi- 
mentally (for instance, [59]). For the latter cases, in principle, there is no reason 
to expect that the swapping [5] of atom-photon to atom-atom entanglement can 
only begin to occur when one atom enter into the light cone of the other. 

The possibility of entanglement generation between spacelike separated atoms 
is of both theoretical and practical interest, and was addressed from different 
points of view in [10j|9l[12] and in the previous section of this Thesis. In section 
3.1, we analyze this issue perturbatively in a simple model of a pair of two- level 
atoms interacting locally with the electromagnetic field, initially in the vacuum 
state. Tracing over the field states, the atoms are only classically correlated, but 
applying | n)( n\ (n — 0, 1, 2 being the number of photons up to second order in 
perturbation theory), the atoms get entangled. For n = the entanglement is 
generated by the interaction term and therefore is only relevant when one atom 
enter into the light cone of the other, despite of the finiteness of the Feynman 
propagator beyond that region. But for n = 1,2 entanglement may be sizeable, 
although small, if the interatomic distance is short enough. In [S], the trace 
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over the field states was considered in a model with a pair of two-level detectors 
coupled to a scalar field. The detectors may get entangled if a suitable time de- 
pendent coupling is introduced, and this was applied to a linear ion trap in [12] . 
In [TU], only the vacuum case when t — > was analyzed, and no entanglement 
measures were considered. 

As commented in chapter 1, there are two possible interpretations for these 
effects: as a transfer of preexisting entanglement of the vacuum [3 \T% or as 
a consequence of the propagation of virtual quanta outside the light cone [10] . 
Both are compared and discussed in [TO] , 

In this section we will go one step further and consider that the photons are 
detected with definite momenta and polarizations. We show that, in principle, a 
high degree of entanglement, even maximal, can be generated between spacelike 
separated atoms if a Bell state of the emitted photons is detected. We will 
consider a pair of neutral two-level atoms separated by a fixed and arbitrary 
distance and study the evolution of an initially uncorrclated state under local 
interaction with the electromagnetic field. We focus on the two-photon emission 
which, although has a smaller probability of success, shows a larger fidelity of 
the projected state with the desired state and has a entanglement robust to 
atomic recoil [60]. The photons pass through a partial Bell-state analyzer [BTj . 
and we use entanglement measures to study the evolution of entanglement in the 
projected atomic states after detection of the different photonic Bell states. The 
results show that interaction times must be short, but interatomic distances can 
be as large as desired. The interaction time is independent of the photodetection 
time, which is only related with the distance from the atoms to the detectors. 
That distance can be such that the photodetection can occur while the atoms 
remain spacelike separated. 

The results can be interpreted as a transfer of part of the vacuum entan- 
glement after a post-selection process. If no measurement were performed the 
atoms would have classical correlations transferred by the vacuum. In [9] the 
classical correlations may become entanglement with a suitable time dependent 
coupling. The post-selection process can be seen as an alternative way to achieve 
the entanglement transference. While the results in [101 IS] are mainly theoret- 
ical, these could be probed experimentally, and would show for the first time 
the possibility of transfer entanglement from the vacuum state of the quantum 
field to spacelike separated atoms. 

3.2.2 Entanglement swapping between spacelike separated 
atoms 



In what follows we choose a system given initially by the product state, 

|^(0)) = |ee0) (3.50) 

in which atoms A and B are in the excited state | e ) and the field in the vacuum 
state | ) . As in the previous section, the system then evolves under the effect 
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Figure 3.10: Schematic setup for the entanglement swapping described in the 
text. The atoms A and B are at (y,z) — (=Fr/2, 0). The emitted photons pass 
through a 50 : 50 BS at (0, r/2) and two PBS at (± d/2y/2, L/2 + d/2\/2), and 
there are four single photon detectors at the outport ports of the two PBS, at 
(±d/y/2,r/2 + d/V2) and (±d/V2,r/2). Taking into account that |*~) and 
|$ _ ) are forbidden in our model, a |\& + ) is detected when there are coincidence 
clicks in two detectors and |<J> + ) when there is a double click in one detector. 



Then the atoms are projected onto the atomic part of the state (3.52). 



of the interaction (3.3) during a lapse of time t, and, up to order e 2 , 0, 1 or 2 
photons may be emitted. If after that a two-photon state is detected, 

|*) = |photo ni ,photon 2 ) = c fefe',A,A' l^ A > fc ~' A ') ( 3 - 51 ) 

fe,fe',A,A' 

(being hk, Hk' momenta and A, A' polarizations), the projected state, up to 
order e 2 , can be given in the interaction picture as 

|photons,ato mi , atom 2 ) * = | ^ ( ^i ee ) + G*\gg) ^ (3 52) 



where 



F* = ^\T(S+S A + S+S B )\0), = (^\T(S B S A )\0) (3.53) 
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and N = y/\ | 2 + | G* | 2 . The objects S± were denned in p.10 ). 

Here, describes single photon emission by both atoms, while Fq, corre- 
sponds to two photon emission by a single atom. In Quantum Optics, Fy is 
usually neglected by the introduction of a rotating wave approximation (RWA), 
but as we will see later, for very short interaction times Fy and G$ may be of 
similar magnitude. Actually, a proper analysis of this model can be performed 
only beyond the RWA, as we have commented in the previous section. Without 
RWA vacuum entanglement cannot be transferred to the atoms with this par- 
ticular post-selection process. In that case, a one photon post-selection process 
would entangle the atoms. 



(3.53) can be written as: 



= ^(ii-i 2 )(*|5+(t 1 )5 A (t 2 )+5+(t 1 )5 B (t 2 ))|0>, 
G* = (*\S B (h)S A (t 2 ))\0) (3.54) 

The photons pass through a partial Bell-state analyzer |61j consisting in a 
beam splitter (BS) and two polarization beam splitters (PBS) with four single 
photon detectors at their output ports. If two detectors, one at one output port 
of one PBS and one at an output port of the other, click at the same time, a 
state \^~) is detected, while if the two clicks are in the two output ports of only 
one PBS, the state is | V I /+ ). If one of the four detectors emits a double click, the 
state can be |$ + ) or l^ - ). Taking into account momenta and symmctrization, 
the Bell states can be written as 

= _L[|fc|, Pt) + |fc't, fc|> 

± (|fct, k'l) + \k'l, fct»] (3.55) 
|$±) = ^=[\ki,9i) + \i?i,ki) 

± (|fct, fc't) + |fc't, fct»] 
where j" and J, are the photon polarizations, with polarization vectors 

e(fc,t) = ^=(e(M) + e(M)) (3.56) 

and 

e(k,i) = -L( e (fc,l)- e (fc,2)), (3.57) 

where 

e(fc,l) = (cos 6>fc cos 0fc , cos sin 0fc , — sin6>fc) (3.58) 

and 

e(k,2) = (-sin0 fc ,cos0 fe ,O). (3.59) 

Here 

\kX,k'X')=a{ x a{, x ,\0). (3.60) 
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We will use the concurrence (2.33 1 C to compute the entanglement of the atomic 
states when the different Bell states are detected. The concurrence of the atomic 



part of a state like (3.52) is just given by 



2\ 

N~ 2 



(3.61) 



We assume that the atoms A, B are along the y axis, at y = TL/2 respectively, 
and the dipoles are parallel along the z axis, corresponding to an experimental 
set up in which the dipoles are induced by suitable external fields [10] . We also 
take 

\k\ = | A? | = fi/c (3.62) 

Under that conditions, the first remarkable thing is that for and |$ _ ), 

we have = G^ = 0. Therefore, at least while only El transitions are 
considered, in this model the Bell-sate analyzer is complete: if two different 
detectors click the state is | \E r_l " } , while if one detector clicks twice the state is 
|$ + ). First, we focus on \^ + ). Considering (3.54) and (3.55), with the mode 



expansion for the electric field and the commutation relation for the creation 
and annihilation operators, a standard computation leads to: 



Gq/ 



-i(fif)oos(-MM)) 



K(Q,t,6) cos (-MM)), 



with 



K(n,t,t 



ca\d\W 
2ir 2 e 2 



sin v k sin Ok' 



(a being the fine structure constant), 

| - 1 + e 2im (l-2iflt)\ 



and 



(nty 



(sin k sin 4> k ± sin 9 k > sin (j> k , ) , 



(3.63) 



(3.64) 



(3.65) 



(3.66) 



6 k , <t>k corresponding to k and k ', 4>k' to k', and z = Qr/c. Notice that j(Clt) 
decreases as t grows, and eventually vanish as t — » oo as required by energy 
conservation. The r dependence is a result of the individual dependence on the 
position of each atom, not on the relative distance between them. 
Taking into account ( |3.61 1 and (3.63) the concurrence is given by: 



4| cos ( 



lh + (9, 0))cos(§M0, 0))l 



cos- 



(§M0, 0)j(fi*) 



cos- 



{%h-{6, 0)) 



(3.67) 



j(nt) 



Now, we assume that the 50:50 BS is at (y,z) = (0,r/2), the two PBS at 
(±d/2V2,r/2 + d/2V 2) an d the f our detectors at (±d/v / 2,r/2 + d/y/2) and 
(±r/\/2, /2) (see Fig. 3.10). (3.63) will not depend on the value of d, which is 
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Figure 3.11: Concurrence for the atomic state when a Bell state \^ + ) or 
of the photons is detected, as a function of x = r/ct for z = flr/c =1 (solid), 
5 (dashed), 10 (dotted). The light cone is at x < 1. For x > 1 the interaction 
time is short enough to have a significative amount of entanglement. 



the distance traveled by the photon to any detector aft er lea ving the BS. Notice 



3.11 



we represent (3.201 



that, with this setup, h- = and h+ = v2- In Fig 
under that conditions as a function of £ = r/ct for three different values of 
z (different values of r). Notice that a high degree of entanglement, maximal 
for x large enough (short enough interaction times t), can be achieved in all 
cases when one atom is beyond the light cone of the other (x > 1). As t — > oo 
(x — > 0), the concurrence eventually vanish, in agreement with the fact that 
the only atomic state allowed by energy conservation is just the separable state 
\9,9)- 



In Fig. 3.12 we represent (3.67) as a function of z for three different values of 
zjx = Clt, to give an alternative description. The mutual light cone corresponds 
to the region z < Q t in each case. The concurrence oscillates with the position 
of the atoms, and eventually vanish at z = v2 (n + 1/2) n (n — 0, 1, 2...), as a 
consequence of the vanishing of cos {z/y/2). For a given interaction time t, the 
maximum of the concurrence can be achieved for interatomic distances as large 
as desired. In particular, a maximally enta ngled state is ge nerat ed for fit = 1, 
which corresponds to t ~ 10~ 15 s. In Fig. 3.13 we sketch (3.67) as a function 
of <j) = 4>k = 4>k' for given values of x and z. Notice that the maximum values 
for the entanglement are around <p = nir/2 (n = 0,1,2...), n/2 corresponding 
to the setup of Fig. |3~T0l 
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z 

Figure 3.12: Concurrence for the atomic state when a Bell state \^ + ) of the 
photons is detected, as a function of z = Qr/c for Sit =1 (solid), 4 (dashed), 7 
(dotted) . The light cone for each curve is at z < £1 1 . 

So far, we have focused on |^+), but, in principle, |$ + ) could be detected as 
well. The coefficients and would have opposite sign to those of |^+) and 
therefore the concurrence would be the same. But, due to the interaction times 
considered here, the relaxation time of a single detector must be extremely short 
in order to emit a double click. 

3.2.3 Conclusions 

In this section, we have shown that, in principle, two neutral two-level atoms can 
evolve from an initially uncorrelated state to a highly entangled state in a time 
shorter than the time required for the light to travel between them. At the initial 
time, both atoms are excited in a common electromagnetic vacuum. They are 
allowed to interact with the field due to an induced dipole during a time t and, 
up to second order in perturbation theory, n = 0, 1, 2 photons may be emitted. 
After that, the emitted radiation pass through a partial Bell-state analyzer. For 
interaction times t ~ 10~ 15 s and if a two-photon Bell state \^> + ) or |$ + ) (the 
other two are forbidden in this model) is detected after that, the atoms are 
projected into an entangled state, which may be maximally entangled for short 
enough t. For a given t, the degree of entanglement oscillates periodically with 
the distance and the maximum degree available can be achieved for interatomic 
distances r as large as desired. Notice that the interaction time t, which must be 
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Figure 3.13: Concurrence for the atomic state when a Bell state |* + ) of the 
photons is detected, as a function of <j> for z = Qr/c — 5 and x = rc/t = 2.5. 



t ~ 10 -15 s, is absolutely independent of the time t' at which the photodetection 
takes place. Since the distance traveled by the photons from the atoms to the 
detector is r/^/2+d, d being arbitrary, the photodetection can occur after a time 
t' ^ r/c. A suitable choice of d is necessary in order to ensure that the atoms 
may remain spacelike separated. The degree of entanglement is independent of 
d. 



3.3 Entanglement Sudden Death and Sudden Birth. 
Photon exchange and correlations transfer 
in atom-atom entanglement dynamics 

3.3.1 Introduction 

Entanglement between qubits may disappear in a finite time when the qubits 
interact with a reservoir. This is commonly known as "entanglement sudden 
death" (ESD). After its discovery [S31 [M] , the phenomenon has attracted 
great attention (for instance, (S3 El E3 E3 EH E3 ED E2 E2 ) and has been 
observed experimentally [74] . 

ESD shows up in a variety of systems that can be roughly divided in two sets: 
those in which the qubits interact individually with different reservoirs and those 
in which they interact with a common environment. In particular, in [6711681 [75] 
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a system of a pair of two-level atoms interacting with a common electromagnetic 
vacuum is considered. The dynamics of the system is given in all the cases by the 
Lchmbcrg-Agarwal master equation |76l 177] which is derived with the rotating 
wave approximation (RWA) and the Born-Markov approximation. Recently, 
non-Markovian [69] and non-RWA [7H] effects have been considered in systems 
of qubits coupled individually to different reservoirs. There are good reasons 
for going beyond the Markovian and RWA scenario in the case of a pair of two- 
level atoms in the electromagnetic vacuum. For short enough times non-RWA 
contributions are relevant (see Chapter 3) and a proper analysis of causality 
issues can only be performed if they are taken into account [46, 5U1IST]. Besides, 
as we shall show in this section the death of the entanglement between the atoms 
is related with the birth of entanglement between the atoms and the field, and 
therefore the field is actually a non- Markovian reservoir. This was also the case 
in [66l [711 E2] with different reservoirs. 

In the previous sections, we have applied the formalism of perturbative quan- 
tum electrodynamics (QED) to the system of a pair of neutral two-level atoms 
interacting locally with the electromagnetic field, and for initially separable 
states analyzed the generation of entanglement. This is a non-Markovian, non- 
RWA approach. The use of the Lehmberg-Agarwal master equation can be seen 
as a coarse-grained in time approximation to the perturbative treatment |79] . 
The first goal of this section is to apply also the QED formalism to analyze 
the ESD in these systems for initially entangled atomic states, comparing the 
results with the previously obtained |67L 168) with master equations. We will 
focus mainly on the range r/(ct) ss 1, r being the interatomic distance and t 
the interaction time, in order to investigate the role of locality. We will also con- 
sider for the first time in these systems the rest of pairwise concurrences, namely 
the entanglement of each atom with the field, and multipartite entanglement, 
following the spirit of [SSI IZH EH EI] • While the mentioned papers deal with a 
four qubit model, our model here consists in two qubits (the atoms) and a qutrit 
(the electromagnetic field, which may have 0, 1 or 2 photons). We shall show 
that the phenomenon of revival of entanglement after the ESD [58] can occur 
for r > ct, and therefore is not related with photon exchange as is usually be- 
lieved. We will see that atom-atom disentaglement is connected with the growth 
of atom-field entanglement and viceversa. Similar relationship will be obtained 
among the "atom-(atom+field)" and "field- (atom+atom)" entanglements. 

The reminder of this section is organized as follows. In subsection 3.3.2 
we will describe the Hamiltonian and the time evolution from the initial state 
of the system. In section 3.3.3 we will obtain the reduced state of the atoms 
and analyze the behavior of its entanglement. In section 3.3.4 the same will 
be performed with the reduced state of each atom and the field, comparing the 
entanglement cycle with the one obtained in the previous section. Tripartite 
entanglement will be considered in section 3.3.5 in terms of the entanglement 
of all the different bi-partitions of the system, and we conclude in section 3.3.6 
with a summary of our results. 
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3.3.2 Hamiltonian and state evolution 

In what follows we choose a system given initially by an atomic entangled state, 
with the field in the vacuum state | 0) : 

|^(0)) =(a|ee)+ 0\gg))\O). (3.68) 



The system then evolves under the effect of the interaction (3.3) during a lapse 
of time t into a state: 

| ^ ( t )) = T( e - 1 /o* dt ' H ' (*')/») | -0>o, (3.69) 



Up to second order in perturbation theory, (3.69 1 can be given in the interaction 
picture as 

|atoml, atom2, field(i)) = a | e e (i)) + f3 \ g g (i)) (3.70) 

where 

|ee0(t)> = ((l+A / )|ee)+X'| 55 ))|0) + (C7 A | 5 e) + C7 fl |e 5 ))|l) + 
(F'|ee)+G'| 5fl ))|2) (3.71) 

and 

\ggO(t))= ((l + A ")\gg)+X"\ee))\0) + (V A \eg)+V B \ge))\l) + 
(F"\gg)+G"\ee))\2) (3.72) 



A' = ~e(t 1 -t 2 )(0|«S+(t 1 )^(t 2 )+«S+(t 1 )^(t 2 )|0) 
^ = ^(ti-t 2 )(0|^(i 1 )5+(i 2 )+5 B (< 1 )5+(< 2 ))|0) 



where 

1 

2 
1 

2 

X' = (0|T(5 B ^)|0),X" = (0|T(5+5+)|0), 
(7 A = (1\S2\0),V A = (1|5+|0) (3.73) 
C/ B = (1|5 B |0), V B = (1|S+|0) 

= ^(*i-*2)(2|«Sj(*i)5i(i a )+5+(ti)5i(ta)|0) 

= ^(*i-*2)(2|5j(ti)5+(i a )+5B(*i)5+(t 2 ))|0), 
G' = (2|r( l S B 5 A )|0),G"-(2|T( l S+5+)|0) 

with the definitions and conventions of Section 3.1. Here, A' and A" describe 
intra-atomic radiative corrections, V a {Ub) and Va (Vb) single photon emission 
by atom A (B), and G' and G" by both atoms, while F' and F" correspond 
to two photon emission by a single atom. Only X' and X" correspond to 
interaction between both atoms. A', X' , X" , Va, Vg, F', F" and G' are non- 
RWA terms. As in the previous sections of this chapter, in our calculations 
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we will take (fi|d|/ec) = 5 • 1CP 3 , which is of the same order as the Is — » 2p 
transition in the hydrogen atom, consider Qt > 1, and focus mainly on the 
cases (r/ct) ~ 1. In this case, we will take into account | e) is actually a triply 
degenerate state | e , m) with m = 0, ±1 and we will average over two different 
independent possibilities for dipole orientations: = = d = du z for 
transitions with Am = [5U] and d = d (u x ± iu y )/^/2 [5U] for transitions with 
Am = ±1. 

3.3.3 Sudden death and revival of atom-atom entangle- 
ment 




x 



Figure 3.14: Concurrence C(pab) m front of x = r/ct for p = 0.98 and z — 
fir/c = 2 • 10 6 (solid line), 5 • 10 6 (dashed line) and 2 • 10 7 (dotted line). In the 
latter case sudden death and revival of entanglement occur for x > 1. 



After tracing over all the states of the field, the density matrix of the atomic 
state pab takes the form (in the basis {|ee), \eg), \ge), \g g)})- 



PAB 



1 

N 



( Pn 



V Pu 



P u\ 

P12 P23 

P*23 P33 

p 44 ) 



(3.74) 
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Figure 3.15: Concurrence C(pab) in front of x — r/ct for z — Or/c = 2 • 10 7 
and p — 0.97 (solid line), p = 0.98 (dashed line) and p = 0.99 (dotted line). In 
the first case, entanglement decreases as t grows up to a minimum value and 
begin to grow since then. This behavior becomes entanglement sudden death 
and revival when the minimum value is for higher values of p. C(pab) tends 
to as x — > oo and p — > 1. 



where 



where 



p n = | a (1 + A') + (3 X" \ 2 + \ a F' + (3 G" | 2 , 

P22 = P33 = \a\ 2 \U\ 2 + \f3\ 2 \V\ 2 + 2Re(af3*l*) 

p 44 = \ab + /3(1 + A") | 2 + \aG' + /?F"| 2 , (3.75) 

p 14 = \a\ 2 ((l + A')X'* +F'G'*) + \/3\ 2 ((l + A")* X" 

+ G" F"*)+a/3*((l + A'){l + A") + F' F"*) 

+ f3a*(X"x'* + G"G'*) 

P23 = \a\ 2 U B U* A + \p\ 2 V A V^ + 2Re(a/3*UV*) 

N = Pn + P22 + P33 + P44 



\U\ 2 = \U A \ 2 = |(7 B | 2 , 
\V\ 2 = \V A \ 2 = \v B \ 2 , 

l = U A V* B = U B VX (3.76) 
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and 



U V* =U* A Vt = U B V X 



(3.77) 



The computation of A' , X' , etc. can be performed following the lines given in 
the Appendix A, where they are computed for the initial state | eg) and only 
for Am = 0. In terms of z — Vtrjc and x = r/ct, being r the interatomic 
distance, we find: 



A' = 
A" = 



AiK z 3 

3 x 
-UK z 3 
3 x 



(In Cl- 
in (1 



-max 



Z 

^max 
Z 



) + tTr), 
) 



X' = =X"* = Q ^(-V%-+V 1 V J )7, 



(3.78) 



with K = a | d | 2 /(e 2 r 2 ) and 1 = 1+ + where: 
J± = 



l -^—-\±2cos{-)e ±lz Ei{Tiz) 
2 z x 



-iz(i±±: 



Ei(iz(l±-)) - e** (1± -'^(-i«(l± -))] 



(3.79) 



for x > 1, having the additional term — 27rie* z ' 1 otherwise. 

1 17 | 2 , | V | 2 , I, U B U%, V a V% and U V* have been computed in Section 3.1. 
Besides: 



G' = U B U' A + U A U' B , G" = V A V' B + V B V' a 

f' = e{h - t 2 )(v A (ti) u' A (t 2 ) + c/a (ti) v' A (t 2 ) 

+ V B {ti)U' B (t 2 ) + u B {h)v B (t 2 )) 

F" = 6(t! - t 2 )(U A (tx) VJi (t 2 ) + C/^ (t t ) V A (t 2 ) 

+ u B (h)v B (t 2 ) + u' B (h)v B (t 2 )) 



(3.80) 



where the primes in the U's and V^'s are introduced to discriminate between the 
two single photons. 



We will use the concurrence C(p) [38 



for a X-state like (3.74) is given by (3.33). If we take 



to compute the entanglement, which 



p, 



(3.81) 



we find that ESD appears at a range of values of p that decreases with increasing 
r, in agreement with [68]. Although this would suggest that ESD disappear for 
r large enough, we find that there are high values of p for which ESD exists for 
arbitrary large r. In Fig. 3.14 we represent C(p AB ) in front of x for different 
values of z and p = 0.98. ESD occurs at z/x = fit of the order of 10 7 . Thus, as 
z (that is r) grows, ESD is shifted to higher values of x. It is also interesting to 
analyze the phenomenon of entanglement revival, discovered in these systems 
in [68 . We find that the dark periods [68] between death and revival has 
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Figure 3.16: Negativity N(pbf) = ^(paf) m front of x = r/ct for p = 0.98 and 
z = Qr/c = 2 ■ 10 6 (solid line), 5 • 10 6 (dashed line) and 2 • 10 7 (dotted line). 
Entanglement increases from at x —¥ oo up to a maximum value and then 
decreases and vanishes eventually. 



larger time durations for increasing z. Besides, although in [68] the revival is 
described as a consequence of the photon exchange, for r sufficiently large both 
the ESD and the revival can occur for x > 1, where photon exchange is not 
allowed. We think that the explanation for entanglement revival is closer to 
the spirit of |66[ I72j where entanglement revival between noninteracting atoms 
is interpreted as coming from entanglement transfer between different parts of 
the system. We shall discuss this point in the following sections. In Fig. |3.15| 
we sketch the dependence with p. Although sudden death and revivals appear 
in a very restricted range of the parameter, they are only a particular case of 
the generic behavior of entanglement observed in a wider range, which can be 
described as disentanglement up to a minimum value and growth of quantum 
correlations since then. 



3.3.4 Atom-field entanglement 



Tracing (3.70) over states of atom A (B) the reduced atom-field density matrix 
Pbf {paf) is obtained. Taking the basis {|e0),|el),|e2),|g0),|gl),|g , 2)} ) 
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Figure 3.17: Negativity N(pbf) = ^(paf) (solid line) and concurrence C(pab) 
(dashed line) in front of x = r/ct for p = 0.98 and z — VLrjc = 5 • 10 6 . 
Entanglement atom-atom cycle is clearly correlated with the atom-field cycle, 
although the sum is not a conserved quantity. Although atom-field entanglement 
may change while the other remains zero, both entanglements cannot increase 
or decrease at the same time. 



we have: 



Pbf = Paf 
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Figure 3.18: I concurrence Ca-bf = ^b-af m front of x = r/ct for p = 0.98 
and z = Qr/c = 2 ■ 10 5 (solid line), 5 • 10 5 (dashed line) and 1 • 10 6 (dotted line). 
Entanglement disappears faster than the entanglement between the atoms (Fig. 



3.14) and remains since then. 



with 



p'n = 


| a (1 + A') + (3X"\ 2 , p' 22 = p ' 55 = P22 


P33 = 


\aF' 


+ l3G"\ 2 ,p' 44 = \aX' +j9(l + A") | 2 


Pee = 


\aG' 


+ PF"\ 2 , ^3 = (o(l + A')+l3X")(aF' + f3G" 


P15 = 


(o(l 


+■ A')+l3X")(aU B + PV A )* 


P2i = 


(aU A 


+ PV B )(P(1 + A") + aX 1 )* 


P26 = 


(alf A 


+ pV B ){aG' + PF"Y 


P35 = 


{aF' 


+ PG"){aU B + I3Va)* 


P46 = 


(aX' 


+ 13(1 +A"))(aG' + f3F")* 


N' = 


P'11 + 


P22 + P33 + Paa + P55 + Pm 



(3.83) 



There are no operational generalizations of concurrence for mixed states in 2 x 3 
dimensions like the ones in Eq. (3.82). We will use the negativity [51] N(p), 
which is the absolute value of the sum of the negative eigenvalues of the partial 
transposes of a state p. For the 2x2 and 2x3 cases N(p) > is a necessary 
and sufficient condition for p to be entangled. 

Up to second order in perturbation theory, we have that N' = N and that 
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Figure 3.19: I concurrence Ca-bf — ^b-af in front of x = r/ct for z = 
ftr/c = 2 ■ 10 6 and p = 0.50 (solid line), p = 0.75 (dashed line) and p = 0.98 
(dotted line). Entanglement sudden death occurs for a wider range than the 
entanglement between the atoms (Fig. 3.15). 



the nonzero eigenvalues of the partial transposes of both pbf and paf & re 



p' n +p' 55 ± Vip'n-pUr + W' 
2N' 



and 



A'+ 



P44 



+ p' 5 5 ± V(PU~ P2 2 ) 2 +4|pi 5 | 5 



2N' 



(3.84) 



(3.85) 



if 



being zero the other two. In Eqs. 
second order are retained. Therefore. 
IP15I 2 > P22P44 then X'_ < 0. In Fig. 13.16 we represent N(pbf) = ^(paf) in 



(3.84| and (3.85) only the terms up to 
2 > p' n then A_ < and if 



£24 1 



front of x for same values of p and z of Fig. 1. We see that the negativity grows 
from at x — ¥ 00 [t = 0) to its maximum value and then starts to decrease and 
eventually vanishes, following the opposite cycle to the entanglement of pab- 

Although it would be interesting to look for conservation rules of entangle- 
ment like the ones in [66l [721 E3] , this search is beyond the focus of this paper 
since in our study we are using different entanglement measures in Hilbert spaces 
of different dimensions. Besides, except for the concurrence between atoms A 
and B, the rest of the concurrences in the mentioned papers have not obvious 
counterparts in our case. But it is clear that in general the entanglement cycle 
between atoms is correlated with the entanglement cycle between each atom and 
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Figure 3.20: I concurrence Cf-ab in front of x = r/ct for p = 0.98 and 
z = flr/c = 2 ■ 10 5 (solid line), 5 • 10 5 (dashed line) and 1 • 10 6 (dotted line). 
Entanglement grows from to its maximum value at x ss 0.1 and then decreases. 



the field, as can be seen in Fig. |3.17| in a particular case. Although atom-field 
entanglement may change while the other remains zero, both entanglements 
cannot increase or decrease at the same time. 



3.3.5 Tripartite entanglement 

Tripartite entanglement has been widely studied in terms of the entanglement 
of the different bipartitions A - BC, B - AC, C — AB in the system [8T 1 l82 l 183] . 
where A, B and C stand for the three parties. Here, we will compute the I 
concurrences [84 'Ca-bf, 'Cb-af, ^f-ab, where 



Cj_ 



J-KL 



^2(1 -Trp 2 j), (3.86) 



where J runs form A to F and KL from BF to AB respectively, being pj the 
reduced density matrix of J. A and B stand for the atoms, and F for the field. 



Tracing (3.9| over BF (AF), we find the following density matrices pa (pb)' 



(3 - 87) 

where 

PAU = Pii + p'm + P22 (3.88) 
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Figure 3.21: I concurrence Cf-ab (solid line) and C 
line) in front of x = r/ct for z = flr/c — 2 • 10 5 and p 
cannot increase or decrease at the same time. 



\-BF 

-- 0.98. 



= Cb-af (dashed 
Both magnitudes 



and 



PA22 = P44 + Pee + P22 



and N A = p A n + PA22- In Fig. 



(3.89) 

3.18| we sketch the behavior of Ca-bf and 
Entanglement vanishes before 



'Cb-af in front of x for different values of z 
the death of the entanglement between A and B, and does not have a revival 
Besides, ESD appears in a wider range of p, as can be seen in Fig. |3.19| Now 
tracing (3.70) over AB we obtain the reduced density matrix of the field pf- 



PF 




PFll PF13 

PF22 
P*F13 PF33 



(3.90) 



where 



Pfii = p'n + p'441 PF22 — 2p 2 2, 

PF33 = P33 + ^66> PFQ2 = P13 + P46' 



(3.91) 



In Fig. 



and N F = Pfii + PF22 + PF33 
x for the same values of z and p as in Fig 
to a maximum value at x 



3.20 we represent Cf-ab in front of 
|3.18| Entanglement grows from 
0.1 and then decreases. The growth of Cf-ab 
is correlated with the decrease of Ca-bf and Cb-af in the same way as the 
magnitudes analyzed in the previous section, as can be seen in Fig. 
particular case. 



3.21 for a 
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3.3.6 Conclusions 



We have analyzed in a previously unexplored spacetime region the entanglement 
dynamics of a system consisting in a pair of neutral two-level atoms A and B 
interacting with a common electromagnetic field F . At t — atoms are in 
the Bell state ^Jp \ e e ) + ^1 — p \gg) and the field in the vacuum state. The 
evolution of this state has been considered within the non-Markovian, non-RWA 
approach of quantum electrodynamics up to second order in perturbation theory. 
We find ESD and revival of entanglement in the reduced state of the atoms, 
in a range of p that decreases with the interatomic distance r, in agreement 
with the results obtained with master equations [55]. For r large enough, we 
find that the revival of entanglement can occur with r > ct and therefore 
is not a consequence of photon exchange between the atoms. We find that 
this phenomenon is strongly related to the transfer of entanglement between 
the different subsystems of two parties that coexist in the entire system: we 
obtain sort of entanglement cycle for the atom-field reduced states opposite to 
the atom-atom one. We have considered also the different bi-partitions of the 
system, namely A — BF, B — AF and F — AB, finding similar relationships 
between their entanglement cycles. 
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Chapter 4 



Circuit QED as a quantum 
simulation of 

matter-radiation interaction 
beyond RWA 

"GALILEO: I improved it. LUDOVICO: Yes, sir. I am beginning 
to understand science. [...] Yes, a pretty red. When I saw it first it 
was covered in green." (Bertold Brecht, Galileo.) 

The experimental implementation of the theoretical results introduced in Chap- 
ter 3 is extremely challenging, mostly because of the impossibility of turning 
the interaction on and off at will, among other reasons. In this chapter we will 
move to the framework of circuit QED, which can be understood as 1-D ver- 
sion of matter-radiation interaction with artificial atoms and photons, with the 
advantage of experimental amenability. An important feature of circuit QED is 
that stronger values of the coupling strength can be achieved, entering into the 
so-called "ultrastrong" coupling regime in which effects beyond RWA become 
accessible to experiment. An important consequence of this will be analyzed in 
section 4.1, in which a possible experimental protocol to detect with certainty 
ground state qubit self-excitations is introduced. In section 4.2 we present an 
experiment proposal to test the entanglement generation between superconduct- 
ing qubits coupled to a quantum field initially in the vacuum state, that is the 
circuit QED version of some of the results presented in chapter 3. 
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4.1 Detecting ground state qubit self-excitations 
in circuit QED: 
slow quantum anti-Zeno effect 

4.1.1 Introduction 

The model of a two-level system interacting with one or more harmonic oscil- 
lators can be implemented in circuit QED combining a superconducting qubit 
with a microwave resonator or a transmission line [34 , 35 , 85j . Compared to ex- 
periments in Quantum Optics with microwave cavities [86, [87] or with trapped 
ions |88j . the superconducting circuit experiments have one important advan- 
tage: the strength of the qubit-photon coupling. The fact that superconducting 
resonators and superconducting qubits follow essentially the same physical laws 
makes it possible not only to reach the strong coupling regime [35 ,155], in which 
multiple Rabi oscillations are possible within the decoherence of the cavity or 
the qubit, g ^> k,j, but also entering the ultrastrong coupling regime, g ~ lj, 
in which the internal and interaction energies become similar [891 190] . In this 
new regime the dynamics is very fast and the usual approximations such as the 
Rotating Wave Approximation (RWA) in the Jaynes-Cummings model break 

down E2E1E3]. 

One of the most astounding predictions of the ultrastrong coupling regime is 
that a single qubit can distort its electromagnetic environment, giving rise to a 
ground state in which the qubit is dressed with photons. As we will show in the 
following sections, in the case of a qubit and a single harmonic oscillator, this 
translates into a state which is a superposition of a desexcited qubit and a vac- 
uum, with other states in which the qubit, the oscillator or both are populated 
with excitations and photons, respectively (92] [94| |95] . This is a completely non- 
RWA effect which requires large values of the coupling to be observed. More 
precisely, the excitation probability grows approximately as p e oc (g/u) 2 and g 
has to become comparable to the energies of a photon, huj, or of a qubit, huiQ, 
making the interaction dynamics both very strong and very fast. From the ex- 
perimental point of view it would thus seem unfeasible to probe a physics that 
takes place at speeds of lu ~ 1 — 10 GHz, while the typical measurement appara- 
tus in circuit-QED have response times which are much slower, of about 50 ns. 
There are four routes to escape this problem: making the ultrastrong coupling 
switchable by design [96] . dynamically turning it off by external drivings |97j . 
engineering faster measurement apparatus or looking for new ways to extract 
information out of slow measurement devices. 

In this section we take the slow route, showing that is possible to extract 
valuable information from the fast dynamics of the system with current mea- 
surement technologies. We will study what happens to an ultra-strongly coupled 
qubit-cavity system when the qubit is subject to repeated measurements by a 
detector with a slow repetition rate that is only capable of performing weak 
measurements of the state of the qubit. The main goal is to detect the qubit in 
its excited state starting from the ground state of the system. The first measure- 
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ment has already a small probability of success, as commented in the previous 
paragraph. In case of failure the system is projected to a non-equilibrium state 
which rapidly exhibits a dynamics with an oscillatory probability of excitation, 
mainly due to non-RWA transitions from the ground state of the qubit \g) to 
the excited one |e). By means of performing repeated measurements, we will 
show that the detector is able to probe these usually considered as "virtual" 
excitations of the qubit and the cavity and at the same time reveal information 
of the interaction model. More precisely, the repeated measurements accumu- 
late information exponentially fast and behave like an anti-Zeno effect [98] in 
which the qubit is projected onto its excited state, revealing those ground-state 
excitations that we were looking for. We show that this anti-Zeno "decay" 
\g) — > |e) , is very efficient and does only require a short number of repeated 
measurements with a repetition rate which is much slower than in the standard 
anti-Zeno effect. 

Like other proposals for probing the ultrastrong coupling limit [M] , (see also 
section 4.2), the anti-Zeno dynamics in this work is supported by the counter- 
rotating terms in the qubit-resonator interaction, using as seed the ground state 
excitations of these systems. The phenomenon is absent in the limit of RWA 
in Jaynes-Cummings models. Let us remark that the non-RWA effects are 
being extensively studied not only in the ultrastrong coupling regime of circuit- 
QED but also in other fields like Quantum Optics [§9] • Models of repeated 
measurements on superconducting qubits were considered for instance in Rcf. 
jlOOj and have been implemented in the lab [1011 1102j . 

The structure of this section is as follows. In 4. 1. 2 we will show that the 
eigenstates of the hamiltonian, and in particular the ground state of a qubit- 
cavity system in the ultrastrong coupling regime are not separable, \g, 2n) or 
|e, In + 1) , but linear combinations of these vacua and excitations. More pre- 
cisely, the qubit-resonator ground state contains a contribution of |e, 1) which 
grows with the coupling strength and becomes relevant in the ultrastrong cou- 
pling regime, g ~ to. We will see that after a few ideal periodic projective mea- 
surements of the qubit state, the probability of finding that it is in the state \g) 
tends quickly to 0, even if an uncertainty in the time taken by the measurement 
is considered. In section 4.1.2.5 we will consider a realistic model of measure- 
ment in which large amounts of errors are allowed, showing the robustness of 
our method. Section 4.1.2.6 is devoted to the analysis of the role of relaxation 
and dephasing. We conclude in section 4.1.3 with a summary of our results. 

4.1.2 Detecting ground state qubit self-excitations 
4.1.2.1 The Rabi model 

We will consider the following Hamiltonian, corresponding to a qubit-cavity 
system 

H = H Q + gHi = hioa)a + + hga x (a + a f ), (4.1) 
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where Hloq is the energy splitting between the two levels of the qubit |e) and 
\g) , u) the frequency of the photons in the cavity or resonator field and g the 
coupling strength. 

In the weak and strong coupling regimes, in which the coupling g <C too is 
only compared to the decay rates of the cavity and the qubit, one may treat Hj 
as a small perturbation on top of the bare qubit and resonator states. In this 
limit the counter-rotating terms a^er + , aa~ average out, and the total Hamilto- 
nian becomes equivalent to the Jaynes-Cummings model, whose ground state is 
a separable combination of the qubit ground state and a cavity vacuum, \g, 0) . 




Figure 4.1: (Color online) Probability of excitation for of a qubit p e (blue, solid) 



vs. the dimensionless ratio g/ui, for a qubit-resonator system [Eq. (4.1)] in the 
ground state of a ultrastrong coupled limit oj — ujq = 1 GHz. This line is 
undistinguishable from a quadratic fit (green, dashed). 

In this work we are interested however in the ultrastrong coupling regime, 
in which g approaches the qubit and photon frequencies, lo and loq. In this 
case it is more convenient to look at the state space in the language of parity 
subspaces [M), and treat Hq and Hj on equal footing. Within this picture, the 
Hilbert space splits up in two different chains of states coupled by Hi, and in 
particular the ground state of the system becomes a linear combination of states 
in the even parity sector 

\G) = co \g0) + ci |el) + c 2 \g2) + c 3 |e3) + . . . (4.2) 
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where the coefficients Cj depend on g,(j,u}Q. 



4.1.2.2 Detecting excitations with one measurement 

One of the goals of this paper is design a protocol for measuring the tiny ex- 



citations in the ground states — |ci| 2 + |c 3 | 2 + |c 5 | 2 + . . . in Eq. ( |4.2[ ) — . Let 
us assume for now that we have a good measurement apparatus and that we 
perform a single measurement of the qubit in the ground \G) of the system. 



In Fig. 4.1 we plot the probability of finding the qubit excited after just one 



measurement 

Pe = (Pe) G =■■ (|e> (e|) G (4.3) 

against different values of the coupling strength, assuming always uj — loq and 
g/uj < 1. For the strongest couplings the values of p e are sizable. Moreover, we 
have: 

2 

Pe = A^, ( u = uo, i- <1) (4.4) 

This quadratic behavior comes as no surprise. The main contribution to p e is 
| cx | 2 - If we think of \G) as the free vacuum \g, 0) dressed by the interaction Hi, 
then |ci| may be computed from perturbation theory in interaction picture, the 
leading term being proportional to |(e, l|ii/|<?, 0) | 2 . It is interesting to see how 
these contributions quickly grow as g approaches w, but that at the same time 
the signal in current experiments with 10% coupling strengths, might have a 
too small excitation signal to be accurately detected. 

This work is born from the idea that perfect projective measurements in c- 
QED might be too difficult, as existing measurement apparatus may be too slow 
or not have enough sensitivity to capture those excitations. The constraint of 
time is found, for instance, in flux qubit measurement devices based on SQUIDs, 
which roughly work as follows: A very short current pulse is sent to the SQUID, 
instantaneously changing its potential from a periodic function to a washboard 
potential. In this brief period of time, one of the flux qubit states which is sitting 
inside the SQUID may provide, through its intrinsic current and flux, enough 
additional energy for the SQUID to tunnel into a voltage state. This stochastic 
process is random in time and does no have a 100% success rate. Moreover, it 
requires an additional sustained current that keeps the SQUID in that voltage 
state during an integration time large enough for the electronics to realize that 
the measurement succeeded. Adding the excitation and integration phases, the 
best experimental setups bring the detection time down to tens of nanoseconds, 
which is still slower than the qubit-resonator dynamics -1.6 ns for a 600 MHz 
coupling, and much shorter for the qubit and resonator periods, 1/ui. 

An additional complication of the ultrastrong coupling limit is that an ar- 
bitrary measurement device might not have enough good coupling to either the 
qubit or the resonator in the ultrastrong coupling regime. If we assume that 
both quantum systems interact so strongly that their eigenstates are highly 
entangled states with large energy gaps, ~ g,u),Uo, the detector could have 
problems coupling to those states and breaking their energy level structure. In 
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other words, the measurement device couples through an operator, a z , which 
typically represents a perturbation of the qubit-resonator model, and if that per- 



turbation, which aims at breaking the linear combinations (4.2 1, is not strong 



enough, it might not extract any information from the system, or the amount 
of information might be reduced, becoming an off-resonant, weak dispersive 
measurement. 

All these considerations brought us to the idea of using more than one mea- 
surement steps in the same experiment, with the aim of increasing the amount 
of information that it is extracted from the same state. This can be done be- 
cause the kind of measurements done in experiments are non-destructive: the 
same qubit can be continued to be measured at another time. It is true, how- 
ever, that the interval between measurements might carry a strong, fast and 
almost chaotic dynamics [Fig. 4.2 , arising from the fact that the measurement 



brings the system into a non-equilibrium state, even if it did not produce any 
information. We will show that this is not a limitation, but a plus, and that the 
repeated measurements may characterize the intermediate dynamics. 




Figure 4.2: (Color online) After the qubit has been measured once, the qubit in 
the qubit-resonator system is left in a non-equilibrium state, |^) . Here we plot 
the probability of excitation for the qubit p\ e as a function of the dimensionless 
time wt, shortly after that measurement. We show three situations, to = luq = 1 
GHz and g/ui = 1/3 (dashed), g/ui = 2/3 (solid) and g/ui = 1 (dotted), which 
exhibit fast dynamics. 
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4.1.2.3 Repeated measurements: survival probability 

If we measure the qubit once, the measurement apparatus does not click and 
we are working with a perfect projective measurement, we conclude that the 
qubit-resonator system has been projected onto the state 

I*) =]Tc 2 „| 5 ,2n), (4.5) 

n 

which is a (normalized) linear combination of deexcited qubits and some photons 
in the cavity. By measuring the ground state in an improper basis, we have 
created a non-stationary state that will evolve very quickly, with frequencies 
that are close to g,w and ojq- Lacking any other relaxation mechanism than 
the cavity and qubit decoherence times, these oscillations will be sustained for 
a large period of time, causing the qubit to get reexcited multiple times. The 
excitation probability 

pi e (t,0) = (*(t) | A| *(*)), (4-6) 
may be computed from the initially measured state as 

|*(t)> oc e- im (l - P e ) \G) . (4.7) 



As Fig. 4.2 shows, p\ e exhibits very fast oscillations, but also average to a 
nonzero value, which is always close to the ground state excitation probability of 
the qubit, p e = J2 n \ c 2n+i\ 2 - Consequently, if we perform a second measurement 
at a later time t\ we will have again a certain probability of success pi e (ti,0) 
of detecting the state |e) , and a certain probability of failure p\ g (ti,0) = 1 — 
Pi e (ti,0). In the latter case the system is projected to a new state with a new 
time dependent probability P2e(t2,ti), and so on. After a few measurements 
we can define the survival probability as the probability that we have never 
detected a state |e) in the qubit 

Pg = Pgplg(tl,0)p 2 g(t2,h)...p N g(t N ,t N -l). (4.8) 

A key idea in the interpretation of this formula is the fact that the intermediate 
probabilities p ng are on average very similar, and almost independent of the 
timespan among measurements. For the range of couplings that are within 
intermediate reach in experiments, g/u> ~ 0.1 — 1, we have verified numerically 
and perturbatively that this probability is well approximated by a quadratic 
law 

g 2 

Png^l-Xn^ (4.9) 

with minor differences among realizations, x„. The accumulation of products in 



Eq. (4.8 1 leads to an approximately exponential decrease of the survival proba- 
bility 

*? -fli 1 - Xn£) ~ cxp (-Nx£) , (4.10) 



71=1 
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as long as << 1- This exponential behavior is typical of the so called anti- 
Zeno effect, in which repeated measurements of a quantum system accelerate 
the transition of a quantum system between two states. In our case the repeated 
measurements are rather creating a non-unitary evolution that excites the qubit 
from \g) to |e) using as seed the nonzero excitation probability p\ e = J2 n \ c 2n+i | 2 
which is present in the equilibrium state of the qubit-resonator system. This last 
point is particularly important because this anti-Zeno evolution is impossible 
when the ground state of the qubit and the resonator is the vacuum |<7,0) . In 
this case g/u is so small, and p\ e so close to zero, that all measurements will 
give no signal at all and the qubit will remain in the state \g) for the duration of 
the experiment. As we will see in the following, there is a key difference between 
the effects described in this section and the standard anti-Zeno effect: we need 
only a few measurements and they can be widely spaced in time. 

In the following sections we will summarize extensive numerical studies of 
the anti-Zeno dynamics. We have contrasted these with various semi-analytical 
methods, one of which, the use of truncated Hilbert spaces, helps us in un- 
derstanding the reason for this behavior. For the range of couplings of current 
interest, g/u> ~ 0.1 — 1, it suffices to take two photons, and the ground, \G) , plus 
the two excited states \E'), \E") within the same parity subspace. All states 



can be expanded as in Eq. (4.2 1 with coefficients Cj, c' i7 c"., as linear combina- 
tions of \g0), |el), \g2). After the first measurement, the qubit will end up in 
an excited state with probability |ci| 2 and it will remain in the unexcited state 
with | Co | 2 — 1 — | ci | 2 , ending up in a combination 

P g \G) = co\G)+c' \E')+c%\E") + ... (4.H) 

The crudest approximation would be to neglect all excited state contributions 
and assume that after each measurement, either the state |e) is detected, or the 
system ends up in \G) . In this case the survival probability would be exactly 
exponential 

p? =(i-\*\ 2 )f[\<*\* = \<>o\ aN+ '- (4.i2) 

i=i 

In practice, however, the combined system does not end up only on the ground 
state, but gets excited state contributions from \E') , \E") . When we average 
the contributions over the period in which the measurement takes place, we find 
that already after the first measurement step, the excited states add up to the 
total probability, (pi e ) T = \c\\ 2 |c | 2 + \c[\ 2 |c | 2 + |c 1 1 2 |c | 2 + . . . , enhancing 
the original behavior. 

4.1.2.4 Numerical experiments 



We have verified the anti-Zeno dynamics and the exponential law (4.10) by 



means of exact numerical simulations in which we compute the outcome of 



repeated measurements on a qubit-resonator Dicke model (4.1). We will now 
explain the main results of this study. 




Figure 4.3: Survival probability after eight measurements Pg vs. g/cJ, with 
uj = ujq. The measurements are performed with periods Ti, T 2 — V^Ti and 
averaged over 100 values of T\ within the interval 27r[0.1, 5]. Note how the law 
approximates the Gaussian behavior in Eq. (4.10). 



From an experimental point of view it might be interesting to maximize the 
exponent x, optimizing the measurement repetition rate to hit all the maxima 



in the evolution of the excitation probability [See Fig. 4.2 . However we found 
that this is very difficult and demands a lot of precision on the measurement 
apparatus; for small errors or some measurement randomization this procedure 
drives the apparatus into exactly the opposite regime: always hitting the minima 
of excitation. Seeking a more robust, less demanding approach we opted for 
using two incommensurate periods, T\ and T 2 ~ y2T\ , simulating measurement 
at times 

t n € {T u Ti +T 2 ,2T 1 +T 2 ,2T 1 + 2T 2 ,...}, (4.13) 

and at most optimizing the value of Ti . 

With this approach, and exploring different values of Ti, we have studied 
the survival probability and concluded that the exponential laws are really ac- 



curate. As shown in Fig. 4.3 if we fix the total number of measurements to 
be N — 8 and sample various periods, Ti, we recover on average the Gaussian 
behavior exp(— Nxg 2 /<jJ 2 ) deduced in Eq. (4.10). Instead of fixing the number 



of measurements, we can also study the same law and verify the exponential 
decay with respect to N. This is shown in Figs. [Oh-b, where we plot the ac- 



cumulated survival probability, P^ , as a function of time, and fit it against the 
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Figure 4.4: (a) Probability after the n— th single measurement, p ng , (dashed) 
and accumulated survival probability, P^ = Yln=i P n a ( souc vs. dimensionless 
time urt. We use g/ui = 1 and perform measurements with approximate periods 
ujTi = WT2/V2 = 2ir, averaging over random perturbations of the actual mea- 
surement time, t n , within the interval ujt n + [— 0.27T, +0.27r]. (b) Survival proba- 
bility P^ (solid lines) and the corresponding exponential fits (dashed lines) for 
g/uj = 1/3 (circles), 2/3 (crosses) and 1 (squares), (c) Mean value p g (solid) and 
the corresponding quadratic fit (dashed) vs. dimensionless coupling strength ^. 
All plots assume u) = ujq ari d (b,c) use ujT\ — 3n/ 4 



same exponential (4.10). 

It is important to remark that the exponential decay is a robust signature 
that survives even when the measurement does not take place at precise times, 
t„ , from the list given before (4.13). This has been verified by simulating mul- 
tiple runs in which t n is randomly perturbed around its average value, and 
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computing the survival probability. We want to remind the reader the impor- 
tance of this robustness, because some measurement apparatus such as SQUIDs 
behave stochastically and produce a signal at a random time that can not be 
determined a-priori. The fact that the measurement protocol and the resulting 
physical behavior are independent of a precise control is encouraging. 



The accuracy of the exponential law (4.101 suggests that the survival proba- 



bility of a single measurement remains constant throughout a single experiment, 
p g ~ p ng Vn. This is qualitatively confirmed by Fig. 4.4 1, where we show that 
these values oscillate around a mean one that is close to the average population 
of \g) in the ground state, i. e. to p g . This suggests us to consider average values 
and approximate 

2 

EPng -, - 9 



(4.14) 



which has the expected quadratic behavior. This estimate is confirmed by 



Fig. 4.4;, where the quadratic fit is almost undistinguishable from the actual 
behavior. The final question which remains to be answered is whether the expo- 




Figure 4.5: Survival probability Pg vs. t/Ti, with g = lu = loq GHz and 
uTi — 7r (blue, circles), 2tt (red, crosses) and 37r (green, squares). Each marker 
corresponds to a measurement. 



nent x depends on the frequency of the measurements or not. For that we have 
fixed the coupling strength and explored three values of the period, T\ , studying 
the average exponential behavior. The result is shown in Fig. |4.5[ collapsing all 
numerical simulations in the dimensionless quantity t/T\, and finding that they 
have very similar slopes. 
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4.1.2.5 Weak measurements 

So far we have considered ideal projective measurements, introducing only some 
stochasticity in the time at which the measurement event is produced. We will 
now add another ingredient to our measurement model, which is the possibil- 
ity that the detector only performs a partial measurement, leaving the state 
"untouched" with a nonzero probability, e. 

We can easily model an imperfect detector using the formalism of completely 
positive maps, operations that transform density matrices into density matrices. 
If p and p' are the states of the qubit-resonator system before and after the 
measurement, we will write, up to normalization 

p' = (l-e)(l-P e )p(l-P e )+ep. (4.15) 

This is read as follows. With probability e the measurement device will do 
nothing, leaving the state untouched. With probability (1 — e) the measurement 
device will detect the state of the qubit. In this case it will either give us a 
positive signal, moment at which we will stop the experiment, or it will not 
produce anything at all, and we will continue with the projected state (1 — 
P e )p(l — P e ), that has the qubit deexcited, \g) . 

This qualitative model describes measurements from a SQUID |101[ 1102] . 
where we place ourselves on the verge of metastability and assume that if the 
qubit is in the excited state, |e) , the SQUID will tunnel to the voltage state 
with probability (1 — e), giving no signal for \g) . Note that with probability 
e the SQUID may not tunnel and then we will gain no information about the 
qubit or the resonator. 

In Fig. |4.6| we analyze the impact of e in our previous results. Even for large 
errors e = 0.2 we retain the exponential behavior observed in Fig. |4.4fr , with 
acceptable error bars that decrease with increasing number of measurements — 
in other words, the qubit is still efficiently projected to the excited state. 



4.1.2.6 Relaxation and dephasing 

Throughout this work we have considered in the numerical simulations the 



model given by the Hamiltonian in Eq. (4.1) which do not include effects like 
relaxation or qubit dephasing, usually included in master equation approaches. 

We want to remark that it is still an open question, both experimentally 
and theoretically, to understand and model the dissipation and decoherence 
processes of quantum circuits in the presence of ultrastrong qubit-cavity cou- 
pling. One popular approach |1031 1104] is to combine the usual photon leakage 
mechanism from quantum optics models, C(p) ~ 2apa) — a) ap — pa^a, with 
the qubit-cavity Hamiltonian. Note that in such a combination, the asymptotic 
states of the dissipation (the vacuum) and of the interaction (populated cavity) 
are incompatible, and one may find excitations induced by the dissipative terms, 
an infinite stream of photons leaking out of the cavity and other controversial 
phenomena. 
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Figure 4.6: (Color online) Survival probability (solid lines) and p ng (dashed 
lines) vs. u>t, with g = lu = loq = 1 GHz , ujT\ = 2?t and e = 0.2,0.1 and 
(crosses, squares and circles, respectively). The probability of a measurement 
at a time uit = ujto is averaged over 20 random values within the interval [uit — 
0.27r,u;i + 0.2tt] 



These effects disappear when one rederives the master equation from first 
principles, using the qubit-cavity eigenstates of the ultrastrong coupling model 
and the usual zero temperature baths. In the resulting models the main re- 
laxation mechanisms are found to be the decay to the ground state \G) and 
a dephasing of the joint cavity-qubit states — in other words, dissipation and 
decoherence in the proper basis — . If we assume this reasonable model, then 
we can conclude that the exponential laws derived in this manuscript are not 
significantly distorted. To begin with, relaxation to the ground state \G) just 
makes the experiment closer to the truncated Hilbert space model considered 
in Sect. 4.1.2.3 and in particular to the exponential law from Eq. (4.12). For 
strong couplings, decoherence amounts to random modulations of the qubit- 
cavity energy levels, without significantly affecting the populations, \ci\ 2 . Since 
this is the most relevant quantity in all the previous discussions, we can also 
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expect that, up to minor changes in the rates, the anti-Zeno effect will also 
survive. 

4.1.3 Conclusions 

We have considered a system consisting in a superconducting qubit coupled to 
a closed transmission line, operating in the ultrastrong coupling regime. The 
ground state in such scheme is not just a product of the ground states of the 
qubit and the cavity, as is the case for weaker couplings. On the contrary, 
the vacuum of the system is dressed by the interaction and so it contains a 
relevant probability of finding the qubit excited. This probability is proportional 
to the square of the coupling strength. We have introduced a protocol for 
detecting that excitations with certainty, maximizing the small probabilities 
that are obtained with only one measurement. 

Our main result is that, after a number of periodic measurements of the 
qubit, the probability of finding it in the ground state in all the measurements 
goes exponentially to zero, even if the measurements are weak and are per- 
formed with a slow repetition rate in comparison with the fast dynamics of the 
interaction. We refer to this as slow quantum anti-Zeno effect. Like the well 
known quantum anti-Zeno effect, the result is the acceleration of a transition, in 
this case the exotic transition \g) — > \e) , which becomes relevant in this regime 
due to the breakdown of the RWA. But this procedure is less experimentally 
demanding, since it requires a smaller number of measurements and a shorter 
duration of the period at which they are performed. We have shown that the 
protocol is robust to large errors in the measurement process, when a realistic 
SQUID readout is considered. 

This is one of the first experimentally accessible consequences of the new 
ultrastrong coupling regime and can only be derived beyond the RWA. The 
physical nature of the ground state qubit self-excitations, commonly considered 
as a virtual process without possible experimental record, seems now to be 
clear. Moreover, although the ultrastrong coupling entails a very fast dynamics, 
we have shown that valuable information of the interaction can be extracted 
efficiently with the current slow and imperfect measurement technologies. 

Finally, we want to remark that strong qubit excitations have also been found 
theoretically in models that combine the full Rabi coupling with traditional 
dissipative contributions [10311104] . However, the form of those dissipative terms 
is questionable in non-RWA setups, and furthermore, there is no justification 
to equate the sparse measurement setup in this work to a particular dissipative 
model. This lack of equivalence between models manifests in the fact that, as 
we have seen numerically, the sparsely repeated measurements can hit certain 
resonances that invalidate the anti-Zeno dynamics. 
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4.2 Entanglement dynamics via propagating mi- 
crowave photons 

4.2.1 Introduction 

Quantum mechanics does not allow us in general to consider two arbitrary dis- 
tant systems as separate [2] . In some cases there exist quantum correlations that 
cannot be generated by local operations and classical communication between 
remote systems. Time enters this picture through two different questions. The 
first one is related to the speed bound of a hypothetical superluminal influence 
which could explain all quantum correlations, estimated to be 10 4 c in a recent 
experiment [105] . The second question is of a more practical nature inside the 
quantum theory [HI H21 HH] (and Chapter 3 of this Thesis) : what is the speed 
at which two distant systems become entangled? 

Quantum field theory (QFT) fulfills the principle of microscopic causality by 
which two space-like separated events cannot influence each other |106j and thus 
cannot be used to transfer information [46J [14] . We may then ask whether mi- 
crocausality also sets a limit on the speed at which entanglement can be created 
between two separate systems. More precisely, can two subsystems, supported 
at regions (x, t) and (x',i'), become entangled while they are still space-like 
separated? Or in simple terms, can finite quantum correlations develop before 
signals arrive? As we have analyzed in Chapter 3, the answer to this far reaching 
question is yes, it is possible. After all, Feynman propagators are finite beyond 
the light cone and even before photon arrival there exist correlations between 
the vacuum fluctuations at any two space-like separated events. 

In this section we demonstrate that circuit QED is arguably one of the most 
suitable fields to study the dynamics of entanglement between distant systems. 
One reason is the existence of various choices of high quality superconducting 
qubits, the so-called artificial atoms [WJ 108, 1091HTD]. Another reason is the 
possibility of coupling those qubits strongly with traveling photons using mi- 
crowave guides and cavities [3H E3 ■ Furthermore, those coupling strengths 
can reach the ultrastrong coupling regime (1111 I112[ I113[ 189"! 190] . where the 
qubit-photon interaction approaches the energies of the qubit and photons. In 
this case, the rotating-wave approximation (RWA) breaks down and a different 
physical structure emerges. Such regimes can be activated and deactivated [96] ■ 
facilitating the creation of a fairly large amount of entanglement in a time- 
dependent way, as we will see in this work. 

We will discuss some of the preceding questions in the framework of a precise 
circuit QED setup, see Fig. |4.7[ consisting on two well separated superconduct- 



ing qubits coupled ultra-strongly to an open transmission line [Fig. 4.7 d]. The 
waveguide provides a continuum of microwave photons propagating with uni- 
form velocity, v, mediating an interaction between the qubits. Given an intitial 
separable state in which only qubit A is excited, we have studied the evolution 
of correlations and related it to the propagation of photons between qubits. The 



main results are: (i) Outside the light cone, that is in region I of Fig. 4.7 1 where 
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Figure 4.7: (a) Qubits that interact via traveling photons with finite velocity v 
can be space-like (I, white) or time-like (II, shaded) separated, depending on the 
value of £ = vt/r. While only in II they are causally connected, entanglement 
may appear already in region I. (b) A possible implementation of these ideas 
consists of flux qubits ultra-strongly coupled to a common transmission line, 
(c) With a slight modification, the coupling of the qubits to the line can be 
dynamically tuned via fast magnetic fluxes, $ (Color online). 



£ = vt/r < 1, the excitation probability of qubit B is independent of the dis- 
tance r to qubit A. (ii) Still in region I, entanglement between the qubits always 
takes a finite value and grows with time, (iii) Once the qubits are time-like sep- 
arated, that is as soon as we cross into region II, entanglement grows faster than 
the excitation probability of qubit B and takes sizeable values. Result (i) is a 
manifestation of the fact that our Quantum Field Theory (QFT) model satisfies 



66 



microscopic causality, which formally translates into the vanishing of commuta- 
tors associated with observables at space-like separations, [Q(x, t), Q'(x', t')} = 
for |x — x'| 2 — c 2 (i — t') 2 > 0. Furthermore, it shows that two qubits which are 
space-like separated cannot be used to communicate superluminal information. 
Result (ii), on the other hand, reveals the fact that correlations between vacuum 
fluctuations at separate points can be established at arbitrarily short times, even 
though they are non-signalling and cannot transmit information. 

It is important to remark that the previous questions have been posed the- 
oretically using model detectors [9] , two-level atoms [10] (and Chapter 3 of this 
Thesis), scalar fields [9J[T2] and photons [TU] (and Chapter 3 of this Thesis), yet 
no experimental test has been accomplished. However in this work, we show 
that the access to the ultrastrong couplings in circuit QED allows us to explore 
these ideas with very advantageous parameter ranges. 



4.2.2 Superconducting qubits coupled to a quantum field 

Our setup consists of two qubits, A and B, interacting via a quantum electro- 
magnetic field. The qubits have two stationary states |e) and \g) separated by 
an energy Ml and interact with a one-dimensional field, which propagates along 
the line connecting them, 

V(x) = i J dk^fm7 k [e lkx a k + H.c.] . (4.16) 

This field is described by a continuum of Fock operators 

[a k ,a\,] = 5 k k', (4-17) 

and a linear spectrum, 

u k =v\k\, (4.18) 

where v is the propagation velocity of the field and plays the role of the speed of 
light. The normalization and the speed of photons depend on the microscopic 
details. In particular 

v = l/Vcl, (4.19) 

where c and I are the capacitance and inductance per unit length. 

We consider qubits that are much smaller than the relevant wavelengths, 
A = v/il, and lay well separated. Under these conditions we can split the 
Hamiltonian, H = H + Hj, into a free part for the qubits and the field 

H = ^Ml(a% +ag) + J dkfno(k)a\a k , (4.20) 

and a point-like interaction between them 

R i = - E t'Vfa)- ( 4 - 21 ) 

J=A,B 
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Here xa and xb are the fixed positions of the atoms, and 

dj = dx a x j (4.22) 

is equivalent to the dipolc moment in the case of atoms interacting with the 
electromagnetic field. 

In what follows we choose the initial state 

|V(0)) = |e<?0>, (4.23) 

where only qubit A has been excited, while both B and the field remain in their 
ground and vacuum states, respectively. In the interaction picture given by the 
"free" Hamiltonian H , the system evolves during a lapse of time t into the state 

= T[e- l ti dt ' Hl(t "> /h ] \eg) ® |0> , (4.24) 

T being the time ordering operator. Up to second order in perturbation theory 
the final state can be written as 

m))= [(1 + A)\eg)+X\ge)]®\0) + 
(U A \gg)+V B \ee))®\l) 

+(F \eg) + G \ge)) ® |2) + 0(d 3 ). (4.25) 

The coefficients for the vacuum, single-photon, and two-photon states, are com- 
puted using the action (J = A, B) 



s+- 1 



( e im ' (ej\da x j\gj) V(xj, t')dt' = -(SJ^ (4.26) 
Jo 

among different photon number states \n) , n = 0, 1, 2 . . ., being 

\n) (n\ = — t [ dfei.... / dk n \h...k n ) (h...k n \ (4.27) 



ni 
and 

\k)=a{\0). (4.28) 
Only one term corresponds to interaction 

X = (0\T(S+S^)\0). (4.29) 

This includes photon exchange only inside the light cone, vt > r, and vacuum 
fluctuations for all values of t and r, being 

r = x B -x A (4.30) 

the distance between the qubits. The remaining terms are 
1 
2 



A = -{Q\T{S\S- A +S B S+)\Q) (4.31) 



U A - (l|5^|0),Vb = <l|5+|0), 
F = ±(2\T(S+S2+SbS+)\0),G = (2\T{S+SX)\0). 
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Here, A describes intra-qubit radiative corrections, while Ua,Vb,F and G cor- 
respond to single-photon emission events by one or more qubits. 

All the above can be understood as the 1-D c-QED version of the formalism 



in chapter 3 The coefficients in Eq. (4.25) can be computed analytically (Ap 



pendix B) as a function of two dimensionless parameters, £ and K. The first 
one, 

£ = vft/r, (4.32) 
was introduced before and it distinguishes the two different spacetime regions 



[Fig. 4.7 1], before and after photons can be exchanged. The second parameter 
is a dimensionless coupling strength 



K 



Ad 2 N 

h 2 v 



(4.33) 



Note that the qubit-line coupling 



9 



(4.34) 



corresponds to the qubit-cavity coupling that appears by taking the same trans- 
mission line and cutting to have a length L = A thus creating a resonator 
Refs. [35l|34]. This formulation has the advantage of being valid both for induc- 
tive and capacitive coupling, the details being hidden in the actual expressions 
for d and N. 

Tracing over the states of the field, we arrive at the following reduced density 
matrix 



Pab = 



( Pn 



V P*u 



p u \_ 

P22 P23 

P*23 P33 

p 44 J 



(4.35) 



representing the two-qubit state in the basis formed by |ee) , \eg) , \ge) , and 
\gg) . The coefficients with the leading order of neglected contributions are 



Pn 

P33 

Pli 

P23 



\v\l 

IX' 2 



0(d 4 ), p 22 
\X{ 2 + \G\ 2 + (D(d 6 ), Pi4 = \U\ A + 0(d 4 ) 
U* A V B +0(d 4 ) : 
X* +0{d i ), 



1 + 2Re(A) + 0(d 4 



\S+S+\0)+O(d 4 ) 



(4.36) 



and the state is normalized c = pa. 

Let us now remark the validity of the perturbative methods applied in this 
work. 



The leading corrections to C(pab) (see (3.33)) come from the leading 
order corrections to P23, Pn, P<a (Eq. ( 4.36| )). In the case of P23 we have 



P2 3 {d 4 ) = (1 + A)X* + FG* 



(4.37) 
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and 

p 23 (d 6 ) = P2 3 (d 4 ) + X 1 + X 2 , (4.38) 

where X\ comes from the interference of one and two photon exchange ampli- 
tudes and X2 comes from the probability amplitude of three photon exchange. 
A rough upper bound for these two terms is given by 2 \X\ 3 . For p\\ and ^44 
they involve a number of photon emissions and re-absorptions by the same atom 
or by the other, giving a term 

PuPaM") = \Ua\ 2 \Vb? + A, + A 2 , (4.39) 

where rough upper bounds to A x and A 2 sue 2\A\\U A \ 2 \V B \ 2 and 2|X||[/ yl | 2 |y B | 2 
respectively . All these products are shown to be small for the regions of interest 
discussed here, £ < 2. The same techniques can be extended to all orders in per- 
turbation theory since the bounds to the different contributions can be grouped 
and treated as power series, giving rise to corrections that remain negligible as 
long as \A\, \X\, \Ua\ 2 and |Vb| 2 are small enough, like in the parameter range 
explored in this work. Finally, note that similar calculations and results can be 
obtained in the case in which the qubits have close but different frequencies. 



4.2.3 Entanglement dynamics and single photons 



We will use the concurrence C to compute the entanglement of the X-state in 
(4.351 which is given by (3.33). 



Since all quantities depend only on two dimcnsionlcss numbers, £ and K, 
we can perform a rather exhaustive study of the dynamics of entanglement 
between both qubits. To cover the widest possible spectrum of experiments, we 



have chosen coupling strengths (4.331 over three orders of magnitude 



A' 

Kr 



= 1,10,100,1000. 



(4.40) 



The smallest value 

A = 1.5-10~ 4 , (4.41) 

which corresponds to g/ir ~ 175 MHz and J7/27T ~ 10 GHz, that is for instance 
a charge qubit in the strong coupling limit with a transmission line |34j . The 
largest value, K = 1000 K corresponds to g ~ 2tt x 500 MHz and fl ~ 2tt x 2 
GHz, and typically corresponds to a flux qubit directly coupled to a transmission 



line [1121 155] . as shown in Fig. 4.7 >c 



In Fig. 4.8 we plot the value of the concurrence for two qubits which are 
separated a distance r = A/8, using the couplings discussed before. Note how 
the entanglement jumps discontinuously to a measurable value right inside the 
light cone (£ > 1), signaling the arrival of photons. Furthermore, even a certain 
amount of entanglement appears outside the light cone, before photons could be 
exchanged. This is best seen for the largest couplings, as Fig. |4.8b illustrates. 

The dynamics looks even more exciting when we go back to lab time and 
space. Fig. |4.9| shows the concurrence and the excitation probability of qubit B, 

PB = \V B \ 2 /c + 0(d 4 ), (4.42) 



70 



0.25 



0.15 




0.05 



Figure 4.8: a) Concurrence vs. dimensionless separation £ for r = irv/Afl ~ A/8 
and couplings X = K , 10K , 100K and 1000-Ko (bottom to top) b) Zoom 
around £ = 1 for the strongest coupling = IOOO-Kq- (Color online). 



for two different separations, r = A/12 and r = A/8. The probability of excita- 
tion appears as independent of the qubit separation. This is exactly the case 
for the lowest order considered here, which only accounts for B self-interaction, 
and at all orders in perturbation theory |46j outside the light cone of this setup 
(region I in Fig. 4.7 1). This is in full agreement with microcausality. However, 



as can be seen in Fig. |4.9[ what was a tiny concurrence jumps to a sizable value 
when crossing the light cone fit = 2n/\2 and fit = 27r/8. In other words, from 
the experimental point of view, it is the entanglement between the qubits and 
not the excitation probability pb what best signals the presence of a light cone 
and a finite propagation speed. 



4.2.4 Experimental implementation 

In order to study the dynamics of quantum correlations between the two su- 
perconducting qubits, one has to perform a partial or full tomography of their 
state. In the first and simpler case, performing measurements in different basis 
should be enough to gather an entanglement quantifier, such as a Bell inequality 
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Figure 4.9: Concurrence (dash) and probability of excitation of atom B (solid) 
vs. dimensionless time, fit. Qubits are separated by r = A/12 (circles) and 
A/8 (crosses) and have a coupling strength K — lOOOiio- Note that, following 
microcausality, the excitation probabilities do not depend on the separation r 
outside the light cone (Color online). 



violation or, as studied in this paper, the concurrence. This has to be repeated 
many times, not only to gather sufficient statistics, but also to resolve different 
of instants of time before and after the light-cone boundary. This may seem a 
daunting task, but thanks to the speed at which quantum circuits operate and 
their fast repetition rate, it will be as demanding as recent experiments realizing 
a controlled-NOT gate [114] or full two-qubit tomography jl 15] . 

The actual experimental challenge, though, arises from the need to perform 
quantum measurements of the qubit state and ensuring that this state is not 
altered by the ongoing dynamics. One possibility is to perform very fast mea- 
surements of the qubits, which means faster than I/O. The typical response of 
measurement apparatus, which in the case of SQUlDs is around a few nanosec- 
onds, sets an upper limit on the qubit and photon frequencies of a few hundreds 
of megahertzs, though we expect this to be improved in the near future. 

Another more reliable approach is to connect and disconnect the coupling 
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between the qubit and the transmission line. In this manner, we could prepare, 
entangle, and finally measure the qubits without interference or decay processes. 
If we work with flux qubits, a simple approach is to apply a very large magnetic 
flux on both qubits, taking the qubit away from its symmetry point. From a 
mathematical point of view, this amounts to adding a large contribution E<j x A B 
to the Hamiltonian. If done very quickly, the field projects the qubit on the same 
basis on which the coupling operates, eliminating the possibility of spontaneous 
emission. One would still need to combine the switching of this flux with short 
pulses that rotate the qubit basis in order to perform a complete set of mea- 
surements. The last and most elegant possibility is to effectively switch off all 
couplings between the qubit and the surrounding field. This can be achieved 
using a direct coupling between the qubit and the transmission line, with an 
scheme that incorporates an intermediate loop [Fig. 4.7;]. As Peropadre et al. 
show in a recent work [96], the result is a coupling that can be rotated and 
completely deactivated in a time of about 0.1 ns, that is the time needed to 
inject flux through the loop. The advantage is that, contrary to the case of a 
large external flux, the influence of the line is completely suppressed and makes 
it possible to easily rotate the qubits to perform all needed measurements. 



4.2.5 Conclusions 

Summing up, in this work we have proposed a circuit-QED experiment to study 
the dynamics of entanglement between two qubits that interact by exchanging 
traveling photons. Our work focuses on the existence of a finite propagation 
speed, the appearance of a light cone, the notion of microcausality and the pos- 
sibility of achieving entanglement both by means of the correlated fluctuations 
of the vacuum and by photon exchange. The resulting predictions have a wide 
interest that goes beyond the assessment of microcausality in the QED of quan- 
tum circuits, demonstrating that the open transmission line is a useful mediator 
of entanglement, much like cavities and zero-dimensional resonators. Further- 
more, the experiment we propose is also among the simplest ones that can probe 
the effective QFT for waveguides, both asserting the existence of propagating 
single photons and probing the dispersion relation at the single-photon level. 
Finally, we have shown that entanglement via traveling photons works better 
for stronger qubit-line couplings, making it one of the first potential applications 
of the ultrastrong coupling regime [5H1 [90] ■ 
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Chapter 5 

Causality in 

matter-radiation 

interactions 



"She was particularly exasperated by the behavior of subatomic par- 
ticles. She wanted the universe to behave sensibly" (Martin Amis, 
The pregnant widow) 

Throughout this thesis we have analyzed several features of nonlocality in Quan- 
tum Mechanics, like the generation and destruction of entanglement between 
qubits outside their mutual light cone. We have remarked many times that 
these effects do not represent a violation of causality of any kind. In section 
5. 1 of this chapter we will show explicitly that this is indeed the case, in the 
context of the so-called Fermi problem and the long-lasting controversy asso- 
ciated to it. Although the results are valid also for the 3-D interaction of real 
atoms and the electromagnetic field, we choose the theoretical framework of 
1-D circuit QED in order to make a realistic experimental proposal. Since this 
proposal involves measurements of a qubit excitation probability at short times 
in a regime at which non-RWA contributions are relevant, in section 5.2 we deal 
with the effect of these contributions in the interpretation of the readouts. 

5.1 The Fermi problem with artificial atoms in 
circuit QED 

5.1.1 Introduction 

Information cannot travel faster than light. But in quantum theory, as we 
have seen in chapters 3 and 4, correlations may be established between space- 
like separated events. We remark again that these facts are not contradictory, 
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since correlations need to be assisted with classical communication in order to 
transmit information. 

The two physical phenomena above arise in a natural fashion in the following 
situation, which is the so-called Fermi problem [15 , originally proposed by Fermi 
to check causality at a microscopic level. At t — a two-level neutral atom A 
is in its excited state and a two-level neutral atom B in its ground state, with 
no photons present. If A and B are separated by a distance r and v is the 
speed of light, can A excite B at times t < r/vl Fermi 's answer was negative 
but his argument had a mathematical flaw. When a proper analysis is carried 
on, fundamental quantum theory questions arise due to the interplay between 
causal signaling and quantum non-local phenomena. 

These issues led to a controversy [HI [TTJ 11161 HB] on the causal behavior of 
the excitation probability of qubit B, whose conclusions were never put to ex- 
perimental test. A notorious claim on causality problems in Fermi's two-atom 
system was given in [IB]. The reply of [T7] was in the abstract language of 
algebraic field theory and the proof of strict causality of [JB] is perturbative, 
although given to all orders in perturbation theory. The Fermi problem is usu- 
ally regarded just as a gedanken experiment, and remains untested, essentially 
because interactions between real atoms cannot be switched on and off. 

In this section we give a complete description of the problem in a physical 
framework in which predictions can be verified. This framework will be circuit 
QED which can be regarded as a 1-D version of Quantum Electrodynamics 
(QED) with two-level (artificial) atoms, a testbed which makes it possible to 
control the interaction and tune the physical parameters. We complete previous 
descriptions made of the problem and explain how there are no real causality 
issues for Fermi's two-atom system. We give an explicit non-perturbative proof 
of strict causality in these setups, showing that the probability of excitation 
of qubit B is completely independent of qubit A for times t < r/v and for 
arbitrary initial states. As a matter of fact, this comes as a manifestation of 
the nonsignaling character of the quantum theory [14]. We also show how this 
is compatible with the existence of nonlocal correlations at times < t < r/v, a 
fact pointed out in various theoretical proposals to entangle qubits at arbitrarily 
short times [HI HO] (and chapters 3 and 4 of this Thesis). More precisely, we give 
a non-perturbative proof of the fact that the probability of B being excited and 
A in the ground state is finite and r-dependent at any time, even for t < r/v. 
We provide also a physical and intuitive explanation of why the conclusions in 
|16| , even if mathematically sound, do not apply to the causality problem. At 
the end of the Section we discuss the time dependence predicted in our model 
for the various excitation probabilities and suggest a feasible experimental test 
of causality using superconducting circuits. 

5.1.2 There are no causality problem in Fermi's two-qubit 
system 

In what follows we focus on a practical setup of circuit-QED, with two qubits, 
A and B, interacting via a quantum field. The qubits have two stationary states 
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|e) and \g) separated by an energy Ml and interact with a one-dimensional field, 
V(x), ( |4.16 1 which propagates along a one-dimensional microwave guide that 



connects them. This field has a continuum of Fock operators [a^, a\,\ — 5{k—k'), 
and a linear spectrum, cj). = w|fc|, where v is the propagation velocity of the 



field. The normalization and the speed of photons, v (4.19) depend on the 
microscopic details such as the capacitance and inductance per unit length, c and 
We will assume qubits that are much smaller than the relevant wavelengths, 
A = v/fl, and are well separated. Under these conditions the Hamiltonian, 



point-like interaction between them (4.21) 



H = H Q + H[, splits into a free part (4.20) for the qubits and the field and a 



The original formulation of the Fermi problem begins with an initial state 

\in) = \e A g B 0) (5.1) 

in which only qubit A has been excited, while B and the field remain in their 
ground and vacuum states, respectively. The total probability of excitation of 
qubit J is the expectation value of the projector onto the excited state 

n = \ej)(ej\. (5.2) 

In the Heisenberg picture 

p eJ = (in\ V e j(t) \in) , J e {A, B}. (5.3) 

We will prove that for vt < r the probability P e s is completely independent of 
the state of qubit A for all initial states. In the Heisenberg picture this amounts 
to showing that there appears no observable of A in the projector V B {t) for 
vt < r. Our proof begins by solving formally the Heisenberg equations for Vj 

r}(t) - 7>}(0) = -jrf* dt'a^t'Mxj, t'). (5.4) 

Integrating the Heisenberg equations of the operators and a\ and insert- 



ing them in Eq. ( |3.3[ ), the total field evaluated at x in Heisenberg picture is 
decomposed 

V(x,t) =V (x,t) + V A (x,t) +V B (x,t) (5.5) 
into the homogenous part of the field 

/>oo 

V (x,t)=i dk^/N^e^-^ak + R.c. (5.6) 



and the back-action of A and B onto the field 



V J (x,t) = ^ j ^x (5.7) 
x / a x j{t') / oj k e lk(x - Xj) - luJk(t - t " ) dkdt' + K.c. 
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Eqs. (5.4) translates into a similar decomposition for the probabilty V B with 



three terms 

n(t)=V e m (t)+V e BB (t)+V c BA (t) (5-8) 
which are proportional to Vq,V b and Va, respectively. The only explicit de- 



pendence on A may come from V BA through V A (x B ,t). Manipulating Eq. (5.7) 
gives 



V A {x B ,t) 



-2nd A N d 
h dr 



t- 



(5.9) 



where the Heaviside function 9 shows that strictly V BA (x B ,t) = for vt < r, 
and no such dependence is possible. We still have to analyze a possible implicit 
dependence on A through V BB , whose expression is 



V BB (t) 



idlN rt 



r, 2 



dt' / dt"a v B {t')a B {t" 

n JO 



dk uike' 



-iwk(t'—t") 



+ H.c. 



(5.10) 



The only implicit dependence could come through the evolution of cr B y (t), but 
again this is not the case. Since 



0. 



(5.11) 



the evolution of u x B does not involve the field in any way, and for a v B (t) we have 
that 

& y B (t) = Qa B (t)/2 - ^V(x B ,t)a%(t) (5.12) 



so using again Eq. (5.5) and Eq. (5.9) we see that the A-dependent part of P BB 
is for vt < r. Thus V B may be finite but is completely independent of qubit 
A for vt < r, as we wanted to show. 

So far, we have demonstrated that although P B {t) is non-zero for vt < r, 
the only non-zero contribution is V Ba , which is not sensitive to the qubit A and 
thus cannot be used to transmit information between the qubits. Now we will 
show that this result is compatible with the existence of correlations for vt < r. 
For instance, we consider the probability of finding qubit B excited and qubit 
A on the ground state P eB ,gA, which is: 



P eB , gA = {in\V e B (t)V 9 A {t)\in) 



where 



V 9 



(5.13) 



(5.14) 



Using Eqs. (5.4), (5.5) and (5.8), we find a term in this probability which is 
proportional to V BB V AA and thus to V B (x B , t) Va{xa, t). From Eq. (5.7) we 
obtain: 

Vj{x J ,t) = 2 ^-j t {a%t)e{t)} (5.15) 
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Therefore, we conclude that in (5.131 there is an unavoidable dependence on 
A at any t > 0, but this is not a causality violation because correlations alone 
cannot transmit information. 

At this point it remains a single question: How can the A-dependent part of 
P e B be zero while the one of P e B, g A is nonzero for vt < r? To better understand 
it we need less formal results that rely on perturbative expansions, but we 
would like to remark here that the conclusions above are valid to all orders in 
perturbation theory. To obtain the total probability of excitation of qubit B 
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Figure 5.1: The terms of order dj, d 2 j and dj contributing to the amplitude for 
exciting qubit B. 

P e B to a given order in perturbation theory, one has to expand to a certain order 



the operators appearing in Eqs. (5.4 1, (5.5), (5.7). The different terms in the 



expansion can be related to the probabilities of the different physical processes 
involved. Fig. |5.1| shows the diagrams of the different amplitudes contributing 
to P b b up to the fourth order in dj. The lowest order amplitude contributing 
to a hnal excited B qubit is of order dj, which means that terms up to order 
d 3 j have to be considered. The only terms leading to this final state will be 

M W =V B , M (2) = X + UaVb, M {3) = A'Vb + UaVaVb + VbU b Vb + 5M {3) 

(5.16) 

where Uj (Vj) represent the amplitude for single photon emission at qubit J 
when the qubit is initially in the ground (excited) state, x is the amplitude for 
photon exchange, Aj are the radiative corrections of qubit J, and finally SA4^ 
is the amplitude for photon exchange accompanied by a single photon emission 
at qubit A. Notice that some of these processes are only possible beyond the 
rotating wave approximation, which breaks down for strongly coupled circuit- 
QED setups |112j as the ones considered later. Keeping only terms up to fourth 
order, we have for the probability to get B excited at a time t 

P eB {t) = \M (1) \ 2 + \M (2) \ 2 + 2Re{M (1) "M {2) } 

+ 2Re{M iiy M {3) } + 0(d 5 ) (5.17) 
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The final states in M.^ are orthogonal to those in M. 1 ^ and to the three 
photon terms in Ai^ . Hence, their interference vanishes. Besides, we are only 
interested in the A-dependent part of the probability, so we can remove the 



r-independent terms left in (5.171, marking the remaining contributions with a 
superscript 

pW(t) = \M {2) f T) + 2Re{M (1) * 6M {3) } + 0(d 5 ). (5.18) 

The first term actually gives P e s,gA up to the fourth order, it is positive and A- 
dependent at all times, as shown in Fig. |5.2fr . The second term is not a projector 
onto any physical state, but an interference term which has the effect of canceling 



out exactly the first term for vt < r but not for vt > r (cf. Fig. 5.2 ;>). In a 
nutshell, interference seems to be the physical mechanism that operates at all 
orders in perturbation theory to give the causal behavior of the total probability 
of excitation that we had previously shown. (See Appendix C for more details 
on the computation of 5A4^ 3 \) 

These perturbative results cast new light on the controversy on the Fermi 
problem and help us understand why our results do not contradict those of 
Hegerfeldt [16]. Hegerfeldt proved mathematically that the expectation value 
of an operator consisting of a sum of projectors cannot be zero for all the 
times vt < r, unless it is zero at any time. Indeed, the expectation value of 
Vg (t) cannot be zero for all vt < r, for it always contains the contribution Vg 



from Eq. (5.6). However, as we showed non-perturbatively, the actual relevant 
question for causality is whether the expectation value of V% A {t) vanishes for 
vt < r or not, since only this part of the probability is sensitive to qubit A 
and could be used to transmit information. Besides, according to our above 
perturbative results to fourth order, the r-dependent part of the probability, 
that is the expectation value of V e BA [t\ is not a mere sum of projectors, but 
also contains interfering terms. Thus, Hegerfeldt 's result does not apply and 
V% A (t) can be zero for vt < r as is actually the case. Both results are in accord 
with a general fact of Relativistic Quantum Field Theory: two global states 
can not be distinguished locally with the aid of a local projector annihilating 
one of the states, since the local observable algebras are Type III von Neumann 
algebras (See [171 1116] for a discussion). 



5.1.3 Experimental proposal 

We will now suggest an experiment to check the causal behavior of P e B- For this 
we need to control the interaction time at will to access the regions at both sides 
of t = r/v. This, which is highly unrealistic with real atoms, becomes feasible 
in circuit-QED. While the ideas are valid for both inductive and capacitivc 
couplings, we will focus on using a pair of three-junction flux qubits [1 171 1118"] . 
Each of the qubits is governed by the Hamiltonian 

Hoj = \eja*j + ±Ajo*j. (5.19) 
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Figure 5.2: (a) Pe B \gA an< ^ O 3 ) ^eS versus v t/ r f° r &r/v = § (blue, crosses), w 
(red, squares), and 27r (green, circles) with Ka.b — 0.0225. For vt < r the qubits 
are spacelike separated, but there are correlations between them and figure (b) 
shows the expected causal behavior, (c) P^ (blue, circles) and lA^ 1 ^ 2 (green, 
squares) vs. vt/r for = 0.20, Kb = 0.04 and a separation of one wavelength 
r = 2ttv/Ha.b- With this data and f2/27r ~ 1GHz we have At ~ Ins. (Color 
online) 



The energy ej = 2I p S^ x j, is approximately linear in the external magnetic flux, 
S$ x j, measured from the degeneracy point, and we assume that the gap Aj is 
fixed. The result is a qubit energy difference 

^j(^j) = ^(2I p 5<$> xJ y + A 2 j. (5.20) 
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The coupling between the qubit and the microwave photons is ruled by the 
dimensionless ratio 

*J= W =2 fo/ n ') 2 ' (5-21) 



Here 



9.1 = - JX f — (5-22) 



is the coupling strength between a qubit and the cavity that would be obtained 
by cutting the transmission line to be perfectly resonant with the qubit tran- 
sition. These numbers enter the qubit excitation probability computed before 



(5.181 through the product 



P^{t)ocK A K B . (5.23) 

Since 

Kj oc (5.24) 

we may use the external fluxes to move from a weakly coupled regime with 
no qubit excitations, £!j <C g, to the maximum coupling strength, £lj ~ 
A j (6$ xJ = 0). 



Let us first discuss how to prepare the initial state (5.1) of the Fermi problem. 
We assume that the system starts in a ground state of the form \g A Qb 0)- This 
is achieved cooling with a large negative value of 5<& x j on both qubits, which 
ensures a small value of g/ilj and Kj. We estimate that couplings g/ftj < 0.15 
and ttj ~ 1.5GHz lower the probability of finding photons in the initial state 
below 5 x 10~ 3 , both for vacuum and thermal excitations. Both magnetic fluxes 
are then raised up linearly in time, 

6<P xJ = ajt, (5.25) 

to prepare the qubits. Using a Landau-Zener analysis [119] of the process we 
conclude that an adiabatic ramp 

a B < nAl/Aklp (5.26) 

of qubit B followed by a diabatic ramp [1201 1121) 

a A > A 2 A /h2I p (5.27) 

of qubit A, leads to the desired state \e A gB 0) with a fidelity that can be close 
to 1, depending only of a A , a B as derived from the Landau-Zener formula: 

Pj = e~^ T1 ^. (5.28) 

Note that the minimum gap has to be large enough to ensure that the 
qubit-line coupling of B remains weak and the qubit does not "dress" the field 
with photons. 
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Once we have the initial state, both magnetic fluxes must take a constant 
value during the desired interaction time. After that, measurements of the 
probability of excitation of qubit B can be performed with a pulsed DC-SQUID 
scheme [12211102] . The timescale of the "jump" of the probability around t = r/v 
for qubit frequencies in the range of GHz and a separation of one wavelength 
r = 2irv/fl is 

At ~ Ins (5.29) 



[Fig. 5.2:]. Although the total measurement of the SQUID may take a few /is, 
the crucial part is the activation pulse (~ 15ns) in which the SQUID approaches 
its critical current and may switch depending on the qubit state. During this 
activation period the SQUID and the qubit are very strongly coupled (g ~ 
GHz), j!22j and the qubit is effectively frozen. The time resolution of the 
measurement is thus determined by the ramp time of the activation pulse, which 
may be below nanoseconds. Among the sources of noise that are expected, 
the short duration of the experiment, well below T\ and T 2 of usual qubits, 
makes the ambient noise and decoherence pretty much irrelevant. Thermal 
excitations of the qubits and the line may be strongly suppressed by using 
larger frequencies (> 1.5GHz). The most challenging aspect is the low accuracy 
of SQUID measurements, which are stochastic, have moderate visibilities |122j 
and will demand a large and careful statistics. 

On the technical side, it is important to choose carefully the coupling regimes. 
If we wish to compare with perturbation theory, we need Kj <C 1. However, 
at the same time the product KaKb must take sizable values for P e s oc 
KaKb to be large. And we need to discriminate the causal signal from the 
r-independent background of the probability of excitation, whose main con- 
tribution is IvW^I 2 oc Kb- Thus, a good strategy would be to work with 
Ka > Kb- In Fig. |5.2| ; we show that it is possible to achieve a regime in which 
the perturbative approximations are still valid and the r-dependent part of P e B 
is comparable to jTW^ 1 ^ 2 in the spacetime region of interest vt ~ r. 



5.1.4 Conclusions 

In this section, we have considered a system of two superconducting qubits 
coupled to a transmission line, which can be suitably described in the framework 
of 1-D QED with two-level (artificial) atoms. Starting from an initial state with 
qubit A excited, qubit B in the ground state and no photons, we have illustrated 
the causal character of the model showing that the probability of excitation of 
qubit B is completely independent of qubit A when vt < r. We have also shown 
that this is in agreement with the existence of nonlocal correlations and we have 
used perturbative computations to see the physical mechanism underlying the 
causal behavior. Finally, we have suggested an experiment feasible with current 
technology that would solve the controversy on the Fermi problem. 
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5.2 Quantum fluctuations and short-time quan- 
tum detection 



5.2.1 Introduction 

A counterintuitive direct consequence of the breakdown of the RWA is that 
a detector in its ground state interacting with the vacuum of the field has a 
certain probability of getting excited and emitting a photon. There is however 
not a widespread consensus on the physical reality of this effect. Introducing 
counterrotating terms is interpreted by some to be a problem as the processes 
described by those terms seem virtual. It seems difficult to accept that a detector 
in a ground state in the vacuum could get excited. As a matter of fact, there 
have been attempts of suggesting effective detector models by imposing this 
phenomenon to be impossible [123] , We should however recall here that these 
peculiar effects should not be that discomforting. They are linked to the fact 
that the initial state considered has not definite energy, since the state "detector 
and field in their ground states" is not an eigenstate of the full Hamiltonian 
beyond RWA. Indeed, in section 4.1 we have devised an experimental protocol 
to detect the trace of this ground-state qubit self-excitations in a simplified-setup 
with a qubit in a cavity. 

In this section we will keep using the theoretical framework applicable to 
circuit QED, without imposing any additional constraint. We will study the 
following setup: a source A initially excited, a detector B initially in the ground 
state and both interacting with the electromagnetic field in its vacuum state. 
If the detector clicks at a given time, does it mean that the source is now in 
the ground state? This problem amounts to compute the probability of decay 
of the source, conditioned to the excitation of the detector. We will show that, 
unlike Glauber's RWA detector in which this conditioned probability would be 
equal to 1 at any time, this circuit QED detector only achieves this value at 
long times due to the impact of non-RWA effects, like the ground-state qubit 
self-excitations described above. We will see how these theoretical results have 
to be taken into account for the interpretation of the readouts of real ongoing 
experiments. 



5.2.2 What does a detector's click mean? 

More precisely, let us consider an initial moment t — where A is excited, B is 
in its ground state and there are no excitations in the transmission line, that is, 



the initial state would be Eq. 5.1 After a certain time t, if we measure qubit 
B and it results excited, that would naively lead us to think A has decayed 
and produced a photon which has then later been absorbed by B. We intend 
to proof otherwise by quantifying what information about the state of A can 
be extracted by knowing qubit B state after a certain time t. For that we will 
compute the probability P g A/eB{t) of A to have decayed at a certain instant t, 
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conditioned we have measured B excited at that same moment: 



P. 



gA/eB 



(*) 



P, 



eB,gA 



P, 



B 



(5.30) 



PeB,gA being the probability of A being in the ground state and B excited and 
P e B the total probability of excitation of B, as they have been defined in the 
previous section 5.1. Within RWA P g A/eB = 1 at any time. But beyond RWA 
and up to the fourth order in perturbation theory, we have seen in the previous 
section that P e s is given by Eq. 5.17 and that 



Pp 



,gA 



\M 



(2)|2 



(5.31) 




Figure 5.3: P g A/eB(t) (5.30) in front of Kt for three different values of the 
coupling strength of K = Ka = Kb = 7.5 • 10~ 3 (solid, blue, circles), 
1.5 10 _2 (dashed, green, crosses), 7.5 • 10~ 2 (dashed, black, squares), 
three cases 2tt^ = 1 and f2/ (2tt) ~ 1 GHz (fl = 0,^ = ^b)- 



In the 



The effect of the non-RWA contributions to the evolution of PgA/eB^) can 



be seen in Figs 5.3 and |5.4| where the consequences of changing the coupling 
and the distance between qubits are considered. The first thing we notice in 
Fig|5.3| is that for short times the information provided by the detector is not 
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Figure 5.4: VgA/eB^) (5.30) in front of t in s for three different values of the 
distance 2nj — 0.5 (dotted, blue, circles), 0.75 (dashed, green, crosses), 1 (solid, 
black, squares). In the three cases the coupling strength K = Ka = Kb = 
1.5 • 10~ 2 and n/(2n) = 1 GHz (fi = il A = Cl B ). 



very much related to the state of the source, that is, self-excitations and other 
non-RWA phenomena dominate over the photon exchange between source and 
detector. For the cases considered, only at interaction times t > 1 ns ~ 1/fl the 
conditioned probability converges to the RWA prediction, that is, the excitation 
of the detector is a reliable way to detect the decay of the source. Since the non- 
RWA contributions are more relevant for large couplings and short distances, 
the convergence is faster as the distance grows and the couplings diminish, as 
can be seen in Figs |5.3| and |5.4| It is convenient to remark here that the ripple 
frequency we see for instance in Fig. |5.3| comes from higher harmonics of the 
qubit frequency 51 (= Qa — &b in our case) and so does not represent the 
emergence of a new time-scale. It can be thought as a process similar to that 
of a Rabi oscillation, where the qubits would be absorbing in cycles (in a self- 
deexcitation fashion) the photons previously emitted in self-excitations. 

The above theoretical results could have an impact in real experiments of 
circuit QED. In particular, a typical setup to measure the internal state of a flux 
qubit coupled to a transmission line consists of a SQUID surrounding the qubit. 
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Although the total measurement process could take up to tens of nanoseconds, 
most of the time the coupling SQUID-qubit is much stronger than K, [122] and 
the dynamics qubit-transmission line is effectively frozen. Thus this dynamics 
is only important during the activation of the SQUID, a process that may be 
in the nanosecond regime. For those measurement times, as we have proved, 
self-excitation effects cannot be disregarded and should manifest themselves. 

Besides, it should, in principle, be possible to prepare experiments in the near 
future to test our predictions directly. We here intend to give just a rough sketch. 
Such experiments would involve the preparation of the system at t = in the 



initial state of Eq. 5.1 the switching of the interaction for a certain time t (in the 
line of previous proposals, as in section 5.1) and then the SQUID-measurement 
of both qubits A and B. By repeating the experiment several times, the result 
frequencies should theoretically match our probability predictions. 
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Chapter 6 



Quantum simulations of 
relativistic dynamics. 
Majorana physics 

"But is this really the poet" I asked. "There are two brothers, I 
know; and both have attained reputation in letters. The Minister I 
believe has written learnedly on the Differential Calculus. He is a 
mathematician, and no poet." 'You are mistaken; I know him well; 
he is both. As poet and mathematician, he would reason well; as 
mere mathematician, he could not have reasoned at all, and thus 
would have been at the mercy of the Prefect." (Edgar Allan Poe, 
The purloined letter) 

"II mare mi ha rifiutato e ritornero domani all'albergo Bologna, viag- 
giando forse con questo stesso foglio." (Ettore Majorana, last letter 
to Prof. Carelli. English free translation: "The sea rejected me and 
I'll be back tomorrow at Hotel Bologna traveling perhaps with this 
sheet.") 

In chapters 4 and 5 circuit QED was considered as a quantum simulation of 
matter-radiation interaction. In this chapter we will deal with different quan- 
tum simulations, in particular quantum simulations of Relativistic Quantum 
Mechanical systems. To this end, we will introduce in section 6.1 a new kind 
of pseudo-Hamiltonian, referred to as "Majorana Hamiltonian" . An example of 
this class is the object appearing in the Majorana equation, that is, the relativis- 
tic equation of a fermion with a Majorana mass term instead of the standard 
mass term of the Dirac equation. We will propose a trapped ion simulation of the 
1-D version of this equation in section 6.2. Finally, in section 6.3 we will show a 
way to extend the quantum simulations of free relativistic quantum-mechanical 
systems to systems under the action of potentials, both for single-particle and 
bipartite systems, by simulating the potentials with free hamiltonians. 
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Note: Throughout this chapter we will use natural units (H = c = 1) 

6.1 On Majorana Hamiltonians 

6.1.1 Introduction 

The existence of spin 1/2 fermions that are their own antiparticles, a possibil- 
ity opened up by Majorana [231 , is under close scrutiny in terms of neutrino 
properties, theoretical schemes beyond the Standard Model or even in solid 
state systems [124] . These particles could have a mass with properties crucially 
different from those of the standard Dirac mass of charged fermions. Their 
very equation -called Majorana Equation- was introduced by Jehle [23] and by 
Case [125] , who also made a thorough study of their behaviour under space-time 
transformations, quantized the field and also analyzed their possible (weak) in- 
teractions. Primers on these topics from a modern perspective can be found in 
Refs. [UnHMj and [Eg]. 

On studying the quantum mechanics of a fermion obeying the Majorana 
equation we have identified a structure of much more general applicability. In 
this section, we present this structure as a very general theoretical framework. 
We introduce a new kind of generalized Hamiltonians with different mathemati- 
cal features of the textbook Hamiltonians: they are neither linear nor antilinear, 
and thus the definition of hermiticity is questionable. We will refer to them as 
Majorana Hamiltonians, since the Majorana equation is a good example, al- 
though their scope is more general. Our main result is that these objects enjoy 
the physically meaningful property of inducing a temporal evolution which con- 
serves the norm, as expected from total probability conservation. However, this 
result comes with a bag of surprises: particles with a Majorana Hamiltonian do 
not have stationary states, probability amplitudes are thus not conserved and 
an initial global phase is observable in the evolution of the expectation values of 
physical magnitudes. Moreover, even the notion of density matrix is no longer 
adequate for a description of this new Majorana dynamics. 

Alternatively, we will also show that all the above physics can be equivalently 
analyzed with a standard Hamiltonian in a real Hilbert space of higher dimen- 
sionality, which is actually the way in which a Majorana Hamiltonian would 
be implemented in the lab, as we will show in section 6.2. As we shall see, 
this complementary focus sheds light on the shocking features of the Majorana 
physics commented above. 

6.1.2 Majorana Hamiltonians. 

In what follows, we will present the general formalism, that will be illustrated 
in the text with some cases of the Majorana equation. We will refer to the 
following object as a Majorana Hamiltonian: 

M v = A + irjBK , (6.1) 



where 77 is a complex number such that 



M = 1, (6-2) 

A and B are linear Hcrmitian operators and the novelty is in the complex 
conjugation operator K, which obeys the following properties: 

K[oal>] = a*i/j* ,K 2 = 1, (6.3) 

and, for all ij) and <fi in the Hilbert space, 

{<t>,kij) = ^,k ( j>), (6.4) 

The general antilinear product operator BK results antiunitary if B is uni- 
tary [129J. We restrict this last operator by the condition 

{K,B} = 0, (6.5) 

i.e. K and B anticommute. We can express this by saying that B is an imaginary 
operator (with respect to k). Thus, is neither linear nor antilinear, but it 
is well defined on (a domain of) the Hilbert space. 

A very important result is that if the evolution in a Hilbert space is governed 
by a Majorana operator, that is by the generalized Schrodinger equation 

idttj) = M v i> = (A + ir]Bk)ip, (6.6) 

then the norm is still conserved, 

d t (m,m) = o- (6.7) 

and we can still think of \ip\ 2 as a probability distribution. The proof is simple 
but subtle. First of all, the term with A leads to norm conservation due to the 
standard hermiticity argument. Next, one can readily check that 

(ip,BKip) = -(i/>,Bkip) = 0. (6.8) 

This last result follows from the Hermiticity of B, the fact that k and B anti- 
commute, and the relation ((p, ktp) = (ip,k^(f>). It is important to remark that 
while the norm is preserved, amplitudes are not. 

The formalism introduced so far appears to respect all the postulates of 
quantum mechanics, but there is one usual convention that will have to be 
dropped: the equivalence of state vectors under global phases. This follows 



from the fact that the constant 77 in Eq. (6.1 1 can be traded for a global phase 



change. In other words, if ip(t) is a solution of the generalized Schrodinger 



equation (6.6) associated to the phase 77, then 

i> 1 (t)=V~ 1/2 i>(t) (6-9) 
will evolve according to a different operator, 

Mi = A + iBK, (6.10) 
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and it will behave, in general, differently. 

This inequivalent evolution of states initially related through a global phase 
entails that the standard description of mixed states with projector density 
matrices is not suitable. A mixed state in the present context will be a collection 
of couples whose elements are unitary vectors ipi and corresponding probabilities 
Pi, with evolution given by 

p M (t) = { Pl ,Mt)}li (6.11) 



with each ipi (t) a solution of Eq. ( 6.6 ) . This gives rise to the following interesting 



point: consider ^i(O) and </>(0) such that 

(p(0) = exp{itp)ip(0). (6.12) 

Then there exists a family of mixed states {(p, ip(0)) , (1 — p, 0(0))} which result 
in identical measurements at time but are nonetheless generically inequivalent 
under evolution. In other words, density matrices can be the adequate tool for 
the description or computation of static quantities; but not for dynamic ones. 
We will illustrate all the above in the Examples section 6.1.4. 



6.1.3 Hamiltonization. 

We shall now provide an alternative description of the previous results in terms 
of a doubled (real) Hilbert space. We introduce the objects 

P+ = i(l + K), P- = \(l-K), (6.13) 

which are not projectors due to the antilinearity of K though satisfy the useful 
relations 

P| = P±. (6.14) 
Using these operators we reconstruct the isomorphism 

between the complex Hilbert space H. and the direct sum HS 2 ^ of two real Hilbert 
spaces, P+H, —i P—H. The norm is conserved under the isomorphism, naturally, 
but it is important to notice that the inner products are not preserved, namely, 

(*,*)= Re [(4>,iP)}- (6-16) 

An alternative way of writing the isomorphism is as [130 

* = ReV><8|O)+Im<0(g)|l), (6.17) 

and demanding that the coefficients be real. Note that the reverse mapping, 
from the doubled real Hilbert space to the initial one is achieved explicitly as 

iP = V^={l n il H ) * = Rc?/> + i\mip. (6.18) 
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In order to compute how the previous isomorphism acts on operators, we have 
to introduce the real and imaginary parts of an arbitrary operator, 

= RcO + iImO { Re0 = U A A¥ M J (619) 

Both the real and imaginary parts of a linear operator are also linear, but while 
the real part is Hermitian, Im O is anti-Hermitian. Using these components the 
Hermitian operator O acting on "H is mapped onto an operator 

6(2) = f Re - Im = Re 6 ® 1 2 - i Im 6 ® (7 2 , (6.20) 
\lmC ReC / 

acting on %( 2 ) that has all components real with respect to the induced complex 
conjugation. 

Using these tools we can find the evolution equa tion for the mapped vector 



€ "%( 2 ), up from the Majorana Hamiltonian (6.6). Taking for simplicity the 



operator without phase, M%, we obtain the Schrodinger-like equation 

.„ T ( ilmA B + iReA\ T - (2 \ T ,„.., 
id t ty = - ; j * = Af (2) *. (6.21) 

VB-iReA ilmA J V ; 

This equation deserves a detailed explanation. First of all, the operator Af( 2 ) 
is designed to act on the doubled complex Hilbert space H^BT-l, which contains 
our real Hilbert space, H^ 2 \ as a subset. Second, note the operator »M' 2 ' is real 
and thus if ^(0) is initially real, then it remains so throughout the evolution. In 
other words, if ^(0) G Ti^ C % ® %, it remains in the real sector at all times, 
preserving our isomorphism. Finally, while the doubled Hamiltonian Af( 2 ) is 
Hermitian and its eigenvalues are real, this only happens when this operator 
is defined on the complex doubled Hilbert space % ® %. As a consequence of 
this and the fact that ikf( 2 ) is purely imaginary, the eigenvectors corresponding 
to nonzero eigenvalues are necessarily nonreal, so they live outside . It 
follows that M% does not have any eigenvectors with nonzero eigenvalues, and 
the evolution under its generalized Schrodinger equation will of necessity be 
either constant or nonstationary. In other words, the Majorana Hamiltonians 



and Eq. (6.6) have no stationary states. 

Now we can understand the role of global phases in Majorana physics in a 
different way. Two states in the original Hilbert space H related by a global 
phase are mapped to states which are no longer related by a global phase in 
7^(2), but by a rotation 

0(0) = exp (i<p)i/i(0) 

t = t (6-22) 

$(0) = exp (-io- 2 <p)*(0). 

Thus clearly $(£) and vE'(i) evolve in a non-equivalent way, as do <fi(t) — VQ(t) 
and 4>{t) = 
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6.1.4 Examples 

We now apply our formalism to a series of cases whose first examples were given 
in the 1+1 Majorana equation [2~4l I125j : 

i dt tp = apip ± i m a 2 Kip, (6.23) 

where a is one of the Pauli matrices <Ji (i = 1, 2, 3), depending on the chosen 
representation of the Clifford algebra (see Appendix D). This equation is derived 
by replacing the ip in the mass term of the usual Dirac equation by its charge 
conjugate ip , as we will see with more detail in the next section and in the 
Appendix D. For simplicity and without loss of generality, we will use now the 
example with p = and m = 1: 

id t ip = ia 2 ip*. (6.24) 



Let us see how the previous statements are realized in (6.24 ). One can readily 
compute the solution, which is given by 

ip(t) = cos t ip(0) + sin t a 2 ip(0)* . (6.25) 

Norm conservation is an immediate consequence of the antisymmetry of a 2 . On 
the other hand, amplitudes are not conserved: consider 

#))=(T) ,<K0)=(T). (6.26) 



It is then the case that 



((f)(t),t(j(t)) =£sin2i. (6.27) 



The family of transformations given by (6.25) is not a family of linear or anti 



linear isometries (131) : we could define an evolution operator 

U(t) = cost 1 + sin* a 2 K, (6.28) 

which is not unitary. 



Let us, as before, study the case (6.24). Consider two initial vectors, 



V'(O) = (Vj , 0(0) = cxpMV(O). (6.29) 

We then obtain 

Furthermore, this phase is observable through the evolution in time of a measur- 
able magnitude. Consider for instance the expectation value of the Hermitian 
operator a 2 as a function of time: 

(a 2 )(t) = sin(2t) cos(2<p). (6.31) 
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Figure 6.1: Von Neumann entropy of the state p(t) Eq. (6.321 for ip — ir/2 in 
front of p and t. In the standard quantum mechanical picture, this would be a 
pure state for any value of p and t and so entropy would be zero. 



We now consider the mixed state defined by ip(t) with probability p and 4>(t) 
with probability 1 — p. If g(t) is the density matrix that would correspond to 
ip{t), then the density matrix p(t), corresponding to the mixed state at time t, 
is 

p(t) = (>(t) + (1 — p) sin tp sin(2t) (cos (pai — sin ya^) ■ (6.32) 



The von Neumann entropy of this mixed state is depicted in Fig. 6.1 as a func- 
tion of time and p for the case of <p = tt/2. In the standard quantum mechanical 
picture, the initial state would be considered as a pure state, and unitary evo- 
lution would conserve entropy. Indeed, if the state were the reduced state of a 
composite system this would mean that the global state is non-entangled. All 
the above suggests that the standard quantum information theory, as it is, does 
not fit to Majorana physics. 

As a general example of the hamiltonization we can consider Eq. ( |6.1[ ) with 
V = l, ' 

B = ma 2 , A — p x a 1 + p y a 2) (6.33) 
which accommodates both the 1+1 and 2+1 Majorana equation. Thus, using 



Eq. (6.21), we obtain the following Hamiltonian in the doubled real Hilbert 
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space: 

jfr<2) = ( Px<Ji i Py <J 2 + ma 2 \ 

\-iPyV2 + m<7 2 PxOl /' 

which is the Hamiltonian of a Dirac equation with 

f3 = ai<g> cr 2 , a 1 = l 2x2 ® ci, a 2 = ~a 2 ® a 2 . (6.35) 

It is straightforward to check that 

7° = P, 7 1 = Pa 1 = -i (Ti <g> (7 3 , 7 2 = /3a 2 = -t <r 3 ® l 2 x 2 (6.36) 

form a Clifford algebra, and also 7 , 7 1 . Thus Majorana equation in 2+1 and 
1+1 dimensions is mapped to a 3+1 Dirac equation with p z — and p a = 
p z = respectively. This result is more general: with this techniques the 3+1 
Majorana equation is mapped to a 7+1 Dirac equation with four components 
of the momentum set to 0. 

As another example, let us consider the Hilbert space of qubits, where a 
general operator A = a^a^ with — (l,(?) and real a M will be decomposed as 
Re A — a,{)l + a\G\ + 03(73, ImA = —ia 2 a 2 . Without losing generality, we can 
choose B = (7 2 . The scale of B would be fixed by the timescale. Thus we can 
explicitly write 

M^=( ° 2CT2 , ^ + *RelA (6 37) 

\a 2 — i~ReA a 2 a 2 J 

whose eigenvalues appear in pairs, 



Ai = |a| + \/l + al, A 2 = -Ai, 



A 3 = |a| - ^l + a 2 ,A 4 = -A 3 . (6.38) 
6.1.5 Conclusions and outlook. 

Putting all together, we have introduced a new kind of generalized Hamiltoni- 



ans Eq. (6.1). We refer to them as Majorana Hamiltonians since the Majorana 
equation is a good example, but the framework is far more general. These ob- 
jects are neither linear nor antilinear and thus there is not a clear definition of 
hermiticity, but they induce a physical dynamics in which the norm is conserved. 
This new physics has very interesting features: for instance, initial global phases 
are relevant since they can be observed in the dynamics of expectation values 
of physical magnitudes. As a consequence, the notion of density matrix is not 
adequate for a dynamical description of the system, and has to be replaced. Al- 
ternatively, all this can be seen within the usual Quantum Mechanics framework, 
since a Majorana Hamiltonian can always be mapped to a standard Hamiltonian 
in a real Hilbert space of double dimensionality. A global phase transformation 
in the original Hilbert space is not mapped to a global phase transformation in 
the new one, which explains the differences in the evolution. This mapping be- 
tween Hilbert spaces is actually the tool to implement Majorana Hamiltonians 
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and also antiunitary operations like time reversal, charge conjugation or partial 
transpose in the lab, as we will see in the next section. A deeper look into the 
physics of real Hilbert spaces will imply also a better understanding of the new 
Majorana physics, for instance the quantum information theory of Majorana 
systems, since it has been noted that quantum information properties in real 
Hilbert spaces are different from the standard ones |132j . [133] . 

6.2 Quantum simulation of Majorana equation 
and unphysical operations 

6.2.1 Introduction 

The Majorana equation [53], [125] is a relativistic wave equation for fermions 
where the mass term contains the charge conjugate of the complex spinor, ip c , 

— mtp c . (6.39) 

Here, $ = 7 M 9^ and 7^ are the Dirac matrices |20j, while the non-Hamiltonian 
character stems from the simultaneous presence of ip and ip c . The significance 
of the Majorana equation rests on the fact that it can be derived from first prin- 
ciples in a similar fashion as the Dirac equation |128j . Both wave equations are 
Lorentz invariant but the former preserves helicity and does not enjoy station- 
ary solutions. The Majorana equation is considered a possible model |126| for 
describing exotic particles in supersymmetric theories -photinos and gluinos-, 
or in grand unified theories, as is the case of neutrinos. Indeed, the discussion 
of whether neutrinos are Dirac or Majorana particles still remains open |124j . 
Relativistic quantum models can be simulated, but they can also emerge as a 
natural description of certain systems, as happens with the 2+1 Dirac equa- 
tion in graphene systems. Nevertheless, note that despite the similar naming, 
this work is neither related to the Majorana fermions (modes) in many-body 
systems, nor to the Majorana fermions (spinors) in the Dirac equation [1 26J . 

In order to simulate physics described by the Majorana equation, we have 
to solve a fundamental problem: the physical implementation of antilinear and 
antiunitary operations in a quantum simulator. In this section, we apply a 
mapping introduced in the previous section by which complex conjugation, an 
unphysical operation, becomes a unitary operation acting on an enlarged Hilbert 
space. The mapping works in arbitrary dimensions and can be immediately 
applied on advanced quantum simulation platforms. As a key application, we 
show how to simulate the Majorana equation in 1+1 dimensions and other 
unphysical operations -time reversal and complex conjugation- using only two 
trapped ions. This is completed with a recipe for measuring relevant observables 
and a roadmap towards more general scenarios, including the combination of 
Majorana and Dirac physics. Finally, we discuss further scopes of quantum 
simulations in the context of fundamental and relativistic quantum physics. 
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6.2.2 Quantum simulation of antiunitary operations 



Spinor 
mapping 



Complex 
spinor 



1 



Unitary operations 



V $ J 

Real bispinor 



P 

I 



Reverse Unphysical 
mapping oper ations 



Vr* = (<r 2 ig)(7 2 )* - 

Vic* = 8 « 




Figure 6.2: Diagram showing the different steps involved in the quantum simu- 
lation of unphysical operations in 1+1 dimensions. 



6.2.3 Simulating antiunitary operations 



There are three discrete symmetries |134j which are central to quantum me- 
chanics and our understanding of particles, fields and their interactions: parity, 
V, time reversal, T, and charge conjugation, C. None of these operations can be 
carried out in the real world: V involves a global change of the whole physical 
space, while C and T are antiunitaries. However, there is no apparent restriction 
for implementing them in a physical system that simulates quantum mechanics. 
We will focus on the study of antiunitary operations, which can be decomposed 
into a product of a unitary, He (for charge conjugation) or Uj- (for time rever- 
sal), and complex conjugation, Kip — ip*. We consider the mapping 6.15 of the 
quantum states of an n-dimcnsional complex Hilbert space, C„, onto an real 
Hilbert space, Ibn- This mapping can be physically implemented by means of 
an auxiliary two-level system, such that M> 2n G %2 (8> H n - In this manner, the 
complex conjugation of the simulated state becomes a local unitary Vjc acting 
solely on the ancillary space, 



jCtp = tp* -> Vic* = (<r z ® 1)*, 



(6.40) 



and thus physically implementable. Furthermore, unitaries and observables can 



be also mapped onto the real space, (6.201 preserving unitarity and Hermiticity. 
In addition to complex conjugation, unitaries and Hermitian operators, the 
proposed simulator also accomodates the antiunitary operations 



and 



C=U C K, 



T = U T JC. 



(6.41) 



(6.42) 



To this end, we have to choose a particular representation (see Appendix D) 
that fixes the unitaries Uc and U-j-, as will be shown below. 
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At this point, we possess the basic tools to simulate the Majorana equation 
in the enlarged space. The expression for the charge conjugate spinor is given 

by 

i> c = 77C 7 °/CV, (6.43) 
with C a unitary matrix satisfying 



c~ Vc = - (Yy 



(6.44) 



We illustrate now with the case of 1+1 dimensions (Appendix D). Here, a suit- 
able representation of charge conjugation is 



that is rjC = icr y , and the Majorana equation reads 



id t ip = a x p x tp - imayip*, 



(6.45) 



(6.46) 



where p x = 
Hamiltonian. 



—ihd x is the momentum operator. Note that Eq. (6.46) is not 
ihdtip = Hip, as is the case of Schrodinger and Dirac equations. 



This is due to the presence of a complex conjugate operation in the right-hand 
side of Eq. (6.46), which is not a linear Hermitian operator. Through our 



mapping (6.15) 



i>2 



(6.47) 



V $ 

the Majorana equation for a complex spinor becomes a 3+1 Dirac equation with 
dimensional reduction, p y ,p z = 0, and a four-component real bispinor 



/ 



ihd t ^ = [(1 ® <J x )Px - ma x ® o y \ 



(6.48) 



Note that, here, the dynamics preserves the reality of the bispinor VP and, in gen- 



eral, cannot be reduced to a single 1+1 Dirac particle. The result of Eq. (6.48) is 



even more general and the complex-to-real map in arbitrary dimensions trans- 
forms always a Majorana equation into a higher dimensional Dirac equation. 



Since Eq. (6.48) is a Hamiltonian wave equation, it can be simulated in a con- 



ventional quantum system while suitably encoding the Majorana dynamics. 

The mapping of wavefunctions into larger spinors may allow us not only to 
implement Majorana equations in the lab, but also to explore exotic symmetries 
and unphysical operations, otherwise impossible in nature. From Eqs. (6.40), 



(6.45), and (6.47), for the 1+1 dimensional case, we can deduce that charge 



conjugation is implemented in the enlarged space via the unitary operation Vc 



(6.49) 



We can do something similar with time reversal, defined as the change t — > (— t). 
In this case, we expect |128j (Appendix D) 



(6.50) 
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Figure 6.3: Scattering of a fermion against a linearly growing potential (inset), 
(a) Ordinary Klein process: a fraction of a Dirac fermion turns into an antipar- 
ticle, entering the potential, (b) At an instant of time we apply the time reversal 
operator T causing the particle to retrace its own trajectory, (c) Similar to (b) 
but now we apply charge conjugation, converting the particle in its antiparticle. 
(d) Scattering of a Majorana fermion, which propagates through the potential. 
Parameters are m = 0.5 and V{x) = x, in dimcnsionlcss units. 

where the time variable r = — t and the modified spinor 

^>'(t) = T0(t). (6.51) 

In order to preserve scalar products and distances, the time reversal operator 
must be an anti-unitary operator and thus decomposable as the product 

T = U T IC. (6.52) 

In 1 + 1 dimensions, imposing that the Hamiltonian be invariant under time 
reversal, 

H' = T~ X HT, (6.53) 
implies that the unitary satisfies 

Wy-(i(j x d x )Ur = -icr x d x , (6.54) 

with a possible choice being 

U T = <J Z . (6.55) 
In other words, in the enlarged simulation space 

Tip = U r Kip V T V = (<J Z ® a g )^l. (6.56) 

See Fig. |6.2| for a scheme of the simulated symmetries. 



Equation ( 6.48 ), a Dirac equation in 3+1 dimensions with dimensional reduc- 
tion p y ,p z — 0, can also be interpreted as a recipe for the quantum simulation of 
the Majorana equation in the laboratory. In a recent experiment, the dynamics 
of a free Dirac particle was simulated using a single trapped ion [22] , a quantum 



platform that has proved instrumental for quantum information implementa- 



tions. Unfortunately, Eq. (6.48) has a more complex structure and a different 
setup is required, which is outlined in the Appendix E. Moreover, the encoded 
Majorana dynamics requires a systematic decoding via a suitable reverse map- 
ping of observables, as we show below. In short, the real bispinor $ € I4 can be 
encoded in the internal state of two ions, while the position and momentum of 
the Majorana particle are mapped onto the quadratures of a collective motional 
mode, e. g. the center-of-mass mode, of the ions |135[ 122 ) . The Hamiltonian of 



Eq. (6.481 can be implemented term-by-term, in principle, in the trapped-ion 
system by a number of lasers coupling the motional and internal states of the 
ions. However, our proposal is valid for a general quantum simulator and we do 
not discard its implementation in other quantum platforms. 

A relevant feature of the Majorana equation in 3+1 dimensions is the con- 
servation of helicity. A reminiscent of the latter in 1 + 1 dimensions is the 
observable called hereafter as pseudo-helicity 

£ = °xPx- (6.57) 



This quantity is conserved in the 1+1 Majorana dynamics of Eq. (6.461 but 
not in the 1+1 Dirac equation. We will use this observable to illustrate mea- 
surements on the Majorana wavefunction. The mapping for operators can be 
simplified if we are only interested in expectation values. Reconstructing the 
complex spinor 

ip = (6.58) 

with matrix 

M = ( 1 il ) (6.59) 
associated with Eqs. ( |6.47 1 and ( 6.48[ ), we can write the following equivalence 



(C% - <-0|O|-0> = (*|MtOM|*) =: (6)*. (6.60) 

According to this, in order to measure the pseudo-helicity S, we have to measure 

£ = M^a x p x M = (1 <g> a x — a y <g> a x ) ® p x (6.61) 

in the enlarged simulation space. In an ion trap implementation, the first term 
of this observable, (1 ® <r x ) ® p x , is measurable with recently developed tech- 
niques i.22 . The second term is a three-operator correlation, {a y ® a x ) ® p x , 
and will require a specific design with measurements involving short-interaction 
times |136j . as explained in the Methods section. 

We want to emphasize that the previous mappings and the implementation 
of discrete symmetries are not only valid for Majorana equations, but also for 
Dirac spinors. Equally interesting is the possibility of combining both Dirac 
and Majorana mass terms in the same equation |126j . 

ifop = niMipc + rriDip, (6.62) 

which still requires only two ions for a 1+1 quantum simulation. It also becomes 
feasible to have CP violating phases in the Dirac mass term, mo exp(i#7 5 ). 
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Furthermore, we could study the dynamics of coupled Majorana neutrinos with 
a term Mip c , where M is now a matrix and ip — ip(xi, x^) is the combination of 
two such particles, simulated with three ions and two vibrational modes. 

So far, we have presented a complete toolbox of unphysical operations, C, T, 
and /C, that are available in the proposed quantum simulator. We can combine 
all these tools to study dynamical properties of the transformed wavefunctions. 
To exemplify the kind of experiments that become available, we have studied 
the scattering of wavepackets against a linearly growing potential, 

V{x) = ax. (6.63) 

It is known that repulsive potentials are partially penetrated by Dirac parti- 
cles [5D], an effect called the Klein paradox [137, 138, 1139] . This is shown in 



Fig. 6.3 1, where a Dirac particle splits into a fraction of a particle, that bounces 
back, and a large antiparticle component that penetrates the barrier. This nu- 
merical experiment has been combined with the discrete symmetries and the 
Majorana equation. In Fig. |6.3b we show a Dirac wavepacket that suffers the 
time reversal operation some time after entering the barrier: all momenta are 
reversed and the wavepacket is refocused, tracing back exactly its original tra- 



jectory. In Fig. 6.3 : we repeat the same procedure but using charge conjugation. 
This operation changes the sign of the charge turning a repulsive electric poten- 
tial into an attractive one, which can be easily penetrated by the antiparticle. 
In our last example, Fig. |6.3p l, we show the scattering of a Majorana particle. 
While the evolution is not so smooth — there are no plane wave solutions in the 
Majorana equation — , we can still identify a wavepacket penetrating the barrier, 
showing a counter-intuitive insensitivity to the presence of the barrier, that will 
be explained in the next section 6.3. 



6.2.4 Conclusions 

In summary, we have introduced a general method to implement the quan- 
tum simulation of unphysical operations and the non-Hamiltonian Majorana 
equation in a Hamiltonian system. To this end, we have designed a suitable 
mapping that enlarges the simulation space by means of an ancillary system to 
allow for complex conjugation, charge conjugation, and time reversal. We have 
exemplified the implementation of the 1+1 dimensional case in the context of 
trapped-ion physics. The flexibility of this protocol allows to explore a novel 
front of quantum simulations, that of unphysical operations and exotic quantum 
rclativistic processes that go beyond ordinary Schrodinger and Dirac quantum 
mechanics. 
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6.3 Quantum simulation of relativist ic poten- 
tials without potentials 

6.3.1 Introduction 

Last years had witnessed an increasing interest in simulating dynamics com- 
ing from the Relativistic Quantum Mechanics realm with physical systems of 
Quantum Optics, such as trapped ions. Striking theoretical predictions like Zit- 
terbewegung and Klein paradox |137j has been observed in these simulations. In 
particular, the proposal of simulation with one trapped ion [135] of the single 
particle free 1-D Dirac equation has been successfully implemented in the lab 
[22] . This is also the case for the single particle 1-D Dirac equation with some 
external potential [138] whose simulation involves two ions [139] . In the last 
section we have proposed that the free Majorana equation and unphysical oper- 
ations like complex conjugation, charge conjugation and time reversal can also 
be simulated with two trapped ions. Besides, two-body Dirac equations have 
been the subject of recent research [140] . 

In this section we will show that the same setups built up for the simulations 
of the free single particle Dirac and Majorana equations can also be used for 
simulations of these equations with the addition of some external potential, for 
a large class of potentials. This is based in the following idea, which is the main 
result of this paper: some states which are solutions of the Dirac or Majorana 
equation with the potential can be related through a phase transformation with 
a solution of the free corresponding equation. So, if we want to simulate the 
dynamics of a given state under certain potential the transformation tells us 
which is the state whose dynamics under the free equation does the job. In 
other words, the potential is codified in the phase involved in the transformation. 
The trick does not work for any potential, but only for potentials belonging to 
a certain class. In some cases, it works for potentials in the Majorana equation 
but not for the same potentials in the Dirac equation, which can be used to 
analyze the different character of Majorana and Dirac dynamics. In general, 
the transformation does not leave the probability density unchanged but this 
happens in some particular cases, showing us an additional amazing feature: 
under certain potentials the particle behaves as a free particle. We extend our 
results also to two-body equations. 

Although our method works in principle in any dimension and any represen- 
tation, we will focus in this work in the 1-D case and the particular representa- 
tions commonly employed in the experiments. 

6.3.2 One particle systems 

Let us come to explain our results with more detail. We will consider the 
following 1-D Dirac equation in natural units: 

ilp — —i cr x ip' + (er 2 m + V(x)) ip (6.64) 
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where ', ' denote time and space partial derivatives respectively. In the same 
representation, Majorana equation looks like: 



iip — —ia x ip' — ioy mip* + V(x)tp 
In 1-D a general potential can be written down as j!38j : 

V = fi(x) + f 2 {x) cr z + f 3 (x) a y + f 4 (x) cr x 



(6.65) 



(6.66) 



We start from the Eq.(6.65) with the potential in Eq.(6.66). We will analyze 
under what conditions a transformation of the form 

ip = JJe _,e ^ a! V ( 6 - 67 ) 

3 

where 



j~J e -iej(a;) _ £ ~i Fi^a* e -i F 2 (x)<j y e -i F 3 {x)<Jz e -i F 4 (x) 



(6.68) 



can convert Eq. (6.65) in the corresponding free Majorana equation for cf>. 
(Please recall that e A e B ^ e A+B unless [A, B] = 0, which in general is not the 
case here.) To this end we first notice that after applying Eq.(6.67) the LHS of 
Eq.(6.65) becomes: 

Y[e~ i& ^ x hi> (6.69) 

3 

while due to the anticommutation properties of the Pauli matrices the first term 
of the RHS transforms as: 



n e ~ iAj {x) h + (II e ~ 4Aj (x) y(-^ 



with 



II' 



-i Fi (x)a x i F 2 (x)a y i F 3 (x)a z —i F 4 (x) 



(6.70) 

(6.71) 



If we choose the F's such that: 



F[{x) = h(x),F^x) = -if 2 (x), 

F^{x) = if 3 (x),Fi(x) = U(x). (6.72) 

the second term of Eq.( 6.70| ) cancels out the one coming from the third term of 
Eq.(6.65), that is, the potential is removed. Finally, if Fx is real and F 2 , F3, F4 
imaginary we have: 



■or 



iQj(x)\* 



r = (II' 



--iAj (a;) 



)a y . 



(6.73) 



Putting all together, we have that applying Eq.(6.67) in Eq.(6.65) we obtain the 
following equation for <$>, provided that F\ is real and the other F's imaginary, 
and that they are related with the f's in the potential through Eq.(6.72|: 



n« 



Y[e- iA iW(ia x ct>' -ia y mj>*). 



(6.74) 
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Then, if F 2 (x) = -^M^O — 0, 4> verifies a free Majorana equation. That is, if 
we want to simulate the dynamics of a particle in the state tp{x,t) under the 
1-D Majorana equation with a potential of the form V(x) — fx(x) + f4,(x)a x , 
we only have to study the state <p(x,t) which is related to it through i/)(x,t) — 
e -iF 1 {x)(T Xe -iF4,{x) a .^ ( j ) ^ x ^_ Since, Fx is real and F 4 imaginary, we have the 
following relationship between the probability densities: 

\^{x,t)\ 2 = <r 2iFi{x) \<t>M\ 2 - (6-75) 

So the probability density observed for cf> can be easily related with the simulated 
probability density for tp. Also, analogous relations can be derived, for instance 
for expectation values of observables. Of particular interest is the case in which 
F A (x) — (then f±{x) = 0) and the probability density is the same for ip and 
<f). Therefore, the probability density of a state under the potential is always 
the same as the one of a free state. We shall illustrate this with an example 
below. The method does not work in general for F 2 (x) ^ 0, F 3 (x) =^ unless 
approximatively in the regions of space in which F 2 (x) ~ 0, F^x) ~ 0, which 
does not necessarily entails f 2 (x) — 0, f${x) ~ 0. 



We now analyze the case of the 1-D Dirac equation Eq.(6.64|. Using the 
same techniques as in the Majorana case, if we consider a potential with fi(x) = 
/s(x) = 0: 

V(x) = f 2 {x)a z + U(x)a x , (6.76) 

the transformation 

1 p = e -iF2(x)a ye -iF 4 (x) ( j ) (g 77 ) 



transforms Eq.(6.64) into: 



g-'^W'sg-'ftW^ _ £ i F 2 {x)<j y e -i F 4 (x) 

(—ia x (f>' + a z m4>), (6.78) 

and similar conclusions as in the Majorana case are reached for this class of 
potentials. But the absence of the complex conjugation in the mass term, pre- 
vents the possibility of extract a phase in the RHS if we include Fx or F 3 , so 
potentials with fx and/or f% cannot be simulated, not even approximatively, 
with free Dirac Hamiltonians. 

6.3.3 Examples 

Now we will illustrate all the above with some examples. First, we will consider 
the 1-D Dirac and Majorana equations with a linear potential 

V(x)=gx, (6.79) 

and the transformation: 

ip(x,t) = e "' 8 2 2 ^ <f>(x,t). (6.80) 

In that case <f> satisfies the free Majorana equation and \ip\ 2 — \<fi\ 2 , so any 
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Figure 6.4: a) Evolution of a fermionic wavepacket under Eq.(6.65) with the 



potential of Eq.(6.79). b) Evolution of the corresponding state with the trans- 
formation in Eq.(6.80) under a free Majorana equation. In both cases: m = 0.5 
and (7 = 1. 



solution of the Majorana equation with potential has the same probability den- 
sity of a solution of the free Majorana equation. The same transformation does 
not work in the Dirac case, shedding light to the different behavior of Majo- 
rana and Dirac particles against such potentials, as can be seen in Fig. |6.4| 
Notice that the Majorana equation does not have stationary solutions and that 
an initial global phase is not conserved as a global phase during the evolution, 
as commented in section 6.1. Thus Eq.(6.80) cannot be written in terms of the 
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evolution of some tp(0) and 0(0). In other words, a solution (f> of the equation 
with potential is always equivalent to some solution tp of the free equation, but 
in general this tf> is a different state at each instant of time, not the evolution of 
one state ip(0). In Fig 6.5 we show that our method works as a good approx- 
imation in a certain region of space in the case of the Dirac equation with a 
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potential: 

V(x)=gxa z , (6.81) 
and the corresponding transformation: 

ifi(x,t) = e 9 2 ° v <j>(x,t). (6.82) 
6.3.4 Bipartite systems 

All the above can be extended to two-particle systems. For instance, it has been 
shown that a Lorentz-invariant two-body Dirac equation with an oscillator-like 
mutual interaction can be written in the center of mass reference frame |141j . 
In 1-D and the particular representation we are using, we have: 

iip — -= {ax — ot2)(ip' + mu}xj3i 2 ip) + (0% + p-zjmij) (6.83) 

v 2 

with ai = a x ® 1, a 2 = 1 ® cr x , fii = a z ® 1, 2 1 <£> <J Z , @ii — a y ® a y , x — 
(x\ — X2)/V%- Thus the corresponding two-body Majorana oscillator equation 
is : 

iip = — (ai - a 2 )(ip + mujx/3i2i>) ~ + ^mip* (6.84) 

v2 

with 0i = <j y ® 1 and 2 — 1 <£) a y . With the techniques explained above, we 
find that the transformation 

ip(x,t) = e~ m ^ 2p " (f>{x,t). (6.85) 



transforms Eq.(6.84) into a quasi-free two-body Majorana equation: 



e 01 *ij> = e m ""' ^ ( — 5= ("i " "2)^' + (/9i + 02)m<t>*) (6.86) 

v2 

Thus, in that case the situation would be similar to the one-particle example of 



Fig 6.5 Eq.(6.84) can be simulated with good approximation with a two-body 
free Majorana equation in a large region of space. Interestingly, the trick does 
not work for the two-body Dirac oscillator. 

6.3.5 Conclusions 

We have shown that Majorana and Dirac potentials can be simulated with free 
Hamiltonians, since there is a mapping between free states and states under 
the action of the potential. The method works better in the Majorana case, in 
which there are potentials that cannot be simulated with a free Dirac Hamilto- 
nian, as has been illustrated with a particular case. This example also exhibits 
the peculiarity of conservation of probability density between free and non-free 
states, illuminating the issue of the different behavior of Majorana and Dirac 
particles against such a potential. In other cases the method works only as a 
good approximation in certain regions of space. We have extended our results 
to two-particle interacting systems. 
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Appendix A 



In this appendix, we will give details on the computations of the quantities of 
interest in Chapter 3. We first start with A and X in (3.9 1. Both are a sum of 
second order transition amplitudes, which can be written as: 

(^) 2 y2(f\ H i\k)(kWl\i) f dtt V dt 2 e l(E ^ E ^ /h e l ^~ E ^/ h (A.l) 
k Jo Jo 

being Ef, E k and E{ the energies of the final | /), intermediate | k) and initial 
| i) states of the system, respectively. The sum over k is a sum over all the 
possible intermediate states of the system which in the case of fixed two-level 
atoms reduces to a sum over all the momenta and polarizations of the emitted 



photon. The time integrations in (A.l) are just 

p i(E k -E,)t/h _ I 



e i(E f -E t )t/h _ 1 

(E f - E k )(E f - Ei) ~ (E f - E k )(E k - 



(A.2) 



The second term in (A.2) is usually neglected, but give rise to a very different 
short time behavior 53J (oc t 4 , not cx t 2 ). Therefore, it is of interest for our 
purposes. 

In order to obtain X we have to sum over the amplitudes for single photon 
emission at atom A (B) followed by absorption at atom B (A). The case where 
a photon is emitted and absorbed by the same atom corresponds to A, that we 
will consider below. Using the mode expansion for the electric field: 



E(x)=J^3 5:/^(e*-e(k,A) 



,-tk 



V (k,A)4 A ), (A.3) 



(with [cifcA, a \'\' ] — <5 3 (k — \s!)8\\>) taking into account (A.l) and (A.2), re- 
calling that 

J2 4 (k, A) e 3 (k, A) = % - h kj , ( A.4) 

A 

and using the tabulated integrals that we list at the end of the appendix, a 
somewhat tedious although straightforward computation leads to: 



X = 



td % d? 



(-V^ii+ViVj)/ 



(A.5) 
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a being the fine structure constant and 



/ = /++/_, (A.6) 
which, in terms of the dimensionless parameters z = VLr/c and x = r/ct are 
1± ' 



u 



-^{e lz [Ei{-iz) - Ei(-iz(l ± l/xj)] + e~ iz [Ei(iz) - Ei(iz(l ± l/x))]\k.7) 



for x > 1, I having the extra term i7r(l — l/x)e~ lz for x < 1. We use the 
conventions of |54j . As noted in |53j . the non-zero contributions for x > 1 come 
from the second term of (A.2). We display here the results of the derivatives 
in (A. 5) only for the particular case where the dipoles are parallel along the z 
axis (d^ = ds = d = du z ) and the atoms are along the x axis, corresponding 
to the physical situation considered previously in, for instance, |10j and in this 
paper. Actually, \E) is a triply degenerate state \E,m) with m = 0, ±1 and 
our scheme holds for a transition with Am = |80j . Another independent 
possibility would be to consider transitions with Am = ±1 that corresponds to 
d = d (u x ± iu y )/V2 [SU]- We find that: 

x =-^4€ 



r{4z(-l + { ~ 2 tflT* ) + e*'[-2 x z 2 Ei(-iz) 



+ (2 + z (-2i + (-1 + x)z)) Eii- *- 1 ^" ) 
+ (-2 + + z + xz)) Ei(- l ^±^)} 



-2e 



Ei( 



z e~ lz [-2 x zSi(iz) + (2i + (-1 + .x)z) 
B»(ii=^) + (-2i + z + xz)E l (^±^)}} 



(A- 



for x > 1, with the additional term zae -iz <i 2 (2 + z (2i + (—1 + x)z)) /(r 2 x) for 
x < 1. 

Now we come to A, which is the sum of the radiative corrections of atoms 
A and B. As can be seen in the main text, A appears in our results only as a 
higher order correction to X. Therefore, instead of finding an exact expression 
for it, we are mainly concerned with removing the divergencies. We followed the 
standard treatment (see, for instance, [79]) which is valid for the times fit > 1 
we are considering. From (A.l), it is possible to arrive at: 



-2ia\ d\ 2 t 
3 7r e 2 c 2 



lim 



dojuj 3 ( 



1 



+ Jo 



n — lj + ie O + cj — ie 



Now, using in (A.9) the identities 

, ,3 



n 3 



(A.9) 



(A.10) 



the first term of (A. 10 1 cancels out the contribution of the Hamiltonian self- 
interaction terms Eq. (2.15), the second is the state-independent contribution 
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that can be absorbed in the definition of the zero of energy [73], the third 
cancels the counterterm coming from the mass renormalization 79] and finally 
the last term has logarithmic divergences and a cut-off, related with the fact 
that we are in the electric dipole representation could be imposed at t min = 
^r- = 1.76 ■ 10~ 19 s. Please notice that the times relevant in our computations 
are of the order of t = ^ ss 4 • 10~ 15 s. Therefore: 



A _ 2ia\d\ 2 z 3 k ( 
3 7r L 2 e 2 x 



1 



(A.ll) 



with z max /z = c/(ftoo). 

Another quantity of interest is what we called I in the main text and it is 
given by ( ]3.24 l and (3.25). Performing the integration in (3.251, we obtained 

M(z,x) = M + (z,x) + M-(z,x), where: 

M ± (z,x) 



/^{sin(z(l±-))[cz(z)-cz(z(l±-))] 

47T Z Z X X 

cos (z(l ± -))[si(z) - si(z(l ± -))]} 

X X 



(A.12) 



The derivatives in (3.24) were performed in the same particular situation as in 
X. 

\U\ 2 = (0|<Sj>S^|0}, and \V\ 2 = (0\SgS%\0), are just the two terms that 
contribute to A without the time ordering and therefore their divergencies are 
removed by the application of (A.9). Taking this into account, we obtain: 



\U\ 2 = 
\V\ 2 = 



2a\d\ 



2,2 



3tt e 2 L 2 
2a\d\ 2 z 2 
3ire 2 L 2 



(-2- 
(2- 



•7T- + 2cos(-) + 2(-)si(-)) 

X XXX 



tt- - 2cos(-) - 2(-)si(-)) 

X XXX 



(A.13) 



F and G can be written in terms of previously computed quantities, taking into 
account that: 



F = 9{h - t 2 ) ( u A (ti)v' A (t 2 ) + v' A (h)u A (t a ) 

+ u B {h)v' B (i a ) + u' B {ti)v B {t 2 ) ) 
G = u A v' B + u' A v B , 



(A.14) 



being V A = ( 1 1 \ ) , and U B = ( 1 1 S B \ ) . The primes are introduced to 
label two different single photons. Therefore, in the computation of | F | 2 , \G\ 2 
and F G* we will only need, besides | U | 2 , | V \ 2 and I, the following: 



V A U* A = V B U* B = 



2a\d\ 2 z 2 
3irL 2 



sin ■ 



V A V* B 



(-V 2 ^+V 4 V,-)/' 



(A.15) 
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being I' = I' + + I'_, with: 

I', = 1 -^-^{e lz Ei(-iz(l± -)) + e- tz Ei(iz(l± -))} 
2 x x 

- e Tlz Ei(±iz)} (A.16) 

and Ua = Va Vg when x > 1, with the additional term — 2-7rsinz (1 — 1/x) 
when x < 1. Again the derivatives were performed as in X and i. 

The following integrals are useful to obtain the results of this appendix [M] : 



DO 



(I 



duj- = -e^^ p Ei{± i^B) (A.17) 



f 

/ dco- = —e ±ij 13 (Ei(=fi~f (3) =p in), 

Jo uj - p 



with a > 0, arg ft < n. 
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Appendix B 



In this appendix we will give further details on the computations of the relevant 
magnitudes |[/a| 2 , |^b| 2 which are necessary to compute the concurrence 
Eq.lpOl in section 4.2. With Eqs. (pL26l- l429| and Eq. Ol: 



X = 



d 2 Nv 



ti 2 



dk\k\ (e lkr I t i+e 



ikr ■. 



h2) 



with 



hi,2 



dt 



2,1 



(B.l) 



(B.2) 



Notice that the term with I t 2 gives the non-RWA probability amplitude associ- 
ated to a single photon emission of qubit B followed by an absorption of qubit 
A. Performing the time integrations, inserting them in Eq. (B.l ) and after some 
algebra, X can be given as a combination of integrals of the form: 

cos( kj) 



dk 
dk 

dk 
dk 



k + p 
cos( &7) 
k - p 

sin( k 7) 
k + p 
sin( k 7) 

k - P 



= — sin(7/3)si(7/3) — cos (7/?) 0(7/?) 

= — sin(7/3)st(7^) — cos (7/3) ci (7/3) 

— 7rsin(7/3) 

= sin(7/3)Ci(7/3) — cos(7/3)si(7/3) 

= - sin(7/3)Cz(7/3) + cos (7/?) si (7/?) 

+ 7rcos(7/3) 



(B.3) 
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with 7, (3 > and the conventions in [S3] for the si and Ci. Putting all together 
we find: 

X = ^-(mp^sm(p) - e^((l + zr_)(C(r_) - S(r_) 

-7rsin(r_)e(l - 0) + (1 - ir+)(-C(r+) - S(r+) 
-7rsin(r + ))+)(-r_ + i)(-5C(r_) + CS(t-) + tt 
cos(t_)) + (-7+ - i)(-SC(T+) +C5(t+) +7rcos(r + )) 

-2) - e^((l - ir_)(-(7(r_) - S'(r-) -7rsin(r_) 

9(£ - 1)) + (1 + it + )(-C(t + ) - S(t+)) + (r_ + i) 

(5C(r_) - CS(t-)) + (t_ - i)(S , C(r+) - OS(t+) + tt 

cos(t+)) - 2) - 2 - 2C(p) - 25(p) - p{-2SC{p) 

+2CS(p)) (B.4) 



where £ has been defined in Eq. (4.32), 



p = Clr/v 



is a dimensionlcss distance, 



and we define 



C(x) = cos(x) Ci(x), S(x) = sin(x) si(x), 
CS{x) — cos(x) si[x), SC(x) = sin(x) Ci(x). 



(B.5) 



(B.6) 



(B.7) 



Notice the dependence with the spacetime region through the factors with the 
Heaviside function 0. 

Now we come to the emission probabilities |C/a| 2 , I^bI 2 , which are given by 

\U A \ 2 = {Q\S+S A |0) , \V B \ 2 = {Q\S B S+ |0) (B.8) 

Following similar techniques we find that 

\u A \ 2 = f+(nt),\v B \ 2 = f.(nt), (b.9) 

where: 

f ± = —(irnt±2(cos(nt) + ntSi{nt) - 1)) (b.io) 

where Si must not be mistaken by si, Si = si + ir/2 as usual. 

Finally, notice that A (Eq.(4.31)) is a sum of two terms like \Ua\ 2 and \Vb\ 2 
with the time ordering operator T and that Ua * V B (Eq.(4.36l) is similar to X 
without T. The treatment of the divergencies has followed the lines sketched in 
Appendix A. 



112 



Appendix C 



In this appendix the task is to compute the term SM ' introduced in section 
5.1. This term is included in the whole third order amplitude 

Mi3) = (^) 3 £</IWHm</3|^> x (c.i) 

a/3 

x f dtie iw f<* u f 1 dt 2 e lw ^ t2 f 2 dt 3 e lw ^ t3 
Jo Jo Jo 

where a, /3 stand for arbitrary sets of quantum numbers and the sum includes 
integrals over momenta and polarization sums. We start computing the time 
integrals: 



dhe™"" 11 f 1 dt 2 e m ^ t2 f 2 dt 3 e iw ^ 3 
( i Jo Jo 

I e iw fi t _ ^ e iw fa t _ i 



(iw/3i)(iw ai ) iwfi 



lWf a 



e iw f0 t 



1 e tWfa 



(iwpi)(iw a p) iwfp 



>Jf a 



(C.2) 



The last two lines in Eq.(C.2) contains all the time dependence along with that 



in the final and exchanged photon frequencies. Among all the terms that can 
interfere with Vb, we select only those that depend on the interatomic distance 
r whose contribution 5A4^ we will now compute. We first consider the three 
cases where the exchange "goes from A to B" , where we display the frequencies 



that enter in Eq.(C.l) and Eq.(C.2| in terms of 



s / = w f 



n, A f = Wf 



flandA 



w 



SI, 
- Si, 



(C.3) 



where Wf,w are the final and exchanged photon frequencies, and SI = (E e — 
E g )/h is two level frequency gap. The primed objects that will appear in the 
equations are associated to the energy E' g of the intermediate atomic ground 
states; we will take E' g — E g — » at the end of the calculations. 
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Wfa = Wf a = —A Wf a = — A 

W a p = -A W al 3 = S'/ = S' 

= A' u;^ = A' u>/jj = A^ 

w ai = A' - A = n' - = 6 

Figure C.l: The three diagrams contributing to in which qubit A emits 

The contribution to Ai^ 3 ' coming from Fig. 2 can be written as 

-^V^A) Ve 9 ^B)Vl ,MA){ /KI{v l i _ y f) {S{H f ,t) 

-Sp f ,t)) + ^(Sm - A,t) - S(Tf f ,t))} - 

l y/ eg (x 5 )K t 9 (xA)^ e (x^){ A/( ^ 1 +A) (5(s] /; t) - 

S(-A,t)) - ^r(S(£/ - A,t) - 5(-A,*))} - 
l K 8 (x B )^ (x A )4 e ( XA ){ * + A) (5( S/ , t) - 
S(-A,i)) - ^O^' - A*) -5(-A,t))} (C.4) 

In the above formula P(x) • • • stands for the contraction of the operators 
P and Q, namely 

Vfi(x)---V} 'Ay) = J dkNd fl ■ ■■d fi ,e ik ^ (C.5) 
and the function S is defined as 

izt _ i 

S(z,t) = (C.6) 
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Finally, we recall some issues about the contractions. First, the frequency that 



appears in Eq.( C.5 ) and in the definitions for A and S is w — v k. Second, we can 



pull out the lone operator V* from the contractions. Third, we will consider real 
dipol e moments, so that V eg — V ge , etc. This allows to factorize from equation 



(C.4) a common factor -(l/h 3 )V^ (x A )V(x B )V^ (x A )I^(w ,t), where lf{w,t) 



contains all the terms within braces. 






Wa/3 = -S' 
W Bi = S 



W fa 

w 8i = £ 



A' 



A', 



W a /3 



-A' 



A', 



Figure C.2: The three diagrams contributing to SAi^ in which qubit A absorbs 

The cases where the exchange "goes from B to A" are depicted in Fig. 3 
below. From arguments similar to those of the previous case, we can factorize 
out the V terms obtaining — (l/h 3 )V^ (x A )V(x A )V 1< (xb)Ib{ w 7 t), where 

IbM = S ( S /_^) ( 5 ( S /' f )- 5 ( S /' f )) + 

— (SCEff -£,*)- SCE' f ,t)) + , 1 a / \ 

(S(Z f ,t) - S(-A',t)) ^r(Sm - A',t) - 

5(-A^ ^ )) + ^ (5] L^ y ( < S(S / , ^ ) 4 S(-A^^))- 

-~(S(Z-A>,t)-S(-A',t)) (C.7) 

Finally, the contribution of third order to the amplitude we are interested in is 
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given by 

SM {3) = -{l/tf){V\x A ) V{x B )V \x A )I*{w,t) + 

V\x A ) V{x A )V Hx B )I*{w,t)} (C.8) 
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Appendix D 



Majorana equation in 1+1 
dimensions 



We start from the Dirac equation in 1+1 dimensions and in covariant form |20j : 

- ij° d Q ip - ij 1 dtip + mifj = (D.l) 

with h = c = 1 and the gamma matrices 7 M , fi = 0, 1 verifying the anticommu- 
tation rules of a Clifford algebra 

{7^,7"} = 2 <T, (D.2) 

being the 1+1 Minkowski metric. Multiplying by 7 and with the usual 
notation 

/? = 7 V = 7 V, (D.3) 
we can write it in the Schrodinger-like form: 

idtip — a pip + mptp (D.4) 

If we want a and /3 to be involutions (a 2 = 1 = 1) they must be Pauli matrices. 
In table |D.1| we spell the six possible cases. 

The M ajorana equa tion is obtained by replacing the \& in the mass term of 
Eq. (D.l I or Eq. (D.4| by the charge conjugate ^ c . Then: 

id t ip = apt/j + m(3ip c (D.5) 

The charge conjugation is defined in the following way |128j : 

= v CiJj* (D.6) 

where the matrix C must verify 

Cl^C- 1 = - 7 „ (D.7) 
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7 1 


a 


H 


Cases 


0~x 




-<?z 


~°zP 


— ima y K 


La 


ia z 


Uy 


(JyP- 


- ima y K 


Lb 


Oy 


ia x 


&Z 




- ima y K 


ILa 


ia z 


~0~x 


-<JxV 


+ imayK 


Il.b 




ia x 


-Oy 


-a y p 


— ima y K 


IILa 


io y 


O x 


&xP - 


- imayK 


IILb 



Table D.l: The six possible representations of the 1+1 Clifford algebra, with 
the corresponding Majorana Hamiltonian associated. 



and 77 is an arbitrary phase. Now, defining 

C = C0 (D.8) 

and using that 

hP = ll (D.9) 
if P is real or purely imaginary, we have: 

CllC = -7 M (D.10) 

and then C must commute with a y and anticonmute with u x ,z, so 

C = a y (D.ll) 

for all representations, and 

C = <7 y /3. (D.12) 
Therefore, the Majorana equation can be written in the following form: 

idtip = apijj + rjm/3ifj* (D.13) 

where (3 is given by 

P = PC = ±o y (D.14) 

and the plus sign corresponds to the representations in which P = a y (Cases II 
in table D.l ) and the minus to the ones with P = a XiZ (Cases I and III). Finally 



we define the antiunitary involution K [129 , given by 

ip*=Kip (D.15) 

and set the global phase to 77 = i by convenience. Thus the Majorana equation 
is 

idtip = apip ± imo~ y Ktp (D.16) 

where now the plus corresponds to II. In chapter 6 we choose the particular 
representation III. b and then 

idtip — o~x pip ~ imoy Kip (D.17) 
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This last equation can alternatively be thought ol as the dimensional reduction 
of one of the equations in which the 3+1 Majorana equation in Weyl represen- 
tation can be decomposed. 

Eq.(D.16) is a Schrodinger-like equation in which 

H = ap± imOyK (D.18) 

plays the role of a Hamiltonian. Hitherto we have considered the "free" Majo- 
rana equation, but we could add a potential V to Eq. (D.16I, obtaining: 

H ' = V + ap ± ima y K (D.19) 
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Appendix E 

Implementation of the 
quantum simulation of the 
Majorana equation in 
trapped ions 



We can simulate Eq. (6.48) using two trapped ions that are subject to two 
dynamical terms, coupling both to the internal states and motion of the ions. 
The kinetic part, cp x (l £§> a x ), is created with a laser tuned to both the blue and 
the red motional sideband of an electronic transition 135, 138J, and focussed 
on ion 2. The spin-spin interaction term, a x ® a y , is derived from detuned red 
and blue sideband excitations acting on each ion in addition to the previous 
one |1421 1143] . The Hamiltonian describing this situation reads 



H 



= h ^ 




T ^2 + 


hua'a + hv r b^b 




+ Ml 


(e i( - qZl ~ 






+ H.c 


+ m 


(e l(qZ2 - 


-U>2t+<f>2) 


_(_ e i(qz2-u'- 2 t+<t>' 2 )^ (J + 


+ H.c 


+ Ml 




-ujt+<t>) _ 


_ e i(qz 2 -ui't+(t>')^ cr + _|_ 


H.c. 



(E.l) 

Here z± t 2 = Z ± | are the ion positions, measured from the center of mass, Z, 
and relative coordinate, z. The phases of the lasers fa for i = 1,2, (fa fa), arc 
controlled to perform the interaction term (kinetic term). The frequencies of 
the center of mass and stretch mode are given by v and v r = s/Zv, while a\ a, 
, and b, are the corresponding creation and annihilation operators. Finally, Q 
and £1 are the Rabi frequencies of the lasers within the bounds of applicability 
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of the rotating-wave approximation. With an appropriate choice of parameters 



ujq + v r — S 
= UJq — v r + S 
L02 = — v T + S 
uj' 2 = uj + u r - 6, 



ui = 

L0> = 



ljj a + V 

0. 



the Hamitonian (E.ll in the interaction picture reads 



H 



hrj r Vl{a x 



iSt 



be 



= tt/2 

= tt/2 

= 

= 0, 



(E.2) 



-iSt\ 



+hr]fl(l (g) a x )i(a i — a) 



(E.3) 



where r\ = ry r 3 1//4 = yJh/Am'v <C 1 is the Lamb-Dicke parameter and m' the 
ion mass. In the limit of large detuning, we have 



S > 77 r fM(&t&), 7y r>|(a T - a )\. 



(E.4) 



We recover Eq. (6.48) with the momentum operator p x = ih{a) — a)/2A and the 
equivalences 

c = 2,Af2, - 



mc = 



(E.5) 



4 m 



is the size of the harmonic oscillator ground state. Note 



Here, A = 

that it is possible to explore all velocity regimes. Introducing the ratio 7 
\mc 2 /(cp x )\, 

2{r lr n/5) 2 



7 



(E.6) 



Ki(ot-o))|(rjn/<y)' 

it is possible to tune the numerator and denominator independently so as to 
preserve the dispersive regime, while exploring simultaneously the range from 
7 ~ (ultrarelativistic limit) to 7 —¥ 00 (nonrelativistic limit). 

We could also consider other implementations that do not require synchro- 
nization of laser phases for different beams. Consider, for example, the Hamil- 
tonian 

H = cl (8) (p x a z ) — mc 2 a y <g) a y , (E.7) 



which is equivalent to that of equation (6.48) up to local unitary rotations. 
Using a detuning of ±^/2 for the blue and red sideband, respectively, in the laser 
focussed onto ion 2, leads to an interaction of the form 8A?yf2 2 l <g) (p x a z )/v and 
the equivalence c = 8A7]fl 2 /u. The spin-spin interaction in the second term 
can be implemented by an additional global bichromatic light field acting on 
the red and blue sidebands of the stretch mode, similar to the case above. In 
the present case, however, the ions experience laser light with the same phase 
4>i = 4>2 = 4>'i — 4>2 — an d a single laser beam can be used. 



Measurement of the pseudo-helicity 
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In experiments using trapped ions, the only observable that can be di- 
rectly measured by fluorescence detection for each ion is a z , and additional 
laser pulses are needed to map other observables onto it. The application of a 
state-dependent displacement operation on ion 2, J7 2 = exp(—ik(l<8><j y )<8)p x /2), 
generated by a resonant blue and red sideband, followed by a measurement of 
1®<j z is equivalent to measuring the observable 



A(k) = U%(l®(T g )U 2 

— cos(kp x )(l ® a z 



sin(fcp a; )(l ® a x ). 



(E. 



Here, k is proporti onal t o the probe time t pro b e 22\. In order to measure 
the first term in Eq. (6.61), we note that -gfe{A(k))\ k _ oc ((1 (g) a x ) ®p x ). 
Therefore, this term can be measured by applying a short probe pulse to the 
ions and measuring the initial slope of the observable A(k) [1361 122] . 

To measure the second term in Eq. (6.61 1, we have to apply a different state- 
dependent displacement operation to ion 1: U\ = exp(—ik(a x (8> 1) <g>p x /2), and 
measure the spin correlation <j z ®<j x , which requires an additional 7r/2 pulse on 
ion 2. Again taking the initial slope in this observable, we have 



d(a z 



dk 



2((a v ® a x ) 



(E.9) 



fc=0 
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Conclusions 



Throughout this Thesis, the specific results have been summarized at the end of 
each section. Now we will outline the most important results and conclusions. 

Dynamics of entanglement 

We have characterized the dynamics of the generation and destruction of entan- 
glement between qubits mediated by a quantum field, focusing in the relation- 
ship with the spacetime region at which the qubits are placed. More precisely, 
we have considered a pair of two-level atoms -real or artificial, that is, super- 
conducting qubits- A and B separated by a fixed distance r interacting through 
a quantum electromagnetic field. 

• (Sections 3.1, 3.2 and 4.2) If the qubits are initially in a separable state, 
two different spacetime regions emerge from our computations. If v t < r 
- t being the time of the interaction and v the propagation velocity of the 
field quanta- entanglement can be generated between the atoms due to the 
non-locality of the Feynman propagator - or in other language, "vacuum 
entanglement" - which give rise to correlations at any time. In the weak 
coupling regime, which is the case with real atoms and the electromagnetic 
field, we show that these correlations are only classical, but entanglement 
-even maximal- can be generated through a measurement of the state of 
the field. For stronger couplings, as is the case in 1-D circuit QED, a small 
amount of entanglement shows up even in the absence of a measurement. 
If v t > r, a very different behavior of entanglement appears due to photon 
exchange, making apparent that entanglement generation may be a good 
signature of single photon propagation, a fact of particular significance in 
the framework of circuit QED. In all the cases the point vt = r plays the 
role of a frontier between two different spacetime regions regarding the 
behavior of entanglement dynamics. We have devised a detailed circuit 
QED proposal of an experiment to test all the above physics. 

• (Section 3.3) If the qubits are initially entangled, we have characterized the 
remarkable phenomenon that a full Entanglement Sudden Death-Sudden 
Birth cycle can take place for vt < r. In this case, these entanglement 
destructions and revivals are related with the revivals and destructions of 



123 



the atom- field entanglement, that is to say, there is sort of an entanglement 
flux between the different pairs of the tripartite system atoms-field. That 
this can happen at times at which photon exchange is not allowed represent 
an additional insight on the nature of the generation and destruction of 
quantum correlations. 

Non-RWA effects 

We have seen that the reasons for going beyond the ubiquitous and celebrated 
Rotating Wave approximation are twofold. On the fundamental side, it is nec- 
essary to include all the non-RWA terms for a complete theoretical short-time 
analysis of the matter-radiation interaction or to properly deal with causal- 
ity questions, although in these systems these effects are beyond experimental 
reach. On a more applied viewpoint, circuit QED in the ultrastrong coupling 
provides a framework in which the non-RWA effects are available to experiment. 
We have explored a striking consequence of working beyond RWA, the fact that 
a qubit in the ground state can get excited and emit a photon, even if the field 
is also in the vacuum state. 

• (Section 4.1) In the simplified model of a circuit QED analog of a cavity 
QED system, we have introduced the "slow quantum anti-Zeno effect", 
providing a full experimental proposal to detect a qubit in its excited 
state with certainty after a few measurements starting from the ground 
state of the whole system. The protocol is based in the certain probability 
of excitation already present in the system for strong enough couplings 
and in the non-equilibrium dynamics after a single measurement. This 
feasible and realistic experiment would shed light onto the physical reality 
of a phenomenon commonly considered as "virtual" . 

• (Section 5.2) In the case of the full multimodc Hamiltonian of circuit QED, 
we have analyzed the implications of the mentioned ground state qubit self 
excitations and in general the non-RWA dynamics in the interpretation 
of the measurements of the probability of excitation. In particular, we 
have seen that only for long interaction times we can say with certainty 
that a qubit detector's click entails the decay of an excited source. For 
short -although experimentally accessible in circuit QED- times, non-RWA 
contributions dominate in the dynamics of the probability of excitation 
conditioned to the decay of the source, meaning that in this regime a click 
is more likely linked with a self-excitation. 

Causality 

We have revisited the old question -Fermi problem- on causality in matter- 
radiation interactions, which we can rephrase as follows: "Is the probability of 
excitation of the atom B causal?" To properly deal with this issue, we have first 
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to pose a different question: "What does causaiity mean?" As we have seen, 
quantum correiations may show up at superluminal rates, thus the probability 
of a given nonlocal state such as one atom in the excited state and the other in 
the ground state becomes non-zero at v t < r. Also we have to take into account 
the mentioned fact that an atom can get self-excited at any time. But all these 
is compatible with causality since cannot be used to transmit information faster 
than light. The relevant question is actually "Is the probability of excitation of 
atom B independent of atom A at times t < r/vV 

• (Section 5.1) We have shown that the answer to the latter question is 
affirmative, by providing a complete non-perturbative proof of the inde- 
pendence of the probability of excitation of atom B with regard to atom 
A. 

• (Section 5.1) We have exploited again the experimental amenability of 
circuit QED to propose a realistic experiment to check the causal behav- 
ior of the probability of excitation, coming true a hitherto considered as 
"gedanken" experiment. The experiment, together with the theoretical 
clarification given above, would close a long-lasting controversy on the 
Fermi problem, confirming the causal behavior of Nature at this level and 
shedding light on the meaning and implications of the notion of causality. 

Quantum simulations 

Although circuit QED can be understood as a closed ID universe and its physics 
is interesting by itself, it can also be understood as a quantum simulator for 
matter-radiation interactions. Besides, we have also introduced some results on 
quantum simulations of rclativistic quantum mechanics with trapped ions. 

• (Sections 6.1 and 6.2) We have introduced a new family of generalized 
hamiltonians, 'Majorana Hamiltonians" which do not follow the standard 
textbook definition of hermiticity but we have shown that induce a norm- 
conserving dynamics and can always be simulated with regular Hamilto- 
nians. In particular, we have proposed a trapped ion quantum simulation 
of one conspicuous member of the family, namely the pseudo-hamiltonian 
appearing in the Majorana equation- the Relativistic Quantum Mechan- 
ics equation of a fermion with a Majorana mass term-. We have explored 
astounding properties of the Majorana Hamiltonians, as the possibility 
of measure initial global phases in the evolution of physical observables. 
Thus, this experiment would be useful to explore fundamental questions 
being at the edges of Quantum Mechanics and Quantum Field Theory. 
But the scope is more general and includes possible simulations of oper- 
ations like time reversal, charge conjugation or partial transpose in the 
lab. 

• (Section 6.3) We have introduced a method to use the current simulations 
of free relativistic Hamiltonians to simulate also Hamiltonians including 
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potentials, opening the way to more complex simulations of interacting 
many-body relativistic equations, in which the results of the first parts of 
this Thesis could be explored within different frameworks. 
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